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An inquiry is made into the nature of the corpuscular entities which possess the properties 
necessary to explain the soft component of the cosmic radiation on the basis of the recent theory 
of W. F. G. Swann. The conclusion is reached that the most likely entities are protons which lose 
energy according to the relation —dE/dx=\E+a, where a represents the energy lost per unit 
path by ionization, and XE the energy which goes into the production of secondaries in the 


amount that is actually observed in the form of showers. 


INCE the discovery of the latitude effect 
which demonstrated that there must be a 
corpuscular component of the cosmic radiation, 
there have been many attempts to find out what 
the nature of the corpuscular entities must be. 
The most complete treatment of this question 
has been given by A. H. Compton.' Compton has 
discussed the problem from the point of view 
that, to a large extent, the observed cosmic-ray 
particles are the primary entities themselves, and 
he has concluded that, at sea level, the hard com- 
ponent consists most probably of protons while 
the soft component may be ascribed to electrons. 
W. F. G. Swann? has discussed corpuscular 
theories from another point of view, namely, that 
the observed cosmic-ray particles are largely 
secondary in character, and that the number of 
secondary particles accompanying a primary one 
varies with the energy of the primary ray. If this 
view is adopted, then the considerations regard- 
ing the nature of the primary entities must be 
modified. It is the purpose of this paper to in- 


1A. H. Compton, Proc. Phys. Soc. 47, 747 (1935); Rev. 
Sci. Inst. 7, 71 (1936). 

?W. F. G. Swann, Phys. Rev. 47, 575 (1935); 48, 641 
(1935) ; 50, 1103 (1936). 


quire into the nature of the primary entities from 
the point of view that the observed cosmic rays 
are largely secondary. 

W. F. G. Swann has shown that the experi- 
mental data are well satisfied if it be assumed 
that the soft component of the cosmic radiation 
is corpuscular in nature and charged, and loses 
energy according tothe relation —dE/dx=E+a, 
where a represents the energy lost in the excita- 
tion and ionization of the atoms per unit of path, 
and XE is the energy which is lost by the produc- 
tion of secondaries. One of the numerical values of 
\ which results from these considerations is 0.3 
per meter of water. It is our purpose to discover 
what the nature of the primary corpuscles must 
be in order that they should lose their energy 
according to this law with this value of \. The 
only candidates which we need consider are pro- 
tons and electrons, since particles of greater mass 
or larger charge will not be able to penetrate to 
any great depths into the earth’s atmosphere.' 
The energy loss a is practically the same for both 
electrons and protons of high energy, and we 
cannot utilize its numerical value as a clue to the 
identity of the primary rays. Besides this con- 
stant energy loss, charged particles may lose 
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energy by three processes: (a) by radiation, 
(b) by the production of pairs of electrons, and 
(c) by the production of showers. The cross 
sections for energy loss by processes (a) and (b) 
have been calculated by Bethe and Heitler,’ and 
Nordheim.‘ The cross section for process (b) is 
only about 1/137 that for process (a). Although a 
theoretical ‘mechanism for the production of 
showers has recently been given by Heisenberg,® 
and Oppenheimer,® the calculations are not as 
yet sufficiently detailed to give a value for the 
energy lost in this way. However, we may utilize 
the observations on the frequency of occurrence 
of showers to derive a value for this cross section 
in the following manner. 

Let F(E)dE be the number of primary rays of 
energy E crossing a unit area at sea level per 
second. The total amount of energy going into 
secondaries of these rays per unit path length is 
then /AEF(E)dE, where the integral is taken 
over all values of the energy. Let e be the average 
energy of a shower ray, and R(N) be the fre- 
quency distribution of showers of N rays. For 
R(N) we may take the experimentally determined 
form R(N) = 7Z*/N°,’ and hence 

ve 
prer@nae- | eNR(N)dN=eyZ?/2. 
“s 
Now the total flux of primary energy through our 
unit area is /EF(E)dE. This energy will even- 
tually all be dissipated in the excitation and 
ionization of atoms below sea level. If J is the 
average energy lost in producing one ion pair, 
and J(x) the observed rate of ionization per unit 
volume caused by cosmic radiation at a distance 
x below sea level, then, neglecting the small 
amount of ionization produced below our unit 
area by secondary rays which originate above 
this area, the flux of primary energy is given by 


f EF(E)dE=J | I(x)dx. 


“0 


Thus we can find \ in terms of quantities which 


3H. Bethe and W. Heitler, Proc. Roy. Soc. A146, 83 
(1934). 

‘L. Nordheim, J. de phys. et rad. 6, 135 (1935). 

5 W. Heisenberg, Zeits. f. Physik 101, 533 (1936). 

6 J. R. Oppenheimer, paper delivered at the Harvard 
Tercentary Conference (1936). 

7™C.G. Montgomery and D. D. Montgomery, Phys. Rev. 
50, 490 (1936), and earlier references there given. 


may be measured experimentally: 


3a 


A=eyZ?/2J | I(x)dx. 
“0 


A value for y may be estimated from experi- 
ments by counters and by ionization methods. 
The counter observations are uncertain, on one 
hand, because of the difficulties in estimating the 
efficiency of a particular counter arrangement, 
and, on the other hand, the evaluation of the 
ionization experiments necessitates an uncertain 
extrapolation to the showers of smallest size. 
However, the values of y obtained from the work 
of various observers’:* vary from 0.7X10-% 
sec.-! mole! to 5X10-* sec.-' mole, and we 
may take for these purposes the value of 2 10-° 
sec.—! mole. J(x) may be taken from the obser- 
vations of R. A. Millikan and his collaborators,’ 
and we find, for the ionization produced by the 
soft component below sea level, the value of 
7.1X10* ion pairs cm~ sec.~'. Assuming J as 
32.2 ev, Z?=52 for air, e=10* ev, we obtain 
\=0.16 per meter of water. 

The values of \ for the processes (a) and (b) 
may be calculated from the known expressions. 
Table I contains the values of \ for these methods 
of energy loss for protons and electrons in units 
of meter of water. The value of \ which is to be 
compared to the 0.3 m~ deduced from the theory 
of W. F. G. Swann is the sum of the }’s for all 
processes of energy loss with the exception of the 
energy lost by ionization and excitation of atoms 
along the path. Thus, even if we postulate that 
electrons are unable to produce showers, the 
total, theoretical \ for electrons is too large. On 
the other hand, if we assume that showers are not 
produced by protons, the theoretical value of \ 
is much too small. Thus it appears that the soft 
component of the cosmic radiation consisis of 


TABLE I. Values of \ in meter of water for energy lost by 
various processes. 





Electrons Protons 


Energy loss by 
Radiation A 2.2/(1840)? 
Pair Production 0.0073 0.0073 /(1840)? 


Showers 0.16 








8 H. Geiger and O. Zeiller, Zeits. f. Physik 97, 300 (1935): 
J. C. Street and R. T. Young, Phys. Pev. 47, 572 (1935): 

*1. S. Bowen, R. A. Millikan and H. V. Neher, Phys: 
Rev. 44, 246 (1933). 
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primary protons which lose their energy according 
to the law —dE/dx=\E+a, where XE represents 
the energy loss by the production of showers. The 
value of 0.16 for \ derived from shower data 
differs from the 0.3 required by Swann’s theory 
only by a factor of two, well within the uncer- 
tainty of the data involved in calculating it. For 
example, if y is taken as 4 10~ sec. mole, or 
eas 2X10 ev, values well within the range of the 
observations, the \ calculated from shower pro- 
duction will agree exactly with that required by 
theory. 

This hypothesis can explain, naturally, an 
east-west effect due to positively charged par- 
ticles, which increases with altitude as does the 
soft component, although the detailed explana- 
tion of this effect is much more complicated than 
is suggested here. The hypothesis is not in con- 
tradiction with the recent experiments attempt- 
ing to detect a proton component of the cosmic 
radiation,'® since the actual number of protons 
passing through an area at sea level will be small 
compared to the number of secondary rays which 
accompany them. Thus, although the earth’s 
magnetic field will affect the soft component of 
the radiation before it enters the atmosphere, and 
so produce the.latitude and azimuthal effects, the 
cosmic-ray particles which will be observed will 
be, for the most part, positive and negative 
electrons: the shower particles which the primary 
protons have produced. 

The recent theory of shower production of 
Heisenberg’ results in the conclusion that about 
half of the energy lost by a primary ray in pro- 
ducing a shower will go into neutrinos, and it 
would seem necessary to take this into account 
by multiplying e, the energy of a shower ray, by 





°C. G. Montgomery, D. D. Montgomery, W. E. 
Ramsey and W. F. G. Swann, Phys. Rev. 50, 403 (1936) ; 
R. B. Brode, H. G. MacPherson, and M. A. Starr, Phys. 
Rev. 50, 581 (1936); C. G. Montgomery and D. D. 
Montgomery, Phys. Rev. 50, 975 (1936). 
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a factor of two. However, this is not necessary, 
since again half, approximately, of the energy 
which passes through a unit area will also go into 
the production of neutrinos and not into the 
eventual production of ionization. Thus to the 
approximation considered here, the presence or 
absence of neutrinos does not influence the result. 

The validity of the above choice of protons for 
the soft component primaries rests, of course, 
upon the assumption that Bethe and Heitler’s 
theoretical formula for the radiation of fast 
electrons is correct for the energy range here 
involved. In the past, considerable doubt has 
been cast upon it, and several modifications" of 
the expression have been suggested. It is not our 
purpose to discuss these here except to point out 
that the recent results of Anderson and Nedder- 
meyer” restore some confidence in the correctness 
of the formula. If the formula is proved to be 
invalid, however, the suitability of protons for 
the soft component is not affected, but we may 
also have to admit the possibility of electrons; 
electrons, however, which do not lose energy pro- 
portionally to their energy, and the corpuscular 
theory would have to be generalized to take this 
into account. 

Thus, with due regard for the possible objec- 
tions outlined above, the most probable entities 
which behave in a manner necessary to satisfy 
the requirements of the corpuscular theory of 
W. F. G. Swann for the soft component of the 
cosmic radiation are protons which produce 
secondaries in the numbers which are actually 
observed as showers. In conclusion, the authors 
wish to express their deep appreciation to Pro- 
fessor Swann for much valuable discussion, 
advice and criticism. 

J. R. Oppenheimer, Phys. Rev. 47, 44 (1935); L. W. 
Nordheim, Phys. Rev. 49, 189 (1936); W. F. G. Swann, 
Phys. Rev. 49, 829 (1936). 


2C, D. Anderson and S. H. Neddermeyer, Phys. Rev. 
50, 263 (1936). 
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F In I we discuss the status of the quantum theoretic 
formulae for pair production and radiation in the domain 
of cosmic-ray energies, and the relevance of these processes 
to an understanding of showers and bursts. In II we give 
a qualitative estimate of the course implied by the theory 
for a shower or burst built up by multiplication from a 
very energetic primary; we then set up the diffusion 
equations for the equilibrium of electrons and gamma-rays, 
and show how these can be simplified. In III we carry 
through the analytic solution of the diffusion equations, 


I 


N nuclear fields, gamma-rays produce pairs, 

and electrons lose energy by radiation. The 
formulae which have been deduced! from the 
quantum theory give for the probability of these 
processes values which, for sufficiently high 
energies, no longer depend upon the energy of the 
radiation. Because of this, the secondaries, pro- 
duced by a photon or electron of very high 
energy, will be nearly as penetrating as the 
primary, so that the primary energy will soon be 
divided over a large number of photons and 
electrons. It is this development and absorption 
of showers which we wish to investigate. 

The finite limiting cross sections for radiative 
loss and for pair production essentially limit the 
penetrating power of electrons and photons; as 
we shall see, 20 cm of Pb should absorb prac- 
tically all such radiation if the primary energies 
are <10° Mev. From this one can conclude, either 
that the theoretical estimates of the probability 
of these processes are inapplicable in the domain 
of cosmic-ray energies, or that the actual pene- 
tration of these rays has to be ascribed to the 
presence of a component other than electrons and 
photons. The second alternative is necessarily 
radical; for cloud chamber and counter experi- 
ments show that particles with the same charge 
as the negative electron belong to the penetrating 
component of the radiation; and if these are not 


1 An account of the results and of the theory upon which 
they are based may be found in Heitler’s book, The 
Quantum Theory of Radiation (Oxford, 1936). 

? 


and find the distribution of electrons and gamma-rays 
as a function of their energy, the primary energy, and the 
thickness and atomic number of the matter traversed. 
We treat the effect of ionization losses on the shower, 
calculate the amount of radiation of low energy to be 
expected, and treat transition effects in passing from one 
substance to another. In IV we discuss the results of the 
calculations, and give a summary of the conclusions to 
which they lead, and the difficulties. 


electrons, they are particles not previously known 
to physics.” 

Direct evidence for the approximate validity 
of the theoretical formulae is provided by the 
latest studies of Anderson and Neddermeyer' on 
the energy loss and pair production of electrons of 
energy up to 400 Mev. This evidence is still 
incomplete; yet it affords absolutely no indica- 
tion of a breakdown of the theoretical formulae. 
Since there is good evidence from the altitude and 
latitude curves of cosmic-ray ionization, as well 
as from the transition curves for showers and 
bursts, of a component in the cosmic rays which 
is strongly absorbed and yet has a very high 
energy, it seems of interest to investigate in detail 
the consequences which the theoretical formulae 
imply for the degradation, multiplication and 
absorption of such radiation. We shall find in this 
way a model for the building up and absorption 
of large showers and bursts which in many impor- 
tant respects agrees with what is found experi- 
mentally.* From this we should like to derive on 
the one hand a further argument for the qualita- 
tive validity of the theoretical formulae, and on 
the other for the often repeated suggestion that 
many showers are built up by a long succession 


2 Thus Williams first suggested that penetrating cosmic 
rays were protons, positive and negative: E. J. Williams, 
Phys. Rev. 45, 729 (1934). Cloud chamber evidence would 
favor a particle of smaller mass. 

’C. D. Anderson and S. Neddermeyer, Phys. Rev. 50, 
263 (1936). 

* An account of the experimental findings on showers and 
bursts and their transition effects may be found in Geiger’s 
article, ‘‘Die Sekundaer Eftekte der Kosmischen Ultra- 
strahlung,’”’ Naturwiss. (Leipzig, 1935), 
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of simple elementary processes, and not by the 
simultaneous ejection of a huge number of par- 
ticles in one elementary act. 

Here a certain caution is necessary. Cloud 
chamber observations have shown the existence 
of two fairly well differentiated types of shower.® 
In one of these, and by far the more common, 
only electrons, positrons and y-rays appear to 
take part; the shower particles are, except for 
those of very low energy, well collimated, with 
transverse momenta of the order of a few million 
volts ; the showers can often be seen to increase in 
passing through matter, and typically, if they are 
large, have no well-defined focus. It is these 
showers for which our calculations will give us 
some understanding. In the other and rarer type 
of shower, transverse momenta of the order of 
100 Mev are common; the shower is usually not 
collimated at all; heavy recoil particles are fre- 
quently seen ; and the total number of particles is 
usually small. It is natural to ascribe these 
showers to the interaction of heavy and light 
particles, and to accept the arguments which 
Heisenberg® has advanced to show that for such 
processes the probability of ejection of several 
particles should be of the same order of magni- 
tude as that of one. It would seem, however, that 
Heisenberg’s attempt to interpret on this basis 
the larger showers and bursts as highly multiple 
elementary processes is without cogent experi- 
mental foundation ; and we believe that in fact it 
rests on an abusive extension of the formalism of 
the theory of the electron neutrino field. 

One can then believe in the applicability of the 
following analysis to cosmic-ray photons and 
electrons and their showers, only if he admits the 
presence of another component to which the 
analysis is not at all applicable, and of other 
types of elementary processes, which essentially 
involve the heavy particles and their coupling 
with electrons, and which find no place in this 
treatment. 


II 


We are interested then in what happens when 
an electron (or positron or photon) of very high 
initial energy passes through matter. We shall 





5’ Anderson and Neddermeyer, see reference 3; Brode, 
MacPherson and Starr, Phys. Rev. 50, 581 (1936). 
®W. Heisenberg, Zeits. f. Physik 101, 533 (1936). 
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consider only the three elementary processes of 
pair production by photons, radiation by elec- 
trons, and ionization losses by electrons, the first 
two because they dominate the multiplication 
which makes the shower, the third because it 
limits the size of the shower and absorbs it. We 
shall not consider relatively rare multiplicative 
processes, such as the production by electrons of 
high energy electronic secondaries, the direct 
production of pairs by electrons, or the Compton 
effect. We shall suppose that all high energy 
particles come off forward, and neglect the 
angular divergence of the shower. We shall not 
try to treat in detail radiation of such low energy 
that these simplifications are invalid. 

To define the probabilities of the elementary 
processes, we shall use simplified formulae which 
approximate closely to the limiting forms which 
the theory gives for high energies.! 

Thus for the probability that an electron (or 
positron) radiate an energy in the range E, E+AE 
in passing through a thickness Ax of matter of 
atomic number Z, nuclear density N, we shall 
take 

KAxAE 4Z*e®N 200 
PAxAE = with K= In . (1) 
E hm?C® Z} 


For the probability that a gamma-ray of energy 
E make a pair of energy E’, E—E’, we take 


E'<E; 


K’'/K=a~ j. (2) 


P’AxAE’' =(K’AxAE’')/E; 


It is convenient to measure length in terms of a 
variable t=xK. For Pb the unit of length is ~} 
cm, for water it is about 0.4 m. For elements as 
heavy as Pb, the constants K and K’ can hardly 
be regarded as known within 20 percent. Formula 
(1) is a good approximation as long as the elec- 
tronic energy is >20 Mev; formula (2) begins to 
give too large results for E<50 Mev, but is off by 
only a factor 3? at E=25 Mev. As for ionization 
losses, we shall suppose them to be independent 
of energy, and evaluate them for those energies, 
where their effect will be important. Thus we 
shall write for electrons 


dE 4rNZe* 8 


=-~, B=—— In —. (3) 
ot KmC? ZRh 








AND 


bdo 
bo 
No 
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For Pb, 8 has the value 6.5 Mev. It is quite 
closely inversely proportional to Z. 

For primary energies which are not too high, 
one can carry through the calculations step by 
step, finding how many gamma-rays are pro- 
duced by the primary, how many of these are 
absorbed by pair production, how many gamma- 
rays these in turn produce. This procedure was 
used in computing the number of pairs theo- 
retically to be expected in the electron traversals 
studied by Anderson and Neddermeyer; it has 
been used by Nordheim, and by Heitler and 
Bhabha,’ though here with neglect of ionization 
losses. It clearly becomes prohibitively laborious 
for thicknesses and energies of the order of those 
involved in large showers and in bursts; and we 
shall try instead to solve the diffusion equations 
implied by (1), (2) and (3), and so obtain an 
insight into the course of the shower for large 
thicknesses and high primary energies. An ex- 
tremely rough indication of what we may expect 
to find we can see from a quite simple argument. 

Both radiation and pair production give two 
rays for one in the shower. Each process of this 
kind has about an even chance of happening in a 
distance t=1 (somewhat less than even for pair 
production, somewhat more for y-radiation). 
Thus the order of magnitude of the total number 
of electrons and y-rays to be expected at a 
thickness ¢ is 2‘. The order of magnitude of the 
energy loss in thickness At is thus ~ }82‘At. When 
the integral of this is equal to the primary 
energy Eo, the shower will be absorbed. Thus if 
T is the distance to which the shower penetrates 


B2T~2E, In 2. (4) 


From this it follows that T will increase log- 
arithmically with Eo and that it will decrease 
slowly with decreasing Z; that the number of 
particles in the shower will increase with Eo 
roughly linearly, and will be roughly proportional 
to Z. For showers of about 30 particles we should 


7 We are indebted to Dr. Nordheim for writing to us of 
his results. See also Heitler and Bhabha, Nature 138, 401 
(1936). We are further indebted to Heitler and Bhabha 
for sending us a manuscript of the paper in which they 
have extended these calculations. Their results differ from 
ours primarily because of their neglect of ionization losses; 
apart from this the agreement between their values and 
ours is excellent. We do not agree with their conclusion 
that these calculations make it possible to ascribe the 
greater part of sea-level cosmic radiation to degraded 
electrons and photons of high initial energy. 
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expect the maximum in the shower to come at 
T ~ —loge 30~5, or a little over two cm of Pb or 
about 7 cm of Fe. These are in fact‘ of the order 
of the distances at which maxima are found for 
the transition curves for relatively small showers. 
For bursts of 1000 particles the transition maxima 
are found® at about twice as great thicknesses, as 
our rough estimate would suggest. 

To set up our diffusion equations, let us write 
y(t, E)AE for the probable number of gamma- 
rays to be found at a thickness ¢ in the energy 
range E,E+AE and @(t, E)AE for the cor- 
responding number of electrons and positrons. 


Then 


OY 1 ;*o : 
—=-—oy+ | P(t, ede, (5) 
ot Fe E 
ely Fo y(t, €) AP 
=20 f —de+-- 
ot E . de 
Fo P(t, €) rE de 
+f ——§de-—O(t, E) —. (6) 
[E] e—E to) € 


The last two terms in (6) can be combined to give 
the clearly finite result 


—s Fo P(t, €) oF de} 
R=lim 7 -de-O(t, E) J - 


8.0 ‘4x43 €E— 5 «| 


E 
= — @(t, E) ln ——— 
Ey-—E 


Fo AP € 
-f —/(t, e) In (<-1)a (7) 
Jr O€ E 


and give the change in pair distribution which 
comes directly from the radiative losses of the 
charged particles. The term in 6 gives the cor- 
responding change from ionization losses, and the 
first term on the right in (6) gives the pair pro- 
duction by y-rays. In (5) the first term gives the 
absorption of y-rays by pair production, and the 
second term their replenishment by radiation. 
These equations are to be solved for the boundary 


conditions 
y(t=0, E)=0; O(t=0, E)=A(E, Ep), (8) 


where A vanishes except when E is in the 
immediate neighborhood of Ep», and has the in- 
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tegral 1. It will be convenient, however, to take 
A(Eo, Eo) =0, and to choose for A, not a delta 
function, but 


a” Eye 1 
A(E, Eo) =lim —(In E/E)" (—) —, 
ano F Eo/ T(n) 


n=4, (9) 


The reasons for this choice will be apparent as we 
develop the solution. 
From (5) we may write 


O(t, E) =e—*'(02z) /(0t0E) with y= —e~*'s/E. (10) 


It is not possible by differentiation to reduce 
(6) to a differential equation, because of the 
occurrence of E in the integrand of (7). We have 
therefore tried to replace the terms (7) by others, 
which would lead to a readily soluble system of 
differential equations, and which would still give 
a good representation of (7). As we shall see, we 
can obtain a solution in terms of integrals of 
elementary functions provided we write in 
place of (7) 


rE) P(t, €) 
Rw | ——“de+- b(t, E) 
“rE € 


9 


ae a0 . 
(t, E)+dE*—+:--- (11) 
dE? 





+ckE ‘ 
if a, b, c, d are constants. Two elementary con- 
ditions now limit their choice. 
1. The number of particles leaving a given 
energy range by radiation must be equal to the 
number entering some other range. 


RdE=0. 


0 


2. The energy lost in Ax must be equal to that 
appearing in the y-rays 


Eo Eo 


REAE= { P(t, E\dE. 
0 0 


If these conditions are to hold for any @, then 
c=d=0, etc. From (1) a+)=0; and from (2) 





a+b/2=1; thus we must have a= —)=2 and 
Eo P(t, e€) 
R=2 de—20(t, E). (11a) 
gE € 


The use of this form for R corresponds to assum- 
ing a uniform distribution of energy losses instead 
of the 1/E law given by (1), and renormalizing 
to give the correct total energy loss. The produc- 
tion of y-rays is, however, treated correctly in 
accordance with (1); it is only the effect of these 
losses on the redistribution of the electrons them- 
selves that is falsified. Since it is difficult to give 
an a priori estimate of the error thus introduced 
into the solution, we have used the solutions 
which we have obtained to compare the values of 
R given by (7) and by (11a). When ¢ is not too 
small, most of the particles have of course an 
energy E<Ep. For this case we found that the 
two values of R agreed within less than 5 percent. 
It is of some interest to ask why this should be. 

The answer is to be found in the circumstance 
that for Epo, our energy distribution curves 
for @ follow quite closely a 1/E* law as Fig. 4 
shows, and as cloud chamber observations sug- 
gest that they should. For this law 


O(t, E)~k(t)/E? 
and for EXEo, both (7) and (11a) agree in giving 
= —k(t)/E?. 


The approximation of replacing (7) by (11a) is 
thus a very good approximation in the range 
where neither ¢ nor In Eo/E are too small ; and we 
shall see that it is only here that our results are of 
direct physical interest, because for small ¢ and 
E~ Eo, the fluctuations to be expected from the 
probable behaviors defined by @® and y are 
all-important. 

Using (11a) for R, and differentiating (6) we 
get, if we use (10), the differential equation for z 


> 


24+ (2—o0)3"+22'/E—202/E?—B2'"=0 (12) 


or, with o= 3, 


2 2e2 
3 +202""+ - 2’ -——- —82'" =0. (12’) 
E E? 


Here primes denote differentiation with respect 
to E and dots with respect to ¢. Eq. (12’) is to be 
solved with the boundary conditions 


2'(0, E)=A(E, Ey). (13) 


2(0, E)=0; 
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Ill 
Let us now write 
A=In (E/E) 
and try to solve (12’) by a Laplace transfor- 
mation: 


z=1/(272) | (dy/y)e"F,(t), 


¢ 


(14) 


where C is some suitably chosen contour in the y 
plane. From (12’) we then find 


1 dy °F, 
| eorwperrh 


Qri cy or? 
OF, ' | 
+[2oyWy+1)—2y }+——2eF, 
ot 
B ¥ OF, 
=—— J evo+0(9-+2)0% . », (15) 
2riky “C at 


where 8/Ep is a very small number. If we neglect 
these terms, then 


Fy, = M,e*'+N,er™*, (16) 
where yu, v are the roots of 
K°+2(06—1/y+1)x—20/y(y+1) =0, - 
(17) 
"- —ot+1/y+1+[(o—1/y+1)?+20/y(y+1) }! 
Vv 


and for R(y)>0, u>v. The constants M,, N,, 
and the contour C must be chosen to satisfy (13). 

It is clear that the integrand of (14) will have 
branch points in the negative half y plane, and 
that it will -0 as |y|—>0. We therefore choose 
for our contour C a closed path, consisting of a 
straight line 6-71“ to 6+7z, and the infinite 
semicircle from 6-+72 through the negative half 
plane to 6—1. From (13) we get 


| ay ye"F,(0) =0, 


e Cc 
| dye" F,/at= —2xiEA(E, Ey). 
co 


If we multiply (18) by e~™, with R(yo) >6, and 
integrate over \ from 0 to ~, we get 
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> dy F,(0) 1 ¢ dy @OF,(0) 
pansies | : 


’ = — $(yo) 
Ycoy—-VYo Y 2riJc y—Vo at 
_ 
with g=lim | e- Wd —_ “Ipaq y 
we « '(n) 
(19) 


=lim (1+ ¥0/a) 


a--co 


Since F, and 0F,/dt-0 as | y|—, and have no 
singularities in the right-hand half-plane, this 
gives 

(20) 


F,(0)=0; o&F,(0)/dt= —¢(y). 


These are the boundary conditions for F. Any 
other function A which gives an F,-0 and 
OF,/dt-0 as them without 
singularities in the half-plane R(y) >0, and whose 
limit gives A(E, Ey) =0 for E+ Eo, could clearly 
be used in place of (9). But for the treatment of 


y|—o, leaves 


ionization losses it is convenient to have the first 
3 derivatives of A vanish as EE. 

Let us now carry through the solution for B=0. 
For this solution we write 20, yo, ®o. From (20 


M,=-—o/(u-v), N,y=¢/(u-v), (21) 
1 . dy e+! — et 
3=- | ery >. (22 
2midc v u—v 
From this we can readily calculate the total 
energy of gamma-rays and electrons: 
»Eo 
| EdE[ @0(t, E)+yo(t, E)] 
“oOo 
(23 
1 edy @ (uvt+l)e*—(vy+1)e” 
= Eye“! | | 
2miv¥c y 1—-y uv 


The integral on the right gives the residue of the 
integrand at y=1, which is just, as a—>~, 


2rie®' since p=o, v=—1 when y=1. 


Thus the total energy remains constant, as it 
should, when ionization losses are neglected. For 
this the condition on the constants in (11) 
a+6/2=1, is essential. 

When ¢ or ) are large, then (22) may most con- 
veniently be evaluated by the saddle point 
method. The approximation here involved is one 
which corresponds closely to the limitations im- 


0) 


e 


2 


ON MULTIPLICATIVE SHOWERS 


posed upon the physical interpretation of our 
solution by the fluctuations which must be ex- 
pected in the actual behavior of the radiation 
from the average behavior given by our diffusion 
equations. The effect of these fluctuations can be 
simply formulated by the physically obvious 
assertion that the addition or subtraction of a 
thickness of matter corresponding to t~1 has an 
even chance of not altering the actual distribu- 
tion of pairs and gamma-rays. The asymptotic 
form of @, y for large ¢ thus gives us results whose 
accuracy corresponds to the applicability of the 
diffusion equations themselves. 

The unique saddle point of the first term of the 
integrand of (22) lies on the positive real axis, at a 
point which moves out to larger y values mono- 
tonically as ¢/X is increased. If we use a subscript 
s to indicate that functions are to be evaluated at 
this saddle point, then we find from this term 


exp (Ay. +s) 


29 = — (2ry,"’t) 
Vs\Ms— ts) 
4 

‘4 . (i+) 
Yo= — Zoe atta Eo, 
Po = ul sVs¥ 0- 

Here "= y/dy’. 
MM é 


From the second term of (22) 


1 dy e 
| 


2miv/¢ 7, PP 


Autre 


we obtain a contribution which remains less than 
Qe-'t, where Q is independent of \ and t, ast. 
For large ¢ this term is thus quite negligible 
compared to (24). 
For \>?, and \<t?, these expressions (24) can 
be evaluated analytically. Thus (writing ¢ = 2/3). 
For >t 


~ 


o= 2x0 exp [M+ 3Inie, 


(25) 
MsVs=2-3-§D-M2h, 
and for \<t, 
So = 2-83 7-134t-1X-1 exp [3-2-IN'4], 
Me¥e= 2-4. (26) 


For intermediate values we give a plot of y, 
against A/t, from which (24) may be evaluated. 
It should be observed that yu,y, gives directly the 


i] 
to 
wn 


ratio of ®o to yo, and that it never exceeds 2/3. At 
any energy and thickness ¢>1, there are always 
more y-rays than electrons. 

The solution we have given can be readily ex- 
tended to diffusion equations with R of the 
general form (11); only the dependence of » and 
v on y becomes more complicated. 

For the treatment of the ionization losses given 
by the terms in 8, it is simplest to return to (14), 
and write instead 


1 dy 3 
= of FW Fy (2) 
ride y E 
B° 
+— Ff (t) (27) 
| Dia 


and leave the contour C unaltered. If for F,(2) 
we take (16), (21), then (12’) gives, for the first 
order terms in 8. 


07 F,@ 1 OF, 2o 
+2 (< — ) = FY 
ot? yt+2 at (y+1)(y+2) 


1 y * dy'(y'+1)(y’+2) 


EE ME | yy 


; ’\*% ; ’ , t 
RT err eH aS 
. 2g) 
A ; ; Q. (Le 
uy )—vly) 


The general solution of this is 


Fy (t)/y=A yee'§ + Byer '§ + Cyew (ute 
+Dyerr+ (29) 


u(ot+uy)(y+1)(v+2) 
where A,=¢ ; : 
2(2o0+yy)(u—v) 

(29a) 
v(o+vy)(y+1)(y+2) 
B,=—-¢- 7 ’ 
2(20+ vy) (u—v) 


and where C, and D, must be chosen to maintain 
the boundary conditions (13), which require 


A yuly) + Byr(y) — »(y+1)(A,+B,) 
v(y+1)—w(y+1) 


(29b) 
A yu(y) +B,r(y) —plyt+ l \(A,+B,) 


D,y= ; 
u(y+1)—v(y+1) 
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For R(y) > —1, C, and D, are analytic, except for 
a simple pole of C, at 


VYp=(3—1/ 73) /8~ —0.69 


and a simple pole of D, at 


Vp =(3+17/73)/8~ 1.44. 


7 | 


The terms of order 8 in (27) may then be written 


B ’ « 
| dye ity TA Ors? ‘+ Bie’ 
qe , 
+Cye" y+ DI+ Die’ y+1 '}. 
For large \ and ¢ only the term 


. 8 ef 
J dye+%A yen = 


Eo2rie Cc 


is important. 
For the second and fourth terms in (30) we can again 


show that they remain less than 


B 
E Qed-2/38t as tw, A> 
* 


and are therefore negligible. 
The term 
B — 
r= | dyCye" 
Eo2mivC 
may be evaluated by deforming the contour C to pass 
through the saddle point at 
Vo = —1+ 4s. 

If we call the path taken along a line parallel to the 
imaginary axis through j~, = then 


_ 
pa _ pf dycyerrvonartne 
2riEo v= 
+—3(Ve— yy Gre l+yp)A+H(1+yp)t 
“0 
where d(x)=0 x>0 
4 x=0 
1 x<0 
and Gyp=lim (eC pte) 


e—0 


Each of the terms in T gives a contribution negligible for 
\—> «© compared to (31). From the saddle point integral 
we get a contribution 
8 
<Q(A, t)—eusd+Het, 
Eo 
where Q is an algebraic function of A and ¢; from the 


residue we get 
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Fic. 1. Plot of \/t against y,. 


with Q, a constant. Both of these are small compared to 
the term 


which we get from the saddle point integral (31). 


From (31), and the application of the saddle 
point method to the integral occurring in it, we 
can readily see by what factors the values of 2», 
Yo, Po given in (24) must be modified: 
2=2o(1—8r,./E); y=yo(1—8r,./E); 

P=P(1—B(1+1 
P=yusy;(1—8r,/y,F), 
y(y+1)(v+2)(o+pny) 
2(20+uy) 

From (24) and (32) we can thus give @, y as 
functions of ¢ and E. The ionization correction 
still further increases, of course, the proportion 
of y-rays in the radiation. 

From (32) we see that our approximate calcu- 
lation of the effect of ionization losses is based 
upon the smallness of 8/E. This parameter is 
about 12 percent for Pb, 35 percent for Fe, at 
E=50 Mev. Clearly the treatment here given is 
limited to energies E>, but within this limit it 
can tell us what the effect of these losses will be. 

In fact (22) and (31) give the first two terms in 
the asymptotic expansion in 8/E of a solution 
of (12’): 


¥e)7t-/E); 
(32) 


where T= 


° a 


1 dy exttry . 
|- i -y) f as —s)%e~* 


4’ oy u-v s 
Buygse 
xaFi( y+2; ytit+g; oc ). (33) 
2cE 
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Fic. 2. Plot of total number of electrons N(t) against ¢, 
for air, Eg = 2500 Mev, computed from (36). The circles are 
from the experimental results of Pfotzer,* on the variation 
of vertical coincidence counting rate in the upper atmos- 
phere, at magnetic latitude ~50°N. At this latitude the 
earth’s field will just admit electrons of energy 2500 Mev. 
The deviations from the curve for {>8 may indicate the 
presence of some electrons of higher energy, as well as a 
penetrating component. 


Here the contour S in the s plane is a simple loop 
from + counterclockwise around the origin 
and back to +; »F; is the hypergeometric 
series and 
g=o(1+y)/(o+uny). 

That (33) is a solution of (12’) may be verified by 
direct substitution. It is not of course the solution 
satisfying the boundary conditions (13); but an 
iteration of the arguments just given in connec- 
tion with (31) indicates that the terms which 
must be added to satisfy (13) give for ¢>1, 
EE, a negligible contribution. 

For heavy elements such as Pb, (31) gives an 
adequate approximation to (33) for all energies 
high enough to make the use of (2) permissible: 
the pair production formulae break down before 
the ionization correction becomes very large. For 
light elements, for the atmosphere, we may apply 
(33) for E<. To do this we use the analytic 
continuation of the hypergeometric series, and 
evaluate the integral over y by the saddle point 
method. For E/8<1 we thus find 

(2ap,"'t)- , 
VE, #) = ehevet(oe-o)t___. 


Vs\Ms— Ve) 
r(i+y.+g:;) 1 
x- 1s — 7 : at 
r(i+y,)r(2+y,)l(1+g,) £ 


(34) 
2cE 





O(E, t)= ¥(E, t) In 
E 


* Zeits. f. Physik 102, 41 (1936), Fig. 1, p. 42. 


bo 
bho 
~ 


Here Ag=In (20Eo)/(Busyegs); and the saddle 
point y, is given as the real positive root of 
d 20k | 


y In- +t) =0 (; 
dy| Buyg 


=) 
uw 


and may -be found from Fig. 1. 

The energy dependence of y and ®&, which for 
E>8 is roughly given by k/E?, becomes, if E <8, 
k,/E, ke|n B/E, respectively. The increase in 7 
and & with decreasing energy will be still further 
reduced at still lower energies by the absorption 
of the y-rays by Compton effect. The total 
number of electrons in the range 0—8 depends 
upon 6 essentially as 

ess 


and its maximum as a function of ¢ is quite 
closely inversely proportional to 8 or proportional 
to Z. 

The total number of electrons is easy to com- 
pute from (33): 


»Eo 
N(t)= | O(E, dE = —e~*'z(E=0, t) 
“o 


pe(2rp,"’t) 


r(i+y.+¢:s) 
7 ys(us—vs) T(1+y,)0(2+y,) (14+ 2.) 


x eravet (us ot (36) 


where y, is given again by (35). The maximum 
value of N(t) occurs for values of ¢ slightly smaller 
than Ag. We give Fig. 2 a plot of N(t) against ¢, for 
water, for which we have taken B=90 Mev, and ¢ 
is measured in units of ~0.4m. The primary 
energy is ~2.5X10* Mev. The maximum of NV 
occurs at about 1.2 meters. It would seem prob- 
able that the latitude sensitive transition effects 
of the cosmic-ray ionization in the upper atmos- 
phere are to be interpreted on this basis. 

It is futile to apply (33) to elements as heavy 
as Pb, because for energies of the order of 8 the 
radiative formulae (1), (2) and the diffusion 
equation based on them become quite wrong. As 
is well known, the actual behavior of radiation in 
the range 1-25 Mev in Pb is extremely compli- 
cated. Since the absorption coefficient for y-rays 
has a minimum value which for Pb lies at about 
3 Mev, we know that for sufficiently great 
thicknesses of matter y-radiation of roughly this 
energy will predominate over y-rays of higher or 
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Fic. 3. A plot of the energy degraded into rauaiation of 
energy <25 Mev per unit ¢, against ¢. The plot is for 
Pb, Eo =2X 104 Mev. The abscissae are in units of } cm, the 
ordinates in units 100 Mev. The number of low energy 
electrons to be expected at ¢ is roughly 0.15 Do;(t—3). 


lower energy and over electrons of all energies. 
A rough estimate of the number of low energy 
electrons we can get by computing the total 
energy degraded into the low energy region per 
unit thickness. Let us call the energy so degraded 
into an energy region 0<E <e, D,(t). This energy 
has ultimately to be absorbed by the ionization 
losses of electrons. The electrons 
which in a distance t= 1 would lose by ionization 
an energy D,(t) is 1/8D,(t); it must be remem- 
bered, however, that because of the relatively 
great penetration of the low energy gamma-rays 
(amounting for Pb to ¢=3-4), the low energy 
electrons actually present in the radiation at 
some fixed thickness t= T must be evaluated, not 
from D,(T)/8, but from some appropriately re- 
tarded value D,(7 —r)/8, where r gives the effec- 
tive mean penetration of the low energy radia- 


number of 


tion, and may be of the order of 2—3 for Pb. It is 
in this sense that the curves we give for D,(t) are 
to be interpreted. For the interpretation of ob- 
servations made with the high pressure, thick- 
walled ionization chambers which are used in the 
study of bursts, the variation of D,(t) with ¢ 
should itself give a valid estimate of the probable 
variation of the recorded size of the burst with 
absorber thickness. 

The calculation of D,(t) is straightforward: 
The total energy lost by the radiation of energy 
>e per unit ¢ is 

d p* 
D(t)= —- [O(E, t)+y(E, t)JEdE. (37) 


dt/, 
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The part of this which is accounted for by 
ionization losses is 
aEo 
D.(t)= | BOE, NaE. (38) 
Then Dt) =Dt)hY? —D (t)™. (39) 
In this way we find, from (39), (32) and (24) 
(us—o)(1 +595) 
D,(t) = —e~*'Z € 
| l-—y, 
—[(o—ws)[1+(1+y.)us lyst +s s' . (40) 





Here 2 is given by (24), and is to be evaluated, as 
are ws, y, and 7, for E=e, \=In Eo/e, t=1. It will 
be observed that the effect of the terms in 6 in 
(40) is to increase D,(t) when \>4, and to de- 
crease it when \<X<?¢. We give in Fig. 3 a plot of 
Dt) for Pb, e=25 Mev, Eo=2X10* Mev. 

Let us now restrict ourselves to E> 8, and ask 
how the shower will be altered by changes in the 
initial radiation. 

A. If the incident radiation is a y-ray of energy 
E, and not an electron, the boundary conditions 
(20) must be altered to 


F,(0)=$¢y; OF,(0)/at=0. (41) 
This gives us 
1 > ve! — ye”? 
207] =— | dye", (42) 
2rie/ ¢ 


uMavY 


from which it follows that for ¢>1, the number of 
gamma-rays and electrons will be changed by a 
factor — v,y, from the values given by (32). This 
factor varies from 0.6 to 1.6 over the range of 
values of ¢/ for which the shower is appreciable. 
The course and magnitude of the shower will 
thus depend very little on whether it is started 
by an electron or a gamma-ray. The maximum of 
the shower will come for very slightly higher ¢ 
values when it is initiated by a y-ray. 

B. If the incident radiation is a group of m 
electrons of energy Eo/m, then from (24) 


Pofm) =m? Po(A—In m, t) ; 


Yo(m) =m*yo(A—In m, t) 
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and insofar as @(X, ¢) is quite closely of the form 
k(t)e, for small variations of X, 


©, :=0 
O'ofm] = 0'o 


for small m, and the shower will be substantially 
the same as if one electron of energy Ey had 
started it. 

C. As long as we neglect the ionization terms, 
we can get no “transition effects’’ when a shower, 
built up in one substance, passes into another. 
But the terms in 8 do give us such effects, in 
apparent qualitative agreement with what is 
found experimentally. Let the radiation go a 
distance ¢; (measured in units approximate to the 
first substance) in (1), and then a distance fg» (in 
units for (2)) in (2). Then the terms independent 
of B in zg are, for t; +¢.>1. 


1 *dy erhutu(irtte) 
t= — —| $ (44) 
2rivc y uv 
The boundary conditions 
O(E,t:) in (1)=@(E,t4) in (2); 
Y(E,t:) in (1)=y(E,4:) in (2) 


give then, for ¢;>1, for the terms in z propor- 
tional to B: 


1e 
ef dyge* (vt Ba yew (uy) trtea) 


2rie Cc 
u(y) —r(y+1) 
eX (ut 1)t2 


Had seorf 
u(y+1)—v(y+1) 


u(y) —u(y+1) 


ecvvonll, (45) 
u(y+1)—v(y+1) 


which for ¢2=0 is just the ionization correction 
appropriate to (1), and which for tz ap- 
proaches that appropriate to (2). If Z:>Z2; 
8,<B2 and @ is decreased by the transition; if 
Z,>2Z;, @ is increased. The transient terms which 
give the transition fall off, relative to the main 
terms, as 


elu(ust1)—ulys)]t2 ele (ustl)—v (ys) } te, 


and 


Thus the second term will be damped out for t2~3, 
but the first term will persist until tg~1.5—2, the 
transition thickness is therefore of this order of 
magnitude. For Ey=2X10* Mev, 4:=6 the 
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Fic. 4. A logarithmic plot of energy distributions for 
shower electrons, made for Pb, Eo=1.5X10° Mev and 
t=15, ¢=8, t=6, and t=4. Abscissae are A=In E/E; 
ordinates are natural logarithms of @ measured in (Mev)~. 


A straight line of slope 2 would correspond to a 1/E 
distribution law. 


transition from Fe to Pb should increase the 
number of electrons with energy >50 Mev by 
about 35 percent. 

It must be remembered that the actual transi- 
tion effects observed are complicated by the fact 
that the high energy radiation which produces 
the showers is to some extent being regenerated 
in the material in which the showers are formed. 


IV 


In Fig. 4 we give typical energy distribution 
curves for the pairs, for Pb, and a primary energy 
of 1.5105 Mev for t=15, 8, 6, 4. In the neigh- 
borhood of \=8, the ionization corrections are 
appreciable, but for \<6 the curves will hold for 
all substances. It will be observed that these 
curves follow quite closely what we would get 
from a distribution law k/E?; thus except for !=4 
they may all be represented by a law of this form 
in which the exponent of E never differs from —2 
by more than five percent. In fact this law gives a 
very good approximation to the energy distribu- 
tion for ¢>.7\, and for all energies EXE, but 
>. We have already pointed out that it is for 
this reason that the simplification of the diffusion 
equations made in replacing (7) by (11) is per- 
missible. Cloud chamber studies® of the energy 
distribution of shower electrons seem to fit this 
law quite well for E>25 Mev, but the experi- 
~ 8C. D. Anderson and S. Neddermeyer, Int. Conf. 
Physics, London 171 (1934). 
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mental conditions here are clearly not such as to 
make this agreement very significant. 

Because of the form of the distribution law, a 
plot of @E will give a good estimate of the total 
number of electrons to be expected with energies 
=E. For Pb, E=50 Mev, and several values of 
E> varying from 2700 Mev up to 1.1X10® Mev 
we give such plots against ¢ (Fig. 5). The rapid 
increase with Ey of the total number of particles, 
and the slow shift in the position of the maximum, 
confirm the qualitative arguments of §2. We 
would like to interpret in terms of this shift the 
well-known observation that the optimum transi- 
tion thickness for bursts is considerably larger 
than for showers. In fact, the agreement between 
the experimental values of the optimal transition 
thicknesses (~2 cm Pb for showers, ~5 cm Pb 
for bursts) with the position of the maxima of 
these curves, seems to us a strong argument for 
the correctness of the model we are treating, and 
of the validity of the high energy formulae we 
have used. In particular the experiments of Nie,°® 
in which it is shown that a burst generated in a 
suitable layer of matter may sometimes have its 
magnitude very much increased by the inter- 
position of a few cm of Pb find a very natural 
interpretation in terms of these theoretical 
curves. It must be remembered, however, that 
the experimental transition curves, which give 
the dependence on ¢ of the probability of finding 
a shower or burst whose magnitude exceeds a 
lower limit defined by the experimental arrange- 
ment are not strictly comparable with the curves 
of Fig. 5 which give the variation, for fixed Eo of 
the probable number of electrons." For in making 
this comparison the initial distribution of the 
radiation over Ey must clearly be taken into 
account. We want here, too, to point out that the 
actual behavior of the radiation will fluctuate 
about that given by our curves, and that the 
order of magnitude of the probable fluctuations 
can be estimated by shifting the curves by 
At~ +1. 


® Nie, Zeits. f. Physik 99, 776 (1936). 

10 When counters are used to detect showers, the angular 
divergence of the shower rays is necessarily exploited. If 
the theory here developed is at all correct, most of these 
rays must have a relatively low energy, of the order of 
several million volts. In fact only cloud chamber observa- 
tion can tell us much about radiation of much higher 
energy. 
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Fic. 5. Plots against t of @-E for E=50 Mev in Pb. 
(a) Eo=2.7X108, (b) Eo=2X10', (c) Eo=1.5X 105; 
(d) Eo=1.1 10°. Abscissae are in units of 4 cm; ordinates 
in units of 1, 4, 10, and 50 for (a), (b), (c) and (d), re- 
spectively. These plots also give the number of electrons of 
energy 250 Mev to be expected as a function of ¢. 


The application of our methods to Eo as low as 
2700 Mev may seem unjustified. Here our results, 
however, agree reasonably with those computed 
by Bhabha and Heitler™ with neglect of ioniza- 
tion losses and of processes of high order. For 50 
Mev the inclusion of ionization losses reduces 
the number of electrons to be expected by 15-20 
percent. 

It may be helpful to give a brief summary of 
the general results. For any absorber we measure 
length (¢) in units which are proportional, roughly, 
to Z*p/A, where Z is the nuclear charge, A the 
atomic weight, and p the density, and we define 
the characteristic energy 8 (see (3)), which varies 
about like 1/Z. Then, 

(1) The number of electrons per unit energy 
is about inversely proportional to the square of 
the energy, as long as B< EXEo, t>4n. 

(2) For given energy and ¢>1, there are 
always more y-rays than electrons; where the 


1 Reference 8. Heitler and Bhabha give for Ey =2.7 X 10° 
Mev, a value PE=4, maximum for t=4; we get PE=4.5 
at t=4. This discrepancy, and the fact that we find the 
maximum at ¢-~~5, are both in the direction to be expected 
from the effects of higher order processes. 
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ON MULTIPLICATIVE SHOWERS 


shower is near its maximum, this ratio varies 
from 1.5—2. 

(3) For E>8, the distribution curves plotted 
against ¢ are the same for all absorbers. 

(4) The number of particles of energy greater 
than E,>8 passes through a maximum for a 
value of ¢ which increases logarithmically with 
E/E, and is always quite close to Ing Eo/F. 

(5) The maximum number of particles with an 
energy less than some small multiple of 8 is 
attained for values of ¢ which are slightly smaller 
than In Eo/8, and decrease slowly with decreas- 
ing Z. The total number of particles of energy in 
this range is about inversely proportional to 8, 
or proportional to Z. 

(6) The maximum size of the shower is 
limited only by Eo with which it increases not 
quite linearly: thus an increase in Ey by a factor 
of 100 gives an increase in shower size of about 70. 

(7) If the initial energy Eo is in an incident 
y-ray, or is divided among a few electrons and 
gamma-rays, the course of the shower will be 
essentially unaltered. 

(8) Passage of a shower from one material (1) 
to another (2) will increase the size of the shower 
if Z2>Z;, decreases it if Z;>Z>». The transition 
takes place in a thickness t2~1}. All of these 
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results apply only for energies Ey above 10° Mev. 

In this paper we have altogether neglected the 
question of how such high energy electrons and 
y-rays can get down through the atmosphere. 
How serious this difficulty is we can see from (32), 
which tells us that for every electron of energy 
~2X105 Mev which hits the earth vertically, 
only 0.15 electron of energy >550 Mev will 
survive at the earth’s surface. This difficulty is 
made even sharper when we consider the form of 
the shower curves for great thicknesses ¢ > 30, or 
the showers and bursts reported under very great 
thicknesses of absorber. In fact, although when 
we go up far in the atmosphere, the showers, and 
still more markedly, the bursts, increase more 
rapidly than the total cosmic-ray ionization, 
below the atmosphere they do not fall off much 
more rapidly than this ionization. This suggests 
that, in addition to primary electrons and per- 
haps y-rays, which are able to produce multipli- 
cative showers directly, there is another cosmic- 
ray component, slowly absorbed, which is re- 
sponsible for the continuation of the showers 
under thicknesses of absorber to which no elec- 
tron or photon can itself penetrate. Some sug- 
gestions which we think relevant to the solution 
of this problem will be discussed in another 


paper. 
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The y-ray ionization current measurements in air at high pressures made by Erikson and by 
Clay and Van Tijn are shown to be largely modified by volume recombination at collecting 
fields below about 400 volts per cm. The Jaffé theory of columnar ionization does not explain the 
observed variation of ion current when it is applied over a wide range of collecting fields to data 


free from volume recombination. 


N 1932 one of the authors made ionization 
current measurements which were published 

under the above title.! Shortly after this Zanstra 
and Clay’ called attention to earlier measure- 
ments made by Erikson® in 1908. Later Zanstra‘ 
applied the Jaffé theory of columnar ionization 
to explain the curves obtained by Erikson, and 
Clay® and Van Tijn made new measurements to 
which the formula developed by Zanstra was 
applied. 

In the discussion of ion current measurements 
three types of recombination which prevent the 
collection of all ions formed are usually con- 
sidered, namely: (a) preferential, i.e., recombina- 
tion of the electron with the particular ion from 
which it was ejected, (b) columnar (called initial 
by some authors), i.e, recombination with a 
positive ion in the same §-ray track as the ion 
from which the electron was ejected, (c) volume, 
i.e., recombination with any positive ion. Ob- 
viously in processes (a) and (b) the probability 
of recombination of a given electron is inde- 
pendent of the number of tracks formed in the 
ionization chamber. Consequently if the pressure 
is kept constant, so that the number of ions per 
track does not vary, the fraction of ions collected 
is independent of the total ionization. On the 
other hand in process (c), i.e., volume recombina- 
tion, the probability of recombination of a given 
electron is proportional to the number of positive 
ions in the gas. This type of recombination is 
therefore negligible at low rates of ionization, 


11, S. Bowen, Phys. Rev. 41, 24 (1932). 

2H. Zanstra and J. Clay, Phys. Rev. 41, 679 (1932). 

3H. A. Erikson, Phys. Rev. 27, 473 (1908). 

‘H. Zanstra, Physica 2, 817 (1935). See also, B. Gross, 
Zeits. f. Physik 78, 271 (1932) and W. R. Harper, Proc. 
Camb. Phil. Soc. 28, 219 (1932) and 29, 149 (1933). 

5 J. Clay and M. A. Van Tijn, Physica 2, 825 (1935). 
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such as those found in measurements of cosmic- 
ray ionization, but becomes prominent at higher 
rates where it shows its presence by the variation 
in the fraction of the ions collected as the rate of 
ionization is changed. 

The measurements made by Bowen were 
undertaken primarily to elucidate the behavior 
of cosmic-ray electrometers which use gases at 
high pressure. Consequently great care was taken 
to eliminate volume recombination. This was 
done by the use of a very weak source of radiation 
which produced a fairly uniform ionization of 
only 121 ions per cc per sec. per atmosphere. 
Furthermore at the two highest pressures used 
tests for the presence of volume recombination 
were made by taking a complete set of additional 
measurements when the ionization was reduced 
to 23.4 ions per cc per sec. per atmosphere. 
Similar precautions for the elimination of volume 
recombination were taken by Cox‘ in his studies 
of ionization in other gases than air. 

On the other hand Erikson, who did not at- 
tempt to eliminate volume recombination, used 
ionization rates of from 2 X 10° ions per cc per sec. 
per atmosphere at points of the ionization 
chamber farthest from the radium to 10’ ions 
for points nearest it. Clay and Van Tijn, accord- 
ing to a letter to one of the authors correcting 
their Table I, used an average of 2.710‘ ions 
per cc per sec. per atmosphere. 

Neither Erikson nor Clay and Van Tijn tested 
for the presence of volume recombination by 
varying the average ionization used. Fortunately 
their collecting field and pressure ranges over- 
lapped those of Bowen and consequently such a 
test can be made by superimposing their curves 


& E. F. Cox, Phys. Rev. 45, 503 (1934). 
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on Bowen’s curves obtained at low rates of 
ionization. This is done in Fig. 1 where the 
number of ions collected is plotted against the 
logarithm of the collecting field for the pressures 
nearest 95 atmospheres used by each observer. 
All of the curves are reduced to the same scale by 
expressing the number of ions collected as the 
fraction of the number collected at 1000 volts 
per cm. Comparing first Bowen’s curves at 23.4 
ions and 121 ions it is seen that lack of super- 
position, indicating the presence of volume re- 
combination at the higher ionization rate, sets in 
at a collecting field of 10 to 15 volts per cm. 
Likewise the curves obtained by the other ob- 
servers using high rates of ionization definitely 
drop below the low ionization curves at 400 to 
800 volts per cm, thus indicating the presence of 
volume recombination at all collecting fields 
below this. Similar results are obtained at other 
pressures although the collecting field at which 
volume recombination becomes negligible varies 
with the pressure. 

The results of Erikson and Clay and Van Tijn 
are therefore not directly applicable to the prob- 
lem of cosmic-ray ionization chambers which 
normally operate with collecting fields far below 
400 volts per cm. Furthermore one is not justified 
in applying the Jaffé formula for columnar re- 
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Fic. 1. Comparison of ionizations as a function of collecting 
field as obtained by various observers. 
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Fic. 2. Application of the Jaffé theory of columnar ioniza- 
tion using the graphical method suggested by Zanstra. 


combination except to points that fall in the 
range which is free from volume recombination, 
i.e., in Erikson’s data to points from 513 to 1223 
volts per cm at 140 atmospheres and from 365 to 
1288 volts per cm at 40 atmospheres, and in Clay 
and Van Tijn’s data from 817 to 4917 volts per 
cm. On the other hand Bowen's ionization cur- 
rents obtained at an ionization rate of 23.4 ions 
per cc per sec. per atmosphere are quite certainly 
free from volume recombination over the range 
from 6.2 to 1009 volts per cm. Fig. 2 represents 
the application of the Jaffé theory to these results 
in the above ranges. Following Zanstra the 
abscissa is his f(x) while the ordinate is propor- 
tional to the reciprocal of the ion current col- 
lected. As in Fig. 1 the unit of current is taken as 
the current collected with a field of 1000 volts per 
cm. When plotted in this way all points obtained 
at a given pressure should fall on a straight line 
according to the Jaffé theory. The curves are 
apparently linear over the small ranges of the 
collecting fields available in Erikson’s data where 
the field varies by a factor of about 3 and in Clay 
and Van Tijn’s data where the range is a factor 
of 6. On the other hand the linearity breaks down 
completely in the range in the field of a factor of 
160 available in Bowen’s data. This breakdown is 
still more striking if Clay and Van Tijn’s points 
are combined with those of Bowen at the same 
pressure thus extending the range to a factor of 
800 in the collecting field. For values of f(x) 
greater than 5, f(x) may be represented by the 
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equation f(x) =A —4.6 log (collecting field/pres- 
sure) with an error of less than 1 percent. Conse- 
quently it is evident that the linearity cannot be 
improved by any reasonable shift of the constant 
A which is the only quantity not fully fixed by 
the Jaffé theory. 

Another difficulty with the interpretation of 
these curves in terms of the Jaffé formula appears 
even in its application to the short linear sections 
of Clay and Van Tijn’s curves. Thus they found 
from extrapolation by means of the Jaffé formula 
that the saturation currents, i.e., the number of 
ions originally formed, were less than half as 
great per atmosphere at 95 atmospheres as at one 
atmosphere. To explain this they assumed that at 
the lower pressure most of the secondary 6-par- 
ticles, that are the direct cause of the ionization, 
are ejected by the primary y-rays from the walls 
of the vessel while at the higher pressure a large 
part of the 8-rays are ejected from the air itself. 
The lower ionization per atmosphere at the 
higher pressure is then to be expected if the 
efficiency of the ejection of these secondary 
8-particles is greater from the walls, which are of 
relatively high atomic number, than from the air 
of low atomic number. However, Workman’? 
found that the ionization was actually increased 
when the particles were ejected from an element 
of low atomic number such as carbon or alu- 
minum rather than from steel or lead. Thus, if 


7E. J. Workman, Phys. Rev. 43, 859 (1933). 
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anything, more ions per atmosphere should be 
formed at the higher pressures. 

All of these considerations show, therefore, 
that the observed variation of ion current with 
collecting field cannot be explained by the Jaffé 
formula for columnar recombination. 

In conclusion the authors wish to emphasize 
again the importance of using very low rates of 
ionization in measurements at high pressure in 
order to eliminate ordinary volume recombina- 
tion. When volume recombination is present the 
current collected is not proportional to the in- 
tensity of the ionizing radiation. Furthermore the 
curves of ion current vs. collecting field become a 
function not only of the intensity of ionization as 
discussed above but of the dimensions of the 
collecting chamber as well. Consequently the 
significance of these curves is limited chiefly to 
that of a calibration of the ionization chamber 
used. The limiting ionization current that can be 
used before volume recombination becomes ap- 
preciable depends of course on the pressure, 
collecting field and dimensions of the apparatus 
used. In general, however, for pressures over 25 
atmospheres and collecting fields below a few 
hundred volts per cm the ionization current 
should be made as low as is feasible without 
having it masked by the residual ionization 
caused by cosmic rays and local radioactivity. 
Finally, measurements should be made at differ- 
ent rates of ionization to determine whether 
volume recombination has really been eliminated. 
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Depolarization of Raman Lines and Structure of Chlorate, Bromate and Iodate Ions 


S. T. SHEN, Y. T. YAo anp Ta-You Wt 
Department of Physics, National University of Peking, Peking, China 


The Raman spectra of chlorate ClO;, bromate BrO; and iodate IO; in solutions are examined 


under high dispersion and the degrees of depolarization of the Raman lines measured. It is 
shown that the strong line at 930 cm™ of ClO; and that at 800 cm~ of IO; in solutions, which 
are regarded as single lines in the literature, are each really composed of two lines. The presence 
of four lines in each ion in the Raman effect, together with the depolarization measurement, and 
the infrared data, shows that these ions have pyramidal structure. The valence force constants 


of these ions are calculated from the present data according to the equations of Lechner. 


INTRODUCTION 


N the literature! it is assumed that the chlorate 
ClO;, bromate BrQO; and iodate IO; ions have 
symmetric, plane structure like the carbonate ion 
CO;. Asymmetric, plane structure belongs to the 
symmetry class D3, and possesses four vibrations: 
one totally symmetric vibration which is active 
in Raman effect but inactive in infrared ab- 
sorption; one asymmetric with respect to the 
plane of the molecule and active in the infrared 
but inactive in Raman effect; and two doubly 
degenerate vibrations which are active in both 
Raman effect and the infrared. An examination 
of the Raman spectrum of CIO; ions in solution 
revealed, however, that in addition to the two 
lines at 479 cm= and 610 cm™, the strong line 
at 930 cm has a wing on the high frequency 
side (measured from the exciting line), which 
suggests the presence of a weak line at 982 cm™. 
It is of some interest to investigate this further ; 
for the presence of four Raman line would point 
to a pyramidal structure for the ion rather than a 
plane structure. A symmetric pyramidal model 
belongs to the symmetry class C3, and possesses 
two totally symmetric vibrations 7, v2 in which 
the change of the electric moment is along the 
axis of the pyramid, and two doubly degenerate 
vibrations v3, v4 in which the change of the 
electric moment is perpendicular to the axis. All 
four vibrations are active in Raman effect. As 
1C. Schaefer, F. Matossi and H. Aderhold, Zeits. f. 
Physik 65, 289 (1930); C. Schaefer and F. Matossi, Das 
Ultrarote Spektrum; Cabannes, Ann. d. Physik 18, 285 
(1932), remarked in passing that in contrast with COs, 
ClO; has four Raman active frequencies but v3 could not 
be observed. We are indebted to Dr. H. C. Cheng for 
calling our attention to this remark. From x-ray studies 
and wave mechanical considerations, however, it is known 


that these ions ClO;, BrO;, IO; have pyramidal structure. 
See, for example, Slater, Phys. Rev. 38, 325 (1931). 


the line 930 cm= of ClO; cannot be resolved 
from the wing even with high dispersion (12A 
per mm), and as the corresponding line in IO; 
(800 cm“) is only barely resolved into two lines 
at 779 and 826 cm™, it is thought desirable to 
study the depolarization of the Raman lines of 
ClO;, BrO; and IO; with a view (i) to decide 
whether the line 930 cm= in ClO; and its wing 
are really two lines by observing their difference 
in the degrees of depolarization, and (ii) to 
decide between the plane model and the pyrami- 
dal model by comparing the observed spectra and 
the depolarizations with the selection and de- 
polarization rules for the two models. 


EXPERIMENTAL 


The experimental method is of the conventional 
type. A mercury arc 14 cm long irradiated the 
Raman tube at a distance of 2 cm through two 
slits, each 0.8 cm wide, one near the arc and the 
other near the Raman tube The scattered 
radiation was focused by a lens one meter from 
the tube and 43 cm from the slit of a three-prism 
Steinheil spectrograph with glass optics. Arrange- 
ment for exact parallelism of the exciting radia- 
tion was not attempted for it would necessitate 
unduly long exposures. The radiation was 
analyzed by a method due to Hanle.? A thick 
calcite crystal decomposed the radiation into two 
components; one with electric vector parallel to 
that of the exciting radiation and perpendicular 
to the length of the slit of the spectrograph, and a 
depolarized component with electric vector per- 
pendicular to that of the exciting radiation and 
parallel to the length of the spectrograph. A 
second thin, wedge-shaped calcite crystal, cut 


2 W. Hanle, Zeits. f. Instrumentenk. 51, 488 (1931). 
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Raman spectra (a) without polarization, (b) polarized component, (c) depolarized component. 


ClO; 
10 N. NaClO; 
0.06 mm 
250 hr. 


Concentration 
Slit width 
Exposure 


with the optical axis in one face, was set with the 
axis at an angle of 45° with the slit to produce 
two images of depolarized light on the slit, so 
that difference in the loss of light due to reflec- 
tions from the prism surfaces of the two com- 
ponents was eliminated. With this disposition, 
the two components were photographed simul- 
taneously and exposures could be prolonged for 
the weak lines. 

The comparison of the intensities of the two 
components of the Raman lines was made by the 
method of Thomson and Duffendack.* Intensity 
marks were put on the plate by means of a step- 
diaphragm. Intensity measurements were made 
on the microphotometer tracings of the line and 
the blackening-intensity curve at the position of 
the line. 

In this study, excitation of the spectra was by 
means of the 4358 line. The jacket of the double- 
walled Raman tube contained a solution of 
sodium nitrite to filter off the line 4047 so as to 
eliminate the water band excited by it. This is 
necessary, for the broad water band falls in the 
region of the Raman lines of the ions excited by 
the line 4358 and would make intensity measure- 
ments very difficult. The Raman spectra of 
ClO;, BrO; and IO; are represented in Figs. 
24. 

From Fig. 1, it is seen that the wing is more 
strongly depolarized than the line at 930 cm“, 
and it can be reasonably ascribed to a weak line 

*K. B. Thomson and O. S. Duffendack, J. Opt. Soc. Am. 
23, 101 (1933). 


BrO; 
3 N. NaBrO; 
0.20 mm 
230 hr. 


10; 
6.1 N. HIO; 
0.06 mm 
210 hr. 


with center at 980 cm~, though a decomposition 
of the intensity curve into two overlapping lines 
cannot be done entirely without arbitrariness. 
In BrOs, the line at 421 cm™ is observed for 
BrQ3; in solution for the first time. The strong 
line at 806 cm™ is quite symmetrical, showing 
that the frequency v3 corresponding to the wing 
of the line 930 cm™ in CIO; either is very close to 
v,=806 cm™ or is too weak to appear since the 
concentration of BrOx; is low. In obtaining the 
curves } and ¢, the slit was opened up to 0.2 mm 
and the exposure was 210 hr. Even then the two 
components of v2 at 421 cm™! were too weak to 
give an accurate value for the degree of depolari- 
zation. In IQO;,‘ the line at 330 cm~ is broad and 
is due to two lines v4, ve very close together. The 
line at 800 cm™ is seen to be composed of two 
lines v; and v3, the latter being more strongly 
depolarized than the former. Because of the 
overlapping of the two lines, an exact measure- 
ment of the intensity of the weak line »; is 
impossible. However, by assuming a reasonable 
shape of the line »; and by assuming that the 
half-width of the line is the same in the polarized 
and the depolarized components, a decomposition 
of the line into »; and v3 may be performed 


4 The line at 639 cm™ is present in the Raman spectra of 
HIO; crystals but is not found in KIO; crystals (Kohl- 
rausch, Smekal-Raman Effekt, S. 351). By varying the con- 
centration of the HIO; in the solution, we found that the 
intensity of the line 639 cm™ relative to the line 800 cm™ 
decreases with decrease of concentration. These facts 
seem to show that the line 639 cm™ is due to the undis- 
sociated molecule HIO;. Quantitative measurement on the 
variation of the relative intensities of the line 639, 800, 
330 with the concentration will be published elsewhere. 
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STRUCTURE OF ClO;, 


graphically. The degree of depolarization was 
then calculated by means of the log blackening- 
intensity curve. The values of p’s for v3 of IO; and 
ClO; obtained in this manner must be regarded as 
qualitative. 

The results of this study are summarized in 
Table I, which includes the infrared data® for 
comparison. 


DISCUSSION OF RESULTS 


The presence of four lines in the Raman 
spectra of ClO;, BrO; and IO; shows that they 
are inconsistent with a plane model, but fit in 
quite well with a pyramidal model. According 
to the polarization rule of Placzek,® the two 
symmetric vibrations »; and v2 are partially 
depolarized and the values of the degrees of 
depolarization p should be less than 6/7; while 
the two degenerate vibrations v3 and vq are 
completely depolarized and p=6/7 for natural 
light. It is seen from Table I that wherever p has 
been measured, it agrees with the theory, except 
for v3 of IO; in which case p cannot be determined 
accurately because of the closeness of v3 to v1. 

Further evidence’ for the pyramidal structure 
is furnished by the infrared absorption curves of 
these ions in crystals. For the ClO; ion, Schaefer, 
Schubert, Matossi and Laski® observed three 
bands at 23u (434 cm™'), 16u (624 cm) and 

§C. Schaefer and M. Schubert, Zeits. f. Physik 7, 309 
(1921); C. Schaefer and Matossi, Das Ultrarote Spektrum, 
S. 328. 

6 Placzek, Leipziger Vortrdge, 1931. 

7 An indirect reason against the plane model for these 
ions is that even if one neglects the presence of the line » 
in the Raman effect and proceeds to calculate the force 
constants of the ions, on the assumption of a plane model, 
using either the central force system of Dennison-Nielsen, 
Phys. Rev. 32, 773 (1928), or the valence force system of 
Lechner, one is led to inconsistency in the equations giving 
the frequencies. 
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curve of KCIO;. 


10.04 (994 cm). At low dispersion, the band at 
10u shows an asymmetry consisting of a wing on 
the long wave-length side. A comparison between 
the infrared curve and the Raman line drawn on 
the same scale is given in Fig. 4. It is seen that 
the wing in the infrared absorption curve corre- 
sponds to a weak band due to the parallel 
vibration »; which appears strongly in Raman 
effect, and the wing in the Raman curve corre- 
sponds to the perpendicular vibration v3 which is 
strong in infrared absorption. 

For BrO; two infrared bands at 234 (434 cm~) 
and 12.24 (820 cm~') are known, and a third at 
about 36y (280 cm~') is obtained by extrapola- 
tion. For IO3, a band at 12.64% (793 cm™) is 
observed, and one at 27y (357 cm™~) correspond- 
ing to ve was obtained by extrapolation. A 
comparison between the infrared curve for the 
12.64 band and the Raman line at 800 cm™ is 
given in Fig. 5. 

Considering the low dispersion with which the 
infrared absorption curves were obtained, and 


®C. Schaefer and F. Matossi, Das Ultrarote Spektrum, 
S. 328 


TABLE I. The infrared and the Raman frequencies of ClO;, BrOs and IO3. The values in brackets are obtained by extrapolation 
only.> The values marked with an asterisk are taken from crystal data, Kohlrausch, ‘‘Smekal-Raman Effekt.”’ 


mill well va(2) & (2) 4 
RAMAN INFRARED RAMAN INFRARED RAMAN INFRARED RAMAN INFRARED 

ClO; »v 930 cm™ 950 cm™ 610 cm™ 624 cm™ 982 cm™ 994 cm™ 479 cm” 434 cm™ 

p 0.35 0.56 0.87 0.90 
BrO; »v 806 800 421 434 836* 820 356 278) 

p 0.38 0.89 
IO; »v 779 780 390* (357) 826 793 330 

p 0.43 0.67 0.81 








238 ZUMSTEIN, GABEL AND McKAY 
mM TABLE II. Constants of the CLO;, BrO3, IO; tons. 
fx<1075 6d X10 B 
ClO 5.55 6.42 54 
BrO 5.25 3.52 53.5 
10; 5.35 3.19 53.3 
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remembering that the long wave-length bands® 
were obtained by the Reststrahlen method of 
which the accuracy is very low, and that the 
infrared frequencies and the Raman frequencies 
were obtained for the ions in crystals and in 
solutions respectively, one feels that the agree- 
ment between the Raman data and the infrared 
data is as good as one would expect. It is de- 
sirable to examine the 10u band of ClO; and the 
12.6% band of IO; in the infrared with higher 
dispersion, and to determine more accurately the 
positions of the long wave-length bands of these 
ions. 

Assuming a pyramidal model, one can calcu- 
late the force constants of these ions. On the 
basis of a valence force system, Lechner® ob- 
tained the following equations for the four 


9F, Lechner, Graz Dissertation, 1933; Monatshefte f. 
Chemie 61, 385 (1932). 
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frequencies v;=+/A;, 27: 
AiA2= (f/m) (6d/m)2p cos* B/ (1+3 cos? B), 
Ait d.=(f/m)[1+(p—1) cos? B | 
+ (6d/m)[p—(p—1) cos* 8 ]2 cos? B 
(1+3 cos? Bg), 
4dsk4= (f/m) (6d/m)[p+(2—>p) cos* 8 ]2 
(1+3 cos? B), 
2(As+A4) = (f/m) [p+1—(p—1) cos? 8] 
+ (6d/m)[1+(p—1)(1—cos? 8)? 
2(1+ cos? 8) ](1+cos* 8) /(1+3 cos* 8), 


where p= (3m+M)/M, and 1, ve are the single 
frequencies and v3, v4, are the doubly degenerate 
frequencies. 

By substituting the observed frequencies in 
these equations and so adjusting the values of 
the constants f, d and 6 that the two sides of the 
equations are closest to each other, we obtain the 
values of f, d, 8 for the ions given n Table II. 

In conclusion, we wish to acknowledge our 
indebtedness to Dr. Ny Tzi-Zé, of the Institute of 
Physics, Peiping Academy, for the use of the 
microphotometer. 
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The '>*—'> Band System of Ionized Cadmium Deuteride 
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The '5*—z bands have been photographed with a dispersion of about 1.3A/mm. Eight 
bands of the system, between \2630 and \2280 were analyzed. The isotopic coefficient p agrees 
with the value obtained from atomic weights. The electronic isotope displacement relative to 
Cd*H appears to be 4.4 cm~!. However, this value is somewhat uncertain as only three vibra- 
tional levels of the upper state and four of the lower state are available. The lines did not show a 


fine structure from the cadmium isotopes. 


HE band spectrum of the ionized cadmium 
hydride molecule has been investigated by 
Svensson and Tyrén' who give references to 
previous work. We have studied the similar 
spectrum obtained by substituting deuterium for 


(1933). 


hydrogen, thus providing additional experi- 
mental evidence concerning the problem of the 


electronic isotope effect. 
EXPERIMENTAL 


The spectrum was excited in an end-on dis- 
charge tube of fused quartz. The length of the 
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TABLE I. Wave numbers of the ionized cadmium deuteride bands. 


J R J R I J RK I J R P J R I 
0 —0 band 0—1 ba n 0—2 band n 1—2 ban ” ban 
4 | 42758.4* 42715.4 || 15 | 472.6 | 324.5) 24) 178.3 | 947.5|13| 193.9 | 064.1) 3 42068.5 
5} 758.4*| 704.0//16| 461.1 | 302.9|/ 25 159.3 14/ 184.4 | 046.6) 4 059.6 
6} 755.1 | 691.8 || 17} 447.2 | 280.2/|26| 140.7 15| 174.8 | 028.0) 5 049.3 
7| 751.5 | 678.8/||18| 433.3 | 256.9||27| 120.1 16| 164.3 |41007.0 || 6 038.4 
8| 747.0 | 664.0|}19| 417.9 | 231.5 ||— 17, 149.1 |40985.7 || 7 026.6 
9| 741.0 | 648.8 || 20) 401.3 | 206.1 1—0 band 18| 139.3 | 964.9) 8 $2013.4 
10| 734.2 | 632.1 || 21) 383.3 | 180.1 19; 121.3 | 941.2) 9 41999.2 
11 726.1 615.0 || 22 363.9 152.2 4  43636.1* 20 109.6 916.4 || 10 984.2 
12} 716.5 | 595.6||23| 344.5 | 121.3|| 5! 636.1*|43582.7| 21 | 093.8 | 896.2 )| 11 968.2 
13| 705.6 | 575.5||24| 324.5 | 093.8|| 6| 632.9 | 570.5| 22) 077.0 12 950.6 
14; 693.6 | 554.7||25| 302.9 | 064.1 || 7| 6285 | 556.1) 23) 058.5 13 932.2 
15| 680.9 | 531.726! 280.2 | 032.2|| 8! 623.4 | 541.8/) 24) 041.2 14 913.2 
16| 666.5 | 5084/27; 256.9 | 000.1) 9| 6174 | 526.2 15 893.6 
17| 650.3 | 483.7 || 28 40964.9 || 10 609.8 | 509.0 16 872.2 
18| 633.8 | 457.7 || 29 930.3 || 11| 600.6 | 490.6 i—3 band 17 849.7 
19} 615.9 | 430.5 12) 590.3 | 470.9 18 826.9 
20} 596.4 | 401.8 ee 13} 578.9 | 449.9)) 4 40062.3 
21| 575.5 | 373.1 o~3 ane 14| 566.1 | 425.5|| 9 056.3 ; 
22| 554.8 | 342.2 ——|}15| 552.1 | 405.0|| 3 0475 of 
23 531.7 310.2 2 40322.6 || 16 537.0 380.1 4 039.8 
24 508.4 277.5 3 314.6 || 17 521.1 354.6 || 5 030.3 2 39818.0 
25 | 483.6 | 244.0)| 4 305.8 || ! 502.1 327.3 || 6 | 40083.1*| 020.7) 3 811 
26 209.3 || 5 | 40349.9*/ 296.8 || 19 | 483.2 | 299.9|| 7) 083.1*| 009.7 || 4 803.3 
27 172.7 || 6| 349.9*| 286.2 || 20) 462.4 269.9 || § | 081.2 |39998.7 |) 5 794.9 
28 135.0 || 7) 349.9*| 276.3 || 21 441.7 239.3 ||} 9) 077.1 986.3 || 6 | 39846.1*| 784.9 
= 8| 345.5 | 263.5 || 22) 417.6 | 2084) 10) 073.7 | 973.0|| 7| 846.1*| 775.4 
eat tine 9) 342.1 | 250.0) 23) 393.9 | 175.1) 11 | 068.4 | 959.1 || 8| 846.1*| 764.0 
10} 337.2 | 235.9) 24) 369.5 | 140.7)/12|) 063.4 | 944.5) 9| 842.9 | 753.0 
™ 11| 330.0 | 221.4)| 25 105.4 |} 13| 057.3 | 928.8)/|10| 839.0 | 740.2 
3 41503.0 || 12} 326.3 | 205.3 —\} 14} 050.2 | 912.5) 11} 8364 | 730.1 
4|41535.3*| 491.9|/13| 319.2 | 188.5 1—2 band 15| 042.0 | 894.1|/12| 831.3 | 714.1 
5} 535.3*| 481.6||14|) 311.2 | 171.5 16| 033.7 | 876.1 | 13| 825.8 | 699.5 
6| 533.7 | 470.6||15|} 301.9 | 152.6|| 4 | 41231.5*|41186.4|/17| 024.2 | 858.2|14| 819.8 | 684.5 
7| 531.1 | 458.8 |} 16) 291.7 | 132.0]| 5| 231.5*| 177.4]/18] 012.9 | 839.0//15| 813.2 | 668.4 
8| 527.5 | 444.8//17] 281.1 113.7 || 6| 2294 | 165.8 || 19|40001.7 | 818.0 || 16| 804.7 | 651.2 
9} 522.8 | 430.5 ||18| 269.1 | 091.9|) 7| 227.5 | 154.7 || 20|39989.1 | 797.6|17| 797.6 | 634.4 
10} 517.4 | 415.3)|}19| 256.8 | 070.2|| 8| 224.3 | 143.0|/21| 976.0 | 775.4//18| 786.9 | 615.6 
11} 510.7 | 399.6|| 20} 242.6 | 047.5|| 9| 220.6 | 1291/22! 962.8 | 753.0)|19| 777.6 | 596.7 
12} 502.7 | 382.2||21| 228.3 |40024.2||/10| 215.2 | 114.9||23| 947.5 | 730.1||20| 765.0 
13| 494.0 | 363.9||22| 212.0 |39998.7|/11| 209.5 | 098.9|| 24] 932.2 | 704.4/|21| 755.0 
14) 484.0 | 344.5||23| 195.9 | 973.0//12| 200.8 | 082.5/|25| 915.9 22 | 742.0 


* Band heads. 


tube was 20 cm and its diameter 1 cm. Small 
pieces of cadmium metal were scattered along 
the fused quartz tube and the latter placed 
inside a water cooling jacket made of Pyrex 
glass. Since the water did not strike the quartz 
directly we obtained sufficient temperature to 
evaporate the cadmium. The heavy hydrogen, 
prepared by the reaction between calcium metal 
and 99.5 percent deuterium oxide, slowly entered 
the tube through a needle valve. The spectrum 
was excited by a 2300 volt direct current gen- 
erator with sufficient series resistance to limit 
the current to about 2 amperes. After a few 
minutes operation we always observed a rapid 


clean up of the deuterium. This seemed to be 
very important when a spectrum free from the 
light hydrogen is desired. By this method plates 
were obtained showing no trace of the light 
hydrogen. The spectrum was photographed with 
a large Hilger quartz spectrograph having a 
dispersion of about 1.3A per mm. This is about 
the same as the dispersion of a 30,000 lines/inch, 
21-foot grating when used in the first order. 


ANALYSIS 


The measurements of the individual wave 
numbers of the band lines are given in Table I 
directly as measured with a Société Génevoise 
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ase Il. Null lines raB_e V. Isotope ratios Cd*D : Cd*H 
0 1 2 3 } = 7 i 
ar : i Be‘/Be)! 0.7117 0.7108 
O | 42746.9 | 41523.2 | 40334.8 we! wr 0.7113 0.7086 
1 43626.6 41215.7 | 40068.6 
2 42089.5 39829.6 
u/ ui)? 0.7105 
TABLE III. Constants from the rotational structure. 
TABLE VI. Other band heads observed but not analyzed. 
B ry’ | my —D,”’ = — 
; — - > Designation v 
0 2.438 0.000076 3.052 0.000058 . — — —~ 
l 2.412 0.000084 2.966 0.000060 4 1-4 38976.9 
2 2.380 0.000092 2.903 0.000079 2-5 780.0 
3 2.833 0.000086 3-6 616.1 
4 2.777 0.000108 4-7 482.6 
5-8 382.4 
6-9 315.7 
7 72 ? 
TABLE IV. Summary of constants > o904 
9-12 261.2 
] wer st ite upper state 
i 75 ” . ° 
a oo a? In this analysis we have used the energy ex- 
s 0 42930.6 pressions given by Mulliken.* Wherever possible 
B 3.075 2.451 ant anes . : ‘ ‘hans 
4 as squeé — ; ( ( ~ « S *¢ 
oi 00682 0.027 the least squares meth od of Birge ind he i 
D, —0.000048 —0,000072 was used in performing the calculations. The 


comparator. The band heads are indicated by 
starred lines in the R branch. Within any band 
the measurements should have an average 
accuracy of about 0.4 cm~', In addition to these 
errors of setting the microscope, the band as a 
whole may be shifted relative to the comparison 
spectrum by about 0.5 cm~. In the 0—Oand 0—1 
bands, which are favorably situated, the return- 
ing R branch almost coincides with the P branch. 
This entirely hides the null line and other faint 
lines nearby. Fortunately the 0—2,1—3 and 2—4 
sequence presented no such overlapping. The 
null lines were close to the band heads. The con- 
tinuous spectrum of the deuterium molecule was 
always present and this also helped to hide faint 
lines. Thus the assignment of J values de- 
pended largely on the combination differences, 


R(J)—P(J) and R(J—1)—P(J+1). 


B and D constants were also determined by the 
graphical method described by Bengtsson Knave.* 
It was always found that the individual lines of 
the P branch could be represented by a least 
squares parabola within the limit of error of 
the observations. The results of the analysis are 
given in Tables II, III and IV. The system as a 
whole checks with the established isotope theory 
as is shown by Table V. If we take 7, for 
Cd*H as 42,935 we find an electronic isotope 
displacement for Cd*+D of 4.4 cm. It seems to 
us that more bands of the system should be 
analyzed with greater dispersion before this last 
number can be regarded as certain. A fine 
structure of the lines due to the cadmium isotopes 
was not observed. 


?R. S. Mulliken, Rev. Mod. Phys. 2, 66 (1930). 

3R. T. Birge and J. D. Shea, University of California 
Pub. in Math. 2, 93 (1927). 

4 E. Bengtsson Knave, Nova. Acta. Reg. Soc. Upsaliensis 


Ser. IV, Vol. 8, No. 4, p. 12 (1932). 
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Interferometer determinations of the wave-lengths of thirty-eight lines in the spectra of 
deuterium and HD in the region \4860 to 5800 have been made by a crossed etalon and 
Littrow spectrograph. The measurements show a mean probable error of 0.0003A. 


INTRODUCTION 


ITH a spectrum as rich in lines as the 
secondary spectrum of hydrogen, inter- 
ferometric determinations of the wave-lengths of 
certain lines distributed throughout the spectrum 
are extremely valuable as tertiary standards. 
The present paper covers the region from 
44860 to A5800 with De and HD. 


APPARATUS 


The Pyrex glass discharge tube was built in U 
form with platinum terminals at each end, the 
capillary being of 7 mm internal diameter and 36 
cm length, coated inside with platinum bright to 
enhance the molecular spectrum. The tube was 
mounted in a tank in which a continuous flow of 
water was maintained. A 5 kw transformer sup- 
plied to the terminals of the tube a 190 ma cur- 
rent (current density 490 ma/cm*) at 14,000 volts. 
A flask of deuterium gas' was connected with the 
vacuum system through three stopcocks, thus 
rendering the gas flow controllable. 

The general arrangement of apparatus is shown 
in Fig. 1. 

The auxiliary dispersing Littrow spectrograph 
has a lens of 6-inch aperture and 30-foot focal 
length with an Anderson plane grating of 15,000 
lines to the inch and a ruled space of 9.7 12.9 
cm. The first order was used, the dispersion being 
about 1.8A per mm. The etalon was aluminum 
coated and had a separation of about 10 mm. 

When the spectrograph slit was open, the 
image of a diaphragm (not shown in Fig. 1), 
placed as near the etalon as possible and on the 
side away from the source, formed upon the 


! Purchased from Professor Hugh S. Taylor of Princeton 
University. 


grating by the concave mirror N, nearly covered 
the ruled surface of the grating. 


PROCEDURE 


The spectrum of a standard Pfund are was 
photographed on each plate, the arc being ex- 
posed for about one minute every half-hour dur- 
ing the tube exposure which in all cases was nine 
hours. The wave-lengths recommended by the 
International Astronomical Union were used as 
standards.’ 

The Dz spectrum was photographed first by 
allowing a small quantity of gas to pass from the 
reservoir into the tube and then isolating the 
tube by a stopcock at each end. A pressure of 0.3 
mm was found to produce the best spectrum. 


G Pp’ 
Heons’p \. 
L , D\ b] 7, 
=. « fi.f 
XX a fh nl 
- tY - Nu 
‘oD 
C+ 
“A 


Fic. 1. General arrangement of optical train. Not drawn 


to scale. 
G—plane grating 
L—six-inch lens of 30-foot focal length. 
S—slit of Littrow spectrograph. 
P—photographic plate (below slit 
N—concave mirror of 97.3 cm focal length 
E—etalon 
D—diaphragm with opening approximately 1 cm square. 
N’, N’’—plane mirrors 
L’—lens forming image of capillary of tube T upon D when N’” is 


removed 

N’’’—plane mirror used for throwing light from arc upon D, light from 
tube being then intercepted 

7—vacuum tube. 

N’’"’—plane mirror to increase brilliancy of source 

D’—diaphragm with 1 cm vertical aperture to cut out all but central 
portion of arc 

L’"—short focus lens giving 10-fold magnification of source A upon D’ 

A—Pfund arc, 13 to 15 mm in length. 

All mirrors are aluminum coated 

All lenses are achromats 

Mirrors N’ and N”’ were necessary because of the narrowness of the 
room. 


? Transactions of International Astronomical Union III, 
86 (1928). 
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raBLE I, Wave-lengths, in angstroms, of deuterium and HD. 


Deuterium HD 
Wave- Wave Wave- 
Intensity length Intensity lengt! Intensity lengt 
0 4875.671 0! 5460.467 0 4900.298 
08 4913.497 0 5473.457 0 4988.770 
03 4963.200 0 5481.942 1 4990.692 


0 4965.643 0! 5488.596 0 4991.654 
0+ 4968 .567 l 5609.285 0! 5209.520 
08 5024.446 0 5631.729 0 5362.868 
0? 5070.930 0 5639.085 0 5381.903 
m 5100.014 0 5664.392 1 5477.495 
08 5204.451 1 5690.228 0 $525.040 
0+ §254.361 1 5702.698 1 5563.479 
04 5259.619 2 §715.583 1 5582.449 
03 5294.952 1 5796.627 0 5589.146 

0 5624.862 

0! 5667.768 


To obtain the HD spectrum, hydrogen, elec- 
trolytically prepared and dried by passing 
through tubes of calcium chloride and dessi- 
chlora, was passed into the tube containing 
deuterium and the pressure pumped down to the 
correct value. The spectrum was observed visu- 
ally by means of a direct vision spectroscope, this 
serving as a guide to the appropriate mixture of 
hydrogen and deuterium. Eastman spectroscopic 
plates I-D were used throughout. 


RESULTS 


The method of reduction of wave-lengths was 
that employed at Mount Wilson,’ and was 
based on the familiar relations; 2¢=pd, and 
E=D?/(D,2—D,_;*) where ¢ is the thickness of 
the etalon; p the order of interference at the 
center of the rings; \ the wave-length; E the 
fractional part of the order of interference; D the 
diameter of a bright ring; and k the number of 
the ring, counting from the center. 


~ 8 Babcock and St. John, Astrophys. J. 53, 260 (1921). 


R. G. LACOUNT 

Dr. Dieke of Johns Hopkins University kindly 
provided most of the preliminary wave-lengths 
forming the basis of the measurements. These 
were correct to 0.01A; hence only one thickness 
of etalon was necessary. 

The wave-lengths listed in Table I are the 
means of measurements on five plates for Dz and 
five for HD. For only two of the twenty-five D, 
lines measured is the probable error as great as 
0.0006, and for sixteen it is not greater than 
0.0003, the average being 0.0003A. Of the four- 
teen HD lines the probable error for one is 0.0005, 
and for ten it is 0.0003 or less, the average being 
0.0003A. 

The intensities listed were estimated by eye, 
O* indicating that the line is just seen on the 
plate, O being the average intensity and 3 that of 
the strongest lines. 

The major part of the work was done by the 
junior writer of this paper and was carried out in 
the Rumford Committee room in the Eastman 
Spectroscopic Laboratory of the Massachusetts 
Institute of Technology. This entire building is 
insulated from vibration and thermostatically 
controlled, and is therefore well suited for long 
exposures. 

The writers wish to acknowledge their in- 
debtedness to the Institute for the splendid 
facilities placed at their disposal, and to Drs. 
Harrison and Boyce especially for their valuable 
suggestions. 

We are also indebted to Mr. Babcock of Mt. 
Wilson Observatory for his helpful advice, to 
Professor Wood of Johns Hopkins for the loan of 
the grating, and to Dr. Dieke for furnishing some 
of the necessary preliminary wave-lengths. 

A grant from the Rumford Committee of the 
American Academy of Arts and Sciences to the 
senior writer made this research possible. 
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ly The Zeeman Effect of N II 
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J. B. GREEN AND H. N. MAXWELL"* 
oe i Mendenhall Laboratory of Physics, Ohio State University, Columbus, Ohio 
‘SS Received November 26, 1936 
he / Che Zeeman effect of N II has been studied at field strengths of about 36,000 gauss. By using 
id helium as residual gas with a small nitrogen impurity, the spectrum has been strongly excited, 
Do especially in the visible region. Edlén’s corrections to Freeman's classification have been 
as 4 verified. Except in a few important cases, g values are LS within experimental error. The 
™ exceptions are considered in detail. 
r- 
5, 
1g HE spectrum of N IT was first well classified copper electrodes. The uppet electrode was a 
by Fowler and Freeman! and later extended — piece of No. 20 gold wire held in a brass rod and 
e by Freeman? who found quintet combinations in the lower electrode was a gold or copper wire 
1e addition to the singlets and triplets and singlet- wrapped on a brass disk which rotated about an 
of triplet intercombinations found by Fowler and axis so that its periphery was at the center of the 
Freeman. In the meantime, some work by Pretty® magnetic field. The interaction of the current in 
1e had disclosed some new singlet levels, but Edlén' the upper electrode and the magnetic field causes 
in has shown that Pretty’s interpretation of his the arc to be intermittent, the greater the current 
in results is probably incorrect. Thus, the 3p 'P, the longer the period during which the arc 
ts and 3p'Dz assigned by Freeman were inter- remains lit. 
is changed by Edlén. The exposures varied from a few hours up to 48 
ly Studies of the Zeeman effect of N II indicate hours, and the arc current averaged about 0.7 
Te that Edlén’s interpretation is the correct one. ampere. Most of the measurements were made 
The Zeeman effect of N II was first studied by from an exposure between gold and copper elec- 
1- Croze® over twenty years ago. The patterns were,’ trodes which lasted about 38 hours and was 
id in general, unresolved, and were therefore of very spread over a period of about four days. The 
s. little use in helping to classify the spectrum or to temperature of the grating was carefully con- 
le study perturbations in the classified spectrum. trolled by an ether thermostat and no effects due 
The present investigation was carried out along to changing temperature could be observed on 
.. with other work in our systematic study of the some of the finely resolved patterns of copper 
oO Zeeman effect. and gold. 
of The atomic spectra of N are not usually About fifty lines in the spectrum of N II were 
le excited with any great intensity in the spark ™easured in a field of about 36,000 gauss and the 
discharge in air or in partial vacuum. But when a results are tabulated in Table I. Typical Zeeman 
le discharge takes place between metal electrodes in P#! eigen are shown in 's. I. 
le an atmosphere of helium with a small percentage rhe lines are arranged in multiplets and the g 


of nitrogen, the spectrum of N IT is very strongly 
excited. This would indicate that the No molecule 
when ionized by electron impact dissociates into 
two N* ions. In the present investigation, the 
discharge took place between gold or gold and 


* Instructor in mathematics, Hood College, Frederick, 
Maryland. 

1 Fowler and Freeman, Proc. Roy. Soc. Al14, 662 (1927). 

2 Freeman, Proc. Roy. Soc. Al24, 654 (1929). 

’ Pretty, Proc. Phys. Soc. London 41, 442 (1929), 

‘Edlén, Nov. Acta Reg. Soc. Sci. Upsala, Sev. IV, 9, 
No. 6. 

Croze, Ann. de physique 1, 35 (1914 


values have been calculated using the resolved 
patterns first, and then using the g values thus 
obtained from the same multiplet (wherever 


| il i 


(a) (b) 


hic. 1. Typical Zeeman patterns of N If. (a) 45045, 
(b) AA5010.6 and 5007.3, (c) \4994.3 
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FABLE I. Zeeman pattern and g values for lines in the 
spectrum of N Tl. 

ZEEMAN 
\ CLASSIFICATION* PATTERN 
3994.995| 3s IP, —3p 'D. 0) 0.980 | 1.052} 1.004 
4432.71 3d 3P, —4f 3D 0) 1.096 | 1.50 | 1.30 
4447.035| 3p'P, —3d "Dz 0) 0.976 | 1.005 | 0.986 
4530.37 | 3d'!F; —4f 1G4? 0) 1.017 
4601.490| 3s 8P, —3p §P; 0) 1.514 | 1.457} 1.495 
4607.167 sp sp 0) 1.530 . 1.530 
4613.884 P - IP 0) 1.500 | 1.465 | 1.535 
4621.405 | >, 3Py 0) 1.457 | 1.457} 
4630.551 3P, 3p 0) 1.502 | 1.509} 1.495 
pogead 3P, 3P, 0) 1.475 | 1.493 | 1.530 
4779.710| 3p 3D, 3d 3— 0) 0.495 | 0.495 | 0.495 
4778.126 D, Ds 0) 1.140 | 1.167} 1.113 
4803.272 D 3p 0) 1.333 | 1.333 | 1.333 
4895.20 | 2p''D, 3p IP 0) 0.971 | 0.985 | 1.013 
4987.377| 3p 3S; 3d 3Po 0) 2.018 | 2.018 
4994.358 - P; |(0.495) 1.505) 2.010} 1.505 
2.009 
5007.316 . iP» 0) (0.490) | 2.008 | 1.514 
1.012 1.514 
2.00 
5001.469| 3p 8p— 3d 3 Ff 0) 0.976 | 1.166) 1.071 
5005.140 —D Fy 0) 1.125 | 1.333} 1.250 
5025.665 72) i J 0.747 1.333 | 1.086 
5012.026| 3s’°5P; —3p'5P (0) 1.664 | 1.664! 1.664 
5002.692! 3s %Py 3p 3S, 0) 2.022 - | 2.022 
5010.620 P 1.460 | 1.460 | 2.016 
2.016 
5045.098 sP, 0) (0.515) | 1.505 | 2.015 
1.000 1.495 
5175.89 | 3p'°D 3d’ °F, 0) 1.087 
5179.50 sD, I 0) 1.073 
5452.12 | 3p%Py) —3d 8P 0) 1.478 1.478 
5454.26 3P, — 3P, (0) 1.528 | 1.528 - 
5462.62 3P, — iP 0) 1.508 | 1.533 | 1.483 
5478.13 3P, — 8P, 0) 1.486 | 1.528 | 1.500 
5480.10 3P, 3p 0) 1.502 | 1.495 | 1.481 
5495.70 3P, — P» 0) 1.496 | 1.492} 1.500 
5535.39 | 3s'®P3.1—3p' Dg. 0) 1.079 
5666.64 | 3s *P; —3p 3D 0) (0.282) | 1.448 | 1.167 
O.888 1.167 
5676.02 sPy — D 0) 0.487 487 
5679.56 2? —- pM 0) 1.104 | 1.532 | 1.319 
5686.21 P, - 3D, 10.501 (0.935)) 1.442| .501 
1.444 
5747.29 | 3s 'P, 3p 8D» 0) 1.216 | 1.056} 1.163 
5927.82 | 3p3%P,y —3d 3D 0) 0.493 dle 493 
5931.79 3P; — 8D, 0)0.910 | 1.53 | 1.115 
5941.67 3, — %) 0) 1.143 | 1.509 | 1.325 
a 3d'P, —4p 'S 
6167.82 |\3g3F, —4p 3D, | (0) 1.146 | 1.250} 1.319 
6170.16 3F, — sD, 0) 0.777. | 0.682 | 0.50 
6173.40 3F; — 3D, (0) 0.913 ? | 1.086 | 1.259 
6284.30 | 3p'D, —3d 'P 0) 0.987 | 1.000 1.026 
6379.63 | 3s 3P, —3p 'P 0.450 1.455 | 1.005 
6482.07 | 3s '!P, —3p 'P, 0) 1.026 | 1.047 | 1.005 
possible), the g values from the unresolved 


H N. MAXWELI 
ABLE II. Configurations 1s?2s?2 p3s 
TERM g(obs g(LS g (H 
35 8P, 1.502 1.500 1.500 
3s P, 1.051 1.000 1.055 
3p 1.455 1.500 1.445 
2.506 2.500 2.500 ze 
35 38P, 
PaBLe II]. Configu on 1s22s?2 p3 
TERM rs 
+p 1So 
3P, 
3p 'P 1.005 1.000 
Ss; 2.015 2.000 
P 1.530 1.500 
3p, 0.494 0.500 
5.044 5.000 »& 
3p'D 1.002 1.000 
P, 1.497 1.500 
D 1.166 1.167 
3.665 3.667 > 
3p 3p 1.330 1.333 
PABLE IV. Configuration 1s?2s?2 p3 
TERM g(obs 2(LS 
3d 8P, 
3d 1P 1.026 1.000 
P, 1.487 1.500 
dD, 0.494 0.500 
3.007 3.000 De 
3d ‘Deo 0.986 1.000 
3P. 1.504 1.500 
sD, 1.114 1.167 
3F, 7 0.667 
4.333 De 
3d '} 1.000 
3p) 1.329 1.333 
3F; 1.079 1.083 
] $17 >> 
3d *F, 1.250 1.250 


In our previous work, we have stated that the 





patterns were determined. ' 

The g values thus obtained are summarized in 
Tables II, III, IV. The shown are 
weighted averages. 


® See e.g., Shenstone and Blair, Phil. Mag. 8, 765 (1929), 


results 


g values obtained have an accuracy of about one- 
half percent. In the present work, due to the lack 
of sharpness of the lines, we must increase the 
limit of error, and suggest that this accuracy is 
not better than one percent. An inspection of the 





OO ee er | 








+ OPP TIE os Wee pee 





REFRACTIVE IN 


tables shows that wherever it is possible to give 
g sums, that they show an accuracy of consider- 
ably better than 1 percent, but we believe this 
result to be entirely fortuitous. 

The configuration 1s°2s°2p3s is compared with 
Houston's’ calculations for this configuration and 
the two levels with j= 1 show considerable devia- 
tion from LS coupling and gratifying agreement 
between theory and observation. 

Of particular interest in connection with the 
configuration 1s*2s?2p3p, is the line 46379, which 
has a pattern belonging definitely to a transition 
in which Aj=0, thus showing that Fowler and 
Freeman’s and Pretty’s designation of the level 
3p'P, (called by them 3p 'Dz) is incorrect, as 
pointed out by Edlén, on the basis of selection 
rules. 

The configurations 1s?2s°?2p3p and 1s°2s°2p3d 
show a large perturbation of the middle term of 
the middle triplet. 3p *P,; has a g value of 1.530 
and 3d *Dz has a g value of 1.114, the LS g values 
being 1.500 and 1.167, respectively. Both of these 
results are too far different from the LS values to 
be accounted for on the basis of experimental 


7 Houston, Phys. Rev. 33, 297 (1929). 
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error. No other levels of the configurations are 
close enough to produce perturbations of so large 
a magnitude. Nevertheless there are perturba- 
tions in both cases, as is evidenced by the fact 
that the 3p *P separations are in the ratio 1.65 : 1 
instead of 2 : 1; and the 3d *D separations are in 
the ratio 2.52 : 2 instead of 3 : 2 for LS coupling. 
In addition to this, it has been pointed out by 
Condon and Shortley® that the intermultiplet 
separations are also considerably different from 
the LS values for the 2p and 3p configuration for 
the whole isoelectronic sequence C I, N II, OITI, 
being in the ratio of 2.27 
P—D/D-—S. They attribute this to interaction 
between the 2p and 3 configurations. There 
could also be ccafiguration interaction between 


: 2 instead of 3 : 2 for 


either of these configurations and the configura- 
tion 1s°2p*; and between the 3d configuration and 
1s*2s2p*. Perturbations of this type are electro- 
static, however, and their effects on g values 
would come as a result of their spoiling the elec- 
trostatic parameters which would in turn perturb 
the LS ratios and thus secondarily affect the g 
values. 

§Condon and Shortley, Theory of Atomic Spectra 
(Cambridge University Press, 1935). 
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The Refractive Index of Water for Electromagnetic Waves Eight to Twenty-Four 
Centimeters in Length 


THomas T. GOLpsMITH, JR. 
Cornell University, Ithaca, New York 


(Received October 14, 1936) 


Wave-lengths were measured in air and water for continuous waves produced by magnetron 
and positive grid oscillators. Although high in absolute value, the index of refraction thus 
determined shows a decrease with increasing frequency, indicating a start of the drop toward the 


infrared value. 
INTRODUCTION 


ARLY attempts':? with electromagnetic 
waves from spark sources have been made 
to detect experimentally the beginning of the 
dispersion region where the index of refraction 
of water drops from the long wave value of 
u=/e=9 toward the infrared value of ».=1.33. 


This paper reports an investigation in the 


! Nichols and Tear, Nat. Acad. of Sci. Proc. 9, 211 (1923). 
* Glagolewa-Arkadiewa, Nature 113, 640 (1924), 


interval of wave-lengths from eight to twenty- 
four centimeters, making use of the improved 
stability obtainable with the magnetron and 
positive grid oscillators as sources of short 
electromagnetic waves. 

Debye*® and Bernal and Fowler‘ give some 
theoretical discussion indicating probable dis- 
persion in this region. 


3 Debye, Polar Molecules (Chemical Catalog Company, 
1929), p. 85. 
* Bernal and Fowler, J. Chem, Phys. 1, 515 (1933), 
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Fic. 1. Wave-length measuring system. 


EXPERIMENTAL METHOD 


The wave-length measuring system, Fig. 1, isa 
modification of that described by Miesowicz.® 
Polarized radiation from the transmitting an- 
tenna 7, reinforced by the reflector R, passes 
between the screens S and up through the tank 
of water, reaching the crystal detectors C; and C2 
over the indicated paths. It can be seen that each 
crystal receives two beams of energy, one of 
which has been retarded over the other by pas- 
sage through two thicknesses of the water. By 
interference, the energy received at either crystal 
goes through maxima and minima as the water 
thickness is varied. The spacing between suc- 
cessive maxima or minima of this interference 
curve, plotting current from the crvstal detector 
against the water thickness, is one-half the wave- 
length of the waves in water. Similarly, the 
movable reflector shown above the crystal C,; 
produces interference from which the wave- 
length in air can be determined. Fig. 2 shows 
several typical curves for both air and water. 

The transmitters used are similar in general to 
those described adequately in the literature®—"! 


5 Miesowicz, Bull. Acad. Polonaise Sci. et Lettres 3—4A, 
95 (1934). 

6 Potapenko, Phys. Rev. 39, 625; 40, 988; 41, 216 (1932). 

7? Kilgore, Proc. Inst. Rad. Eng. 20, 1741 (1932). 

8 Megaw, Inst. Elect. Eng., London 72, 313 (1933). 

® Cleeton and Williams, Phys. Rev. 45, 234 (1934). 

10 Wolff, Linder and Baden, Proc. Inst. Rad. Eng. 23, 11 
(1935). 

™ McPherson and Ullrich, Inst. Elect. Eng., London 78, 
629 (1936). 
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except that the magnetron tubes were designed 
and built to be useful over a range of an octave 
for each tube. The positive grid tubes were of 
conventional design. 

A more elaborate interferometer system, de- 
signed to measure the absorption coefficient as 
well as the index of refraction, was built and 
tried but proved unsatisfactory because of the 
intense diffracted energy which came around the 
mirrors employed in the interferometer. Its 
mirrors were 12 by 16 inches in size. Satisfactory 
results might be obtained with a much larger 
system, though a greater beam intensity would 
be required. 

The curves with C2 proved more satisfactory 
than those with C,, because the interference with 
Cs is superposed upon a straight zero line, while 
the interference with C,; is superposed on the 
absorption curve for the water, making the peaks 
more difficult to locate. This can be seen by 
inspection of the water curves of Fig. 2. All 
curves were interpreted by taking as many 
independent peak-to-peak and trough-to-trough 
intervals as possible. Several curves were taken 


at each wave-length. 


Tc, 40 Air 
Galv. 2y RANKA, ; _ a Curve 
Div fe) ¥ Vv ¥ 7 st ! 

60 \ Air 

Te, 40 \ o, Curve 
ieee A ie 
Div Neo! aap oe : - 

D 4 & 2 16 20 24 26 32 36 4044 48 
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%® 1: 2 3 4 § 6 
Water Depth (cm) 
Fic. 2. Typical air and water interference curves. 


Curve 1. Air. A\o=8.53 cm. Curve 2. Air. Ao =23.82 cm. 
Curve 3. Water. 3a, interference curve for crystal Cs, 
Fig. 1, 3b, reaction on plate current of positive grid tube 
transmitter, 3c, interference curve for crystal C;, Fig. 1. 
From curves 3a, 3b and 3c, mean \= 1.89 cm, 








, 
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RESULTS 


Table I gives a summary of results for each 
wave-length and Fig. 3 shows the index of 
refraction for water, plotted against frequency. 
The two positive grid tubes had almost the same 
wave-length. They offered the added advantage 
of a third indication of interference, exhibited as 
reaction in the plate circuit of the oscillator itself 
as the water level was varied. This reaction is 
objectionable, however, and may introduce an 
error which could account for the fact that 
neither of the points determined with these tubes 
fall on the curve described by the magnetron 
points in Fig. 3. 

The water temperature, recorded in Table I, 
was not the same for all the runs, but correction 
for this difference” would not shift the points 
appreciably in Fig. 3. The variation of index of 
refraction with temperature for long radio waves 
amounts to 0.23 percent decrease in index per 
degree increase in temperature. The extreme 
points on this curve were taken at the same 
temperature. 


CONCLUSIONS 


Fig. 3 shows a decrease of the index of re- 
fraction of water with increasing frequency, in 
agreement with the suggestions of Debye* and 
Bernal and Fowler.‘ However, the absolute value 
of the index at the low frequency end of the 
curve is greater than that given by p= \/e= 1/81, 
where ¢ is the dielectric constant of water for 
frequencies of the order of a megacycle. An 
examination of curves 1 and 2 of Fig. 2 at once 
suggests that the air wave-lengths for the longer 


TABLE I. Summary for H,0. 


MEAN | MEAN 


RUN | Aofcm A(cm) | Temp _ 0 yp =X 100 
No. | IN AIR | IN H2O ( TUBE USED re do 

I 23.82 | 2.45 | 28 | Magnetron #1 | 9.73) 12.6105 
II | 17.76] 1.88 | 23 | Magnetron #1 | 9.44) 16.9 

III | 17.42] 1.89 23 | Positive grid #1 | 9.23) 17.2 

IV | 17.02} 1.85 | 24 | Positive grid #2 | 9.22) 17.6 

V | 12.65} 1.40 | 21 | Magnetron #1 | 9.04) 23.7 

VI | 11.48] 1.30 | 25 | Magnetron #2 | 8.83) 26.2 

VII} 8.53} 1.02 | 28 | Magnetron #2 | 8.36} 35.1 





2 Lattey, Gatty and Davies, Phil. Mag. 12, 1019 (1931), 
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Fic. 3. Refractive index uw for H.O. 


waves are not so accurate as for the shortet 
waves, due to the experimentally limited motion 
of the wave meter, allowing only a few peaks to 
be obtained for large } while many could be 
obtained for small \. The type of anomalous 
dispersion discussed by Debye?’ for this spectral 
region would lead to a monotonically decreasing 
value of » with increasing frequency as shown in 
the /nternational Critical Tables.‘* However, the 
experimental curve of Fig. 3 indicates a rise in the 
refractive index in the region between low fre- 
quencies and those employed here, and thus 
indicates a departure from the theory. 

It is anticipated that a colleague in this work 
will carry out plans to investigate this method 
further. A beat frequency method of wave-length 
comparison should prove valuable. 

The method employed here should be appli- 
cable to the study of the index of refraction of 
many fluids having low conductivity. 
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13 International Critical Tables 6, 77 (1929). 
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Approximately Relativistic Equations for Nuclear Particles 
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Cosmic-ray showers indicate that at high energies inter- 
action between nuclear particles is concerned with the 
creation and destruction of matter. It may, therefore, be 
expected that a complete relativistic theory of nuclear 
forces will involve explicit reference to the phenomena 
described at present as the electron-neutrino field. Theories 
of this kind are still too incomplete and self-contradictory 
to be reliable in practical work. It is, nevertheless, possible 
to set up equations which are relativistically invariant for 
transformations involving low velocities. Such equations 
form the subject of the present report. They are. restricted 
to energies of relative motion that are small in comparison 
with the rest mass. By means of them it should be possible 
to discuss relativistic effects for ordinary nuclear energy 
levels. Possible forms of classical equations contain an 
interaction energy between two particles in the form given 
by Eq. (13.2). Here a, 6 are arbitrary real constants. The 
vector from particle 1 to particle 2 is r, the velocities of the 
particles are Vi, V2, the velocity of light is c. If a=b=1, one 
obtains a generalization of Darwin's equation, which 
describes the motion of electrically charged particles. For 
a=-—1, b=1 particle 1 acts on particle 2 approximately as 
though it produced a scalar potential field responsible for 
the acceleration of particle 2. The requirement of invariance 


T has been noticed by Inglis! that the apparent 

inversion of doublets in nuclei can be ex- 
plained as a result of taking into account the 
effect introduced by Thomas in the discussion of 
the fine structure in atomic spectra. A nuclear 
particle, such as a neutron, is subjected to 
forces the actual nature of which is not known. 
There is no reason to consider these forces as 
being representable by an electromagnetic field 
and it is more likely that they owe their origin 
to the electron-neutrino phenomena. The direct 
influence of the electric field on the spin which 
is important for the fine structure of atomic 
spectra is thus presumably absent in nuclei. 
Very schematically one may consider a single 
nuclear particle as moving in the combined field 
of the other particles. In this rough approxima- 
tion one may attempt to regard a single nuclear 
particle as exposed to the field of a scalar 
potential due to the remainder of the nucleus. 
This potential will be provisionally supposed to 
be a function of space-time and to be independent 


for wave equations of particles with spin (Pauli types) 
makes it necessary to have spin-orbit coupling which 
should give rise to the fine structure of nuclear levels. 
For ordinary interactions the spin-orbit energy may have 
the form given by Eq. (15.4), where 0 is an arbitrary real 
constant, p is the momentum and @ is Pauli’s spin matrix. 
For 5;;= —1 one obtains the type of coupling taking place 
between extranuclear electrons. If b;;=1 each particle 
interacts only with its own orbit as though it were moving 
in a scalar field. It is the latter hypothesis that is simplest 
and corresponds to a= —1, b=1 of the classical equation. 
Extensions of the above classifications have been made to 
the Majorana [Eqs. (15.7), (15.8) ] and the Heisenberg 
[Eq. (15.9) ] exchange interactions. The simplest type in 
the Majorana case appears to be in satisfactory agreement 
with experiment for Li’ and agrees in order of magnitude 
with other cases. Extensions to Dirac’s types of equations 
have been made. They lead one to expect coupling between 
spins of nuclear particles in apparent qualitative but not 
quantitative agreement with experiment for the deuteron. 
This agreement is not sufficiently good to establish a 
form of interaction energy but indicates a possibility of 
doing so in the future. 


of the time in the coordinate system in which 
the nucleus is at rest. For such a model it 
follows that the Thomas term is the only one 
that exists for a neutron and that for a proton 
it should be combined with the effect of the 
electric field of the nucleus that is expected from 
the Dirac equation for a particle in a central 
field. 

A proof of this has been already given by 
Furry,” who saw the possibility of inserting the 
scalar as an addition to the mass term mc? in 
Dirac’s equation. The same result can be ob- 
tained by considering the transformation proper- 
ties of wave functions. The latter method is 
desirable in order to bring out the connection 
with the treatment of many particles as well as 
in order to establish the interactions inde- 
pendently of special Dirac equations. It will be 
given in Section 1. 

The neutrons and protons in a nucleus have 
approximately equal masses. In this respect the 
nucleus differs qualitatively from an atom which 


248 








=S |= afew se rn kth 


ont OP Oe 





oe 





RELATIVISTIC EQUATIONS 249 


is composed of a heavy nucleus and several 
light electrons. Even in an atom the spin-orbit 
interactions are describable by a central field 
only somewhat accidentally. The correct de- 
scription must be made by considering inter- 
actions between all pairs of particles.* Calcula- 
tion shows® that the naive central field picture 
is approximately correct mainly because electrons 
in inner shells are so firmly bound that their 
action on a valence electron is comparable with 
that of a static distribution of charge. No such 
simplicity is apparent in a nucleus. It is therefore 
desirable to understand what kind of spin-orbit 
and orbit-orbit interactions may be postulated. 
In Dirac’s theory spin effects may usually be 
considered as of the order v*/c? in comparison 
with the kinetic energy, where v is the velocity 
of the particle and c¢ is the velocity of light. 
It is, therefore, reasonable to require that the 
theory of nuclear particles be relativistically 
invariant to this order. It is presumably im- 
possible to have a completely relativistically 
invariant theory for the interaction between 
several particles with a finite range of force and 
without explicit reference to a field. It is super- 
fluous, therefore, and probably fruitless in the 
discussion of the present problem to require 
complete invariance. In Section 2 possible 
classical interactions between particles will be 
discussed and enumerated. It will be found that 
those among them that have special physical 
interest form a two-parameter set. Two types 
deserve special mention. One of them is the 
most immediate extension of Darwin’s equation 
for two charged particles interacting according 
to the laws of classical electrodynamics. The 
other corresponds to each particle setting up a 
scalar field in a reference system in which it is 
at rest. It is this type that is related most closely 
to the interaction discussed in Section 1. In 
Section 3 possible types of wave equations will 
be discussed. It will be seen that: (a) The 
presence of spin-orbit interactions is a general 
consequence of the transformation properties of 
wave functions; (b) Their exact form is not 
uniquely determined by the transformation 
properties, but all possible forms of physical 
interest may be considered as linear combina- 
tions of electromagnetic interactions with pure 
Thomas terms; (c) A possible extension of the 


scalar one-particle picture gives only the Thomas 
terms. The effect of Majorana and Heisenberg 
exchange will be also considered in this section. 


1. ONE PARTICLE 


It will be supposed that the particle in its 
free state is describable by Dirac’s equation. 
The fact that the magnetic moment of nuclear 
particles is not that given by Dirac’s equation 
in its original form does not contradict this 
supposition, because it is possible to introduce in 
Dirac’s equation direct interactions between the 
magnetic field and the magnetic moment in 
the manner proposed by Pauli.‘ The transforma- 
tion properties of Dirac’s equation are the same 
whether a true magnetic moment of the Pauli 
type is introduced or not. It is on the other 
hand questionable whether Dirac’s equation or 
any of its modifications gives an exact account 
of the properties of nuclear particles. In fact 
high energy collisions involving relative velocities 
of nuclear particles comparable with c¢ indicate 
from experiment,® in the light of Heisenberg’s 
theory® of showers, the importance of processes 
involving the creation and destruction of par- 
ticles. For this reason it is preferable to base 
theories of nuclear spin-orbit interactions on 
corrections of the order v?/c? rather than on an 
exact form of a single particle equation. 

The original discussion of Thomas does not 
apply directly to particles whose spin is de- 
scribable by means of a Pauli or a Dirac equa- 
tion. For a Dirac particle the spin is a secondary 
characteristic. As has been pcinted out by 
Schroedinger,’ a Dirac particle is unable to move 
uniformly in a straight line when left to its own 
devices. This property of the particle gives rise 
to the spin. It is, therefore, necessary to consider 
the Thomas term directly from the point of 
view of the wave equation rather than in terms 
of a classical analogy. 

According to Dirac* the Lorentz transfor- 
mation 


x’=x ch 6—ct sh 8, ct’ =ct ch 6—x sh 8, 


po’ = po ch 6— pz sh 8, 
(1) 
sh = (v/c)(1—v*/c*)-1; 


p:' =p; ch O0— po sh 6, 
ch 6=(1—v?/c?)-, 
pPo=E/ce (ch 6=cosh 6, sh 6=sinh @) 
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induces the transformation 


y’ =e 2y, yp’ * = y*er10/2 (2) 


on the wave function. This means that the trans- 
formed wave function y’ satisfies the Dirac 
equation in the transformed system K’ and that 
all physical quantities calculated by means of 
y’ in K’ are related by relativistic formulas to 
corresponding quantities calculated by means of 
y in the original system K. 
The Dirac equations can be written as 
(po+ Mc)®+ (op) ¥ =0, 
(3) 
(po— Mc) ¥+(ep)>=0, 


where the o are Pauli spin matrices, V is the 
two-component Pauli wave function, and ® is 
also a two-component wave function. Both @ and 
WV are column matrices. The Dirac function y is 
a column matrix having the two elements of ® in 
first and second place and the two elements of V 
in third and fourth place. From (3) one obtains 
in sufficient approximation 


= — (ep) ¥/2Mc. (3’) 


This approximation is valid for interactions of 
both the electromagnetic and the scalar types. 
The occurrence of terms in the vector potential 
in (3) leads only to negligible corrections in (3’) 
because these terms can be taken to be of the 
order e*v/c’r, where r is the distance between 
nuclear particles and v is the velocity. This 
quantity is of the order of mv for r~e*?/mce’, 
where m is the electronic mass, while the opera- 
tions involved in —ihd/dx give rise to Mz, 
where M is the nuclear particle mass. Similarly 
the addition of a scalar interaction energy to 
Mc? or of an electrostatic potential to po gives 
rise to insignificant corrections for forces having 
the approximate range e?/mc’. 
Using (3’) in (2) 

W’= {ch (6/2) 

—(2Mc)“(p+i[pXe]). sh (0/2)}W. (4) 
Let now the wave equation be of the form 

jh 0 ‘ Pp p* r 
—+ Me?+—_—-——_—-J+Q;¥=0. (5) 
lia 2M 8Mic2 


Substituting V in terms of ¥’ and expressing all 
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other quantities in K in terms of quantities in K’, 
the wave equation should return to its original 
form. It will be supposed that J is a scalar. 
It will thus transform into itself. The quantity Q 
is introduced into Eq. (5) in order to correct the 
remaining terms for lack of invariance. Instead 
of considering all of Eq. (4), it will be first 


approximated by 


VW=SW’; S=1+iv[pXe]/4Me. (4’ 
This approximation neglects some effects of 
order v*/c?. However, an inspection of Eq. (4 
shows that effects of order v*?/c? are additive and 
no essential error is committed.® Substituting 
Eq. (4’) into Eq. (5) and multiplying by S“ 
from the left one obtains on account of the 
presence of J the term 


J—S“IS=(h/4Mec)[vJ Xe ]v. (6) 


This term must be canceled by the introduction 
of a suitable term into Q. It will be supposed 
that Q contains only terms of order v*/c?. There- 
fore one may take S"'QS=Q. Thus in the 
approximation of Eq. (4’) there is a part Qo of 
Q such that Qo’—Q> must be equal to the right 
side of Eq. (6). This can be accomplished by 


Qo= —(h/4M?C) [9 J Xe Jp. (6’) 


This is the Thomas term which must, therefore, 
be included in order to enforce invariance. 
In addition it is possible to postulate the presence 
of other terms in @ entering Q. Such terms must 
form an invariant to the order v?/c?. It is im- 
possible to form such an invariant by having 
VJ, «, and p entering the expression once. 
There exists, in consequence, no term which can 
be added to Q which has an essential relation 
to Qo. 

Transforming Eq. (5) to variables in K’ one 
finds that 
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The terms containing v have to be compensated 
by Q—Q’. A possible Q, determined in a manner 
analogous to that in Eq. (6’), gives an invariant 
form by substitution into Eq. (5) 


Pp p' 


, p’J+pJp 
2M 8 Mic? 


4.\f?c? 


| -E+Me+ 


h 
= ((pxvJ]—[wsxp]o!¥=0. (8) 
81l*c? 

The term in @ is here symmetrized so as to be 
Hermitean. The terms in p*/ and p/p form a non- 
Hermitean operator. This is due to the fact that 
the total probability of finding the particle is not 
given by the volume integral of | ¥}* but by 


fvrvar= [verte 4\P7c?)Vadr (8’) 


It is possible, however, to transform the wave 
function WV so as to have a Hermitean operator. 
Thus if 

PWY=(14+ p?/8 Mc?) v (9) 


Eq. (8) becomes 


Pp p* 
| _ £4 Me+——— ay 
| 2M 8M 
h : | 
- ([pxvJ]—[wJXp])o}¥=0, (9) 
8 \M2c? 


pS +2pJp+Jp 
Fs. a 


8 M22 


which is Hermitean. The modified function ¥™ 
satisfies the simple normalization condition 


[v.orveoar= fusrvdr=1, (9”’) 


which follows from Eqs. (8’) and (9). An elegant 
treatment of the single particle problem due to 
Furry’ can be carried out introducing the scalar 
—J into the Dirac equation as an addition to 
the rest mass energy Mc*. The equation in the 
Dirac form is 


{pot (ap) +8(Mc—J/c)}y=0. (10) 


Expressing the components ® in terms of V by 
means of Eq. (3’) and substituting into the 
second Eq. (3) modified by Mc—>Mc—J/c, it is 
found that W satisfies Eq. (8). This is a natural 
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result inasmuch as Furry’s Eq. (10) is strictly 
invariant. It should be noted that the transition 
from Eq. (7) to Eq. (8) is not unique. Thus it is 
possible to use instead of p°J+pJp in this 
equation the quantity (3p°J/+Jp*)/2 without 
changing the transformation properties, and 
more generally an arbitrary constant times 
p’J+Jp?>—2pJp can be added to p°>J+pJp. 
There is thus a greater freedom in the forms of 
two-component equations invariant to order 
v/c? than is apparent from Eq. (10). The equa- 
tions obtained in this manner are identical, 
however, from the point of view of classical 
analogy, and in the absence of another clue it is 
impossible to give reasons for preferring one to 
another. The particular form (8) may perhaps be 
advocated on the grounds that it is equivalent 
to Eq. (10). 

It is possible to add to the operators in 
Eqs. (8), (9’) any quantity invariant to order 
v*/c?. Thus J?/Mc* multiplied by a reasonably 
small number may be added. The addition of a 
term in this quantity is, however, only of 
trivial interest in the present case because such 
a term may be considered as being incorporated 
in J. The addition of commutator-like terms 
such as (¢?J+Jp*?—2pJ/p)/Mc* amounts to add- 
ing terms in derivatives of J and also offers no 
essentially new possibility. Terms in p*/ and 
p/p can occur only as in Eq. (8) to within the 
freedom of adding a commutator. Using ~* and 
p* it is impossible to have an invariant except that 
already occurring in the equation. Within the 
above discussed set of the simplest possibilities, 
Eqs. (8) and (9’) are thus essentially unique. 

The parts of the operators occurring in Eqs. 
(8), (9) which do not contain @ can be in- 
terpreted as corresponding to the wave equation 
of a particle without spin, invariant to order 
v’/c*. The verification is analogous to the step 
from Eq. (7) to Eq. (8). A slight difference from 
Eq. (8) is found in the order of factors. Allowing 
the wave function to transform itself by 


Ww" =(1407/4e—vp,/2Me)v — (11) 
it is found that 
P’ p* 


= an 
| 2M 8Mice 


p/P | 
a v’=0 
2M22\ (11’) 


is an invariant form. 
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The transformation of Eq. (11) gives 


Cworwoy’ _ ( ee :) yor? 
2c? 


hv aw? av 
Scena yor — we... 5 

2UMec Ox Ox 
One may, therefore, interpret W*W as 
particle density and use the normalizing con- 


fvowear =1. 


In obtaining Eq. (8) it was supposed that the 
contributions due to time derivatives of J arising 
in the transformations are small. This require- 
ment can be satisfied by supposing that there 
exists a frame of reference in which J is a function 
only of space and does not depend on time. 
This frame may be called the rest frame. The 
other frames considered above are those obtained 
by the application of a small translational 
velocity to the rest frame. For such transforma- 
tions VJ is an invariant since it changes only by 
terms Of order (v?/c?)¥J on transformation. If J 
is only a function of the distance r from a fixed 
point there is no necessity of symmetrizing the 


the 


dition 


(11) 


last term containing o because then 
[pxVJ]=-—[VJxp]. 


If J, instead of transforming as a scalar, 
transforms as the time component of a four- 
vector, one obtains essentially the same results 
as are derivable from the ordinary Dirac 


equation. 


2. CLASSICAL EQUATIONS FOR SEVERAL 
PARTICLES 


It will be supposed that the system can be 
treated by means of a variational equation, 


6 f zat =(). 


The trajectory of the dynamical system in the 
reference system K corresponds to segments of 
world lines in an interval t;<t<t:. The La- 
grangian L will be supposed to be given approxi- 
mately by 


(12) 


Lo= >. M yv;? 2+>- J 


a 


(rij), (12’) 


: ae | 


where M; are the masses, v; the velocities of the 
particles, and r;; are the distances between them. 
In addition to Ly the function LZ will be supposed 
to contain terms of the order (v*/c*?)Lo. These 
terms will now be determined by the requirement 
that Eq. (12) should determine the same world 
lines independently of the coordinate system 
used. For a system of noninteracting particles 
(J;;=0) there is an 
sponding to describing each particle separately. 


obvious solution corre- 


In this case one may use 


Lo= > Mic{1 —(1—0,2/c)*"]. 


= 
t 


Expanding each term of this sum one obtains 


(J=0) 


Lo= >. M,(0;?2/2+0;4/8c?) 
i 

to the required order. For this Lo the integral in 
Eq. (12) is not invariant because: (a) each term 
of the sum contains M,c? fdt which is not in- 
variant. However, the variation of this quantity 
is zero and its inclusion is harmless. (b) In the 
transformation K—K’ it is impossible to deal 
with identical segments of world lines in both 
systems. If the end points of the world line for 
1 are kept the same, those of 7 have to be changed 
in K’ from what they were in K. This circum- 
stance also does not matter because the varia- 
tions of the coordinates of 7 and j are supposed 
to be independent of each other. Therefore the 
variations over that portion of the world line 
of 7 which is common to K and K’ for fixed end 
points of the world line of 7 are automatically 
the same as long as each term of Ld is invariant. 
It is thus seen that 
sufficient for the invariance of the world lines. 
It is obviously not a necessary condition as is 
seen in the case of the addition of terms in /dt. 
The terms J;; substituted into Eq. (12) give 
rise to noninvariant quantities. These will now 
be calculated. It will be supposed that in calcu- 
lating the correction terms it is legitimate to 
neglect the curvature of the world lines. This 
assumption does not follow from considerations 


the invariance of Ldt is 


of invariance. It simplifies the discussion by 
preventing the explicit occurrence of accelera- 
tions in the Lagrangian. For electromagnetic 
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interactions it leads to neglecting radiation reac- 
tions. For nuclear forces it amounts to neg- 
lecting x¢#/c?~J/Mce~20m/M~10~ in com- 
parison with unity. Here m is the electronic 
mass. In terms of this estimate the approximation 
is a good one since the effects neglected are of 
the order of 10-* times the correction terms in 
v- c?. 

For the transformation given by Eq. (1) the 
distances between particles 1,2 in K’ and K are 
related by 


(vr)? (vr) M,(rv2)+1/.(rv;) 
, =o — - 
2c*r = c*r M,+Mo 
and hence 
J(r’)=J(r) (13) 
(vr)? (vr) W,(rve) + Moe(rv,)]dJ 
+|- + — | 
2c? c? M,+ Me rdr 


In order that Ldt be invariant it is necessary that 
L should transform as ds/dt where ds is the 
differential of the four dimensional distance. 
The transformation formulas for ds/dt are 
different for the world lines of 1 and 2. The 
world line of the center of gravity of 1 and 2 
gives a still different transformation of ds/d?. 
At this stage, therefore, the definition of the 
world line is still arbitrary, since one has the 
above three choices, each of which has some 
physical plausibility. This arbitrariness of the 
choice of the world line would destroy the 
possibility of having a sensible description of 
the system by means of Eq. (12) if it were not 
for the fact that the results of Eq. (12) are 
independent of the choice in the following sense. 
Let ,, Mz in Eq. (13) be given arbitrary values 
and let the same values be used for the definition 
of the path of the center of gravity. Let the 
transformation of Z be similar to that of ds/dt 
along this path of the center of gravity. The 
results of Eq. (12) will be seen to be independent 
of the values of 4, and M2. This means that it is 
possible to use either the path of 1 or that of 2 
and that the results of the two procedures are 
consistent. 

For the path of the center of gravity for the 
transformation of Eq. (1) 


ds | v M,(vv;)+ Me me 

== 1 — - _, 
dt’ 2c? (17, +Me)c? dt 
It is now required that J—Q should transform as 
ds/dt. Since Q contains by definition only terms 


of order v*/c?, the transformation formula for Q is 
determined by means of Eqs. (13), (13’). It is 


v= (vv,) + (vv2) 
0-0’ =| ——+ J 
. 2c 2c? 
(vr)? (vr) yd J 
~ ((Vi+Ve) -r) 
2c? 2c? rdr 
M,—Ms. d 7 
~ —[(rv)J], where r=r,— fm. (13"’) 


2c?(M,+ M2) dt 


The last term in Eq. (13’’) may be disregarded 
because in Eq. (12) it gives rise only to a function 
of the end points of the path. The remainder of 
Q—(Q’ is independent of 1, and M, and is equal 
to Q—Q’ for M,=Mz2. It is symmetric in the 
coordinates and velocities of 1 and 2. It is 
natural to look for possible Q which are also 
symmetric in 1 and 2. Unsymmetric Q are either 
reducible to symmetric Q by the addition of a 
complete differential to Ldt or else they imply a 
different action of particle 1 on particle 2 from 
that of 2 on 1. Such a law of interaction is 
improbable and will not be considered here. 
Only symmetric possibilities for Q will be dis- 
cussed. 

A possible J/—Q with Q transforming itself as 
in Eq. (13”’) is 


J-—Q=J—([(vive)J 
—(vir)(ver)dJ/rdr|/4c?, (13.1) 


as is verified by applying the transformation of 
Eq. (1) and observing that for the transformation 
of Q it is sufficient to use Newtonian rather than 
relativistic addition of velocities and to consider 
r as invariant in this connection. All other Q can 
be obtained from Eq. (13.1) by the addition of 
invariants. In this way one obtains among other 
possibilities the following expression for J—Q: 


J—Q=J-— {2a(viv2) + (1 —a@) (012 +22") } J/4c? 
+ {2b(vir)(ver)+(1—)[ (vir)? 
+(ver)? |}dJ/4c*%rdr, (13.2) 
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where a and 2 are arbitrary real constants. For 
a=b=1, Eq. (13.2) agrees with Eq. (13.1). 
In this case the Lagrangian Lo +J—Q is a 
generalization of that obtained by Darwin for 
electromagnetic interactions, as is seen by sub- 
stituting e*®/r for J. It will be recalled that 
Darwin’s approach was that of retarded po- 
tentials and that the Lagrangian was found by: 
(a) working out the equations of motion for 
particle 1 using position and velocity of 2 as 
parameters, (b) finding the Lagrangian for 1, 
and (c) showing that there exists another sym- 
metrical Lagrangian differing from the first one 
by a time derivative. It is desirable to explain 
why Darwin’s equation is of type (13.1) rather 
than of some other type contained in Eq. (13.2). 
Obviously the distinguishing feature of the form 
(13.1) is that of antisymmetry with respect to 
reversal of sign of one of the velocities. All other 
types in Eq. (13.2) have another type of sym- 
metry. In the electromagnetic case corrections 
for the velocity v; enter for two reasons. In the 
first place particles 1 and 2 interact magnetically. 
This interaction is antisymmetric in v, because 
it depends on (v;v2). In the second place cor- 
rections for retardation must be made. These 
corrections can depend on the direction of 
motion of 1, while one discusses the motion of 2 
only because 1 is moving away or towards 2. 
Hence again the addition to the Lagrangian 
must be antisymmetric in Vv}. 

It can be verified by explicit calculation that 
the forms given by Eq. (13.2) give an expression 
for the force which is in agreement with rela- 
tivity kinematics to the order v?/c?. According to 
relativity if P, stands for W,v,(1—v;?/c?)~'? then 
at the same point in space time the rates of 
change of P; and P;’ in K and K’ are con- 
nected by 


P, =P,’+[—v*P,’+v(vPy’) ] 2c? (14) 


dx’ dy’ dz’ 
(—- —=- -0), 
dt’ dt’ dt’ 


where Vv is the transformation velocity, provided 
terms of order higher than v?/c? are neglected and 
provided P,’=0. The forces on 1 are thus 
determinable in terms of the force in that frame 
in which 1 is instantaneously at rest. The 
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verification of Eq. (14) is straightforward but 
somewhat lengthy and is not reproduced here. 
It involves:no approximations except that of 
neglecting all terms of order higher than v?/c?. 

If one sets a=b=0 in Eq. (13.2) a simple form 
for J—Q is obtained : 
orton? 1 dJ 
J-—Q=J—-——J+—[ (vir)? + (ver)? }—. (14.1) 

4c? 4c rdr 


This is somewhat related to Eqs. (9’), (11’) for a 
single particle but does not quite correspond to 
them. Thus if particle 2 has a sufficiently large 
mass to make v2 small at all times, the above 
J—Q is not altogether equivalent to J—v*J/2: 


because 
v°J+(vir)*dJ rdr= (d dt)[ (vir) J ]— (vir)J 


and is not a complete differential. Working out 
the equations of motion for Eq. (14.1) and 
comparing them with the equations of motion 
following from Eqs. (9’), (11), no simple corre- 
spondence is found. 

A simple correspondence is found, however, 
setting a= —1, b=1. In this case 


J—Q=J—(vi—Vv2)?J/2c?— (viW2)J/2c? 
+(v)r)(ver)\dJ/2c?rdr (14.2) 
(a= —1, b=1). 


If ve=0, this agrees with Eqs (9’), (11’). 
Comparing the equations of motion following 
from Eq. (14.2) with those from Eq. (9’), it is 
found in both cases that in the reference system 
in which 2 is at rest the acceleration of 1 is the 
same function of its position with respect to 2 
and its velocity. One may, therefore, regard 
Eq. (14.2) as being a generalization of Eq. (9’ 
and of Eq. (10) if in the latter the effect of the 
spin is neglected. One may think of Eq. (14.2) 
as giving such a motion of the two particles as 
corresponds to particle 2 producing in its refer- 
ence system a scalar field which acts on particle 
1, and the other way around. 

Only two particles are explicitly considered 
above. The arguments which lead to Eq. (13”’) 
for Q—Q’ can be generalized to any number of 
particles having arbitrary interactions in pairs. 
The calculation is more lengthy but the result is 
simply that, to within a term which is a complete 
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time derivative, Q—Q’ is a sum over pairs of 
terms such as are given by Eq. (13”) for a single 
pair. An extension of Eq. (14.2) to many par- 
ticles is obtained by using a sum over pairs of 
such expressions. If particles 2, 3, --- m are 
tied by strong forces into a large mass approxi- 
mately at rest and particle 1 moves in their 
field, this generalization becomes equivalent to 
Eq. (9’) because the time average of Ve, V3, --- Va 
can be considered to be 0. 

The forms used in Eq. (13.2) can be generalized 
by the addition of other invariants. Thus one 
could add terms in (v¥,;—Ve)*7dJ/cdr or in 
((v;—Ve)r)?J/cr*. The use of such terms cannot 
be eliminated by general principles but appears 
to be an unnecessary complication. Similarly one 
can add terms in higher derivatives of J, and it 
is possible to have approximate invariants of 
order of magnitude v/c which in Lorentz trans- 
formation undergo changes only of order v*/c?. 


3. WAVE EQUATIONS FOR MANY PARTICLES 


General methods for setting up relativistic 
equations for composite systems have been given 
by Heisenberg and Pauli, Fermi, Dirac, and 
Dirac, Fock and Podolsky.'! These methods do 
not apply directly to the problem considered 
here because they are concerned with exact 
rather than approximate invariance. The ap- 
proach of Heisenberg and Pauli by means of an 
invariant Lagrange function can nevertheless be 
made use of in a modified form. Instead of 
extending the integration over space time as is 
done by them it is more natural to consider here 
variational integrals over the time and the con- 
figuration space because this is the simplest 
formulation of a nonrelativistic equation. Speak- 
ing of the motion classically, it may be expected 
that the interaction energy J should be treated 
as small. Otherwise the velocity of at least one 
particle will be large and an approximate 
relativistic treatment will have no sense. This 
assumption corresponds to the similar restriction 
of small accelerations made in the classical case 
when Hamilton’s principle was used. Forms of 


if ear=0, dr=dV\dV2---dV,dt_ (15) 


will be looked for. Here dV,, dV, dV, are 
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elements of volume for each particle. If the 
interactions are not of an exchange type and if 
the particles have no spin, it may be expected 
that the resultant wave equations should have a 
possible classical form such as one of the forms 
in Section 2. Only the occurrence of com- 
mutators in the wave equations will, therefore, 
become clearer in such a discussion. The presence 
of terms in the spin-orbit and spin-spin inter- 
actions and of the effect of exchange is, on the 
other hand, obtainable from Eq. (15) and cannot 
be inferred from a classical discussion. 

It will be supposed that if c= in £ and the 
corresponding wave equation, the system will be 
invariant under Galilean transformations. Thus 
for exchange interactions it will be simplest to 
confine oneself to particles of equal masses. In 
discussing Lorentz invariance it is easiest to 
consider 

L=F) £:+)> L,, (15.1) 

i>j 

and to require that (a) the forms £; should yield 
invariant equations for the separate particles and 
(b) all the //£i;dr should be invariant when one 
substitutes into £;; arbitrary linear combinations 
of products of single particle wave functions, 
each single particle wave function being assumed 
here to obey the wave equation of a free particle. 
Assuming Eq. (15.1), the restriction (b) will 
formally give a correct treatment only for very 
small J. Nevertheless, it may be expected to be 
also applicable, to some extent, also in cases of 
reasonably large J, because the largest part of the 
interaction energy will be Galileo invariant. Re- 
garding the interaction energy as a power series 
in 1/c, the first term of it will be thus correct by 
assumption and the remaining term will be 
established only for small J. This is less than one 
might wish to attain but it is not much worse 
than the knowledge of corresponding terms for 
interactions between charged particles. The dis- 
cussion of the classical spinless equations by 
means of their equations of motion showed that 
even though the assumption of small accelera- 
tions was made in using the variational principle, 
the transformation formulas for the force are 
formally invariant to the order v*/c? independ- 
ently of this assumption. Since the classical 
equations are better than expected it is possible 
that the wave equations behave similarly. 
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As in Eq. (4’) one may discuss the spin-orbit 
interactions separately. Corresponding to Eq. (3) 
one now has 2” equations, where 7 is the number 
of particles. Instead of Eq. (2), one has for non- 
interacting particles a transformation in which 
the exponential function of Eq. (2) is replaced by 
a product of similar exponentials. Eq. (3’) is 
replaced by m equations which determine those 
wave function components in which the Dirac 
indices for just one particle are in the small group, 
(1, 2). The extension of Eq. (2) gives similar ex- 
tensions of Eqs. (4), (4’) in which there occur 
sums over [ p; Xo; _]. In this way instead of Eq. (6) 
one has 

Qo’ —Qo= & (h4/4Mic?) 09 Ji; Xe; ]v. (15.2) 


72 
This equation is satisfied by 


Oo=— Dd (4/4 Mic?) 9 Ji; Xe. ):Saunp./M 
‘i k 

Yan=1 (15.3) 
and the a;, are so far arbitrary. Since it was de- 
cided to use for J/—Q only interactions of par- 
ticles in pairs, since the occurrence of p,/M;, in 
the interaction energy is strange, and since finally 
it is simplest to suppose that the interaction 
energy of a pair is symmetric, the following form 
is sufficiently general for most purposes: 


Oo= — DX (h/4e?) | OV Si; Xe: ]Lbipi/M2 
tm] 
+(1—0;;)p;/M:M,;)]+[¥jJi;Xe;] 
-(b:;p;/M?+(1—);,)pi/MiM;]}. (15.4) 


The 0};; may be different for every pair. For 
electromagnetic interactions? );;=—1, so that 
2p2—p:i goes with @;. For 5;;=1 an extension of 
Eq. (9’) is obtained, as is readily verified by 
letting all 17; but one become infinite. Since J;, 
is a function only of r;;, Eq. (15.4) gives a 
Hermitean operator for Qo. 

If the interaction is of the Majorana type, it is 
simplest to take the masses as equal to each other 
because otherwise the operator must be modified 
in order to preserve Galilean invariance. One 
obtains 


Qo’ —Qo= (Vv 8 Mc?)>°[ Ap X (oe; =e ) | (15.5) 


with 


EIT 


Aii= if (px — py) JiitJei(pe— pr) Pir”, (15.6) 


where P,;™“ is the Majorana operator exchanging 
the space coordinates of k and /. Possible forms 
for Qo are given by 


—1 
Ov-¥/ 7 )Cap +(t—a)p [A Xe, |} 
k= 8 2¢2 
| 1a . 
at? laure [pi X JisPir™ pi jor 
<!\4M?*c? 


1(1—a) 


4M a [ pe x JusPus ps Je, (15.7) 
M *c* 


where a is a real constant. In this case J;; rather 
than just dJx:/drx: enters into Qo. Thus for 
a=1 and a=0, respectively, 


h 
(a=1) O=> ——(¥iJx:) [Pi Xox ] 
k21(4.M?c? 
= = Jil pr x prlor | Px™, 
4M2c2 nai 
(15.8) 
h 
(a=0) Q=>r (Wier) [Pe Xor ] 
114) 2¢? 
oad Jip x pion} P ™. 
4M%c2 = 


For a=1 the value of Q) may be thought of as 


h 
(a=1) Qo=——— > [p:Xpil]o; 
4M?c? i 
h _ 
=— | Pi Xx pi je; (15 8’) 
4M?c? i 


and corresponds therefore most closely to the 
expectation on classical theory. From the point 
of view of Eq. (15.8’) the choice a=1 is the 
simplest. It should be noted that with Majorana 
interactions p;/M is only approximately the ac- 
celeration. If classical analogy were followed 
literally, Eq. (15.8’) would be more complicated 
because the operator representing the accelera- 
tion contains many terms. The expression for it 
can be substituted for p;/V/ in Eq. (15.8’). How- 
ever, there appears to be little justification for 
such a literal use of the Thomas term, particularly 
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since p;/.W does not represent the velocity. For 
exchange interactions involving Heisenberg’s 
operator, 
h 
Qo= ——— > [ap.t+(1—a)p,] 
4M?*c? kei 
PO JurXon|Per®. (15.9) 


Possible forms for Q can be discussed analo- 
gously to the classical Eqs. (13.2) introducing the 
momenta and spins into the equations. The dis- 
cussion is somewhat simpler, however, if instead 
one uses equations of Dirac’s type with the 
Hamiltonian given by 
H= — (ap) —C(a@2Po) 

—(6:+62)Mc?—J+Q (16.1) 


as well as 


H= —c(@p;) — C(@2P>2) 
~(8,+8:)Mc?—B,B2.I+0 (16.2) 


in the customary notation. 

It is also possible to combine the two types and 
to have J occur partly with 6:82 and partly by 
itself. It will first be shown that Eq. (16.1) is 
satisfactory with 


O= }(e;a2)J — 3(a;r)(a@er)(dJ/rdr). (16.3) 


Only the contributions due to J—Q to the in- 
tegral in Eq. (15) need be considered since the 
other terms give an invariant contribution. The 
Lorentz frame K is changed to K’. The integral 
in K’ will be evaluated by comparison with the 
integral in K. The world point (1) corresponding 
to (t, X1, ¥1, 21) will be kept fixed temporarily. In 
K,dV2at a world point A corresponds then to the 
same time ¢. In K’, however, the point (t, x, 
Vi, 21:)=(t', x1’, yi’, 2:') occurs in the integral 
together with points B=(t’, x’, ve’, 22’) which, 
therefore, have values of ¢ different from that of 
(1). Let the points B be associated to the points 
A by having (x, y, z)p=(x, y, 2)4. Then 


, 


dx’ = (1 —v?, c?) Wdx,, ys =Vsp, 2p =—ZB; 
tp—ts=(v/c?)(xe—x)) ; 

_ 2 (xe—xy)? dJ 

Kae! —2y) = J”) -—_————, (16.4) 


2c? r dr 


. 


where quantities without suffixes refer to system 
K. Consider a wave function y,,=¢,x, where ¢ 
depends only on 1 and x depends only on 2. Then, 
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using the conservation theorem for particle 
density, it is found, on neglecting all terms of 
order higher than v*/c*, that 


° 


(v'*x')o=| (1 + )orn+ (x*a'x) 
2c? Cc 


v O(x*a*x) 
t-(X2—- x1)— - , Che. 


ri 0x2" { 


ws 


where (a?, a”, a?)=(a', a*, a’) and (x, vy, 2 


=(x', x*, x’). The symbol (x*Yx) stands for 
xu* Y,.x.. From now on the suffix A is omitted 
where no confusion is caused. From Eqs. (16.4) 
and (16.5) it follows that 


J (xu! J(a0!—11))x0!)d Vs" 


wx (ras)dJ v* x* dJ 
=f (x|/- - — - |x)avs 
c vr dr 2 r dr 


and 


F vx (Tae) dJ 
frown = f (ofa Zo 
B ¢ t @ 


v x*dJ v? 
‘ v)a Vo, 


= - -J+ a;7J 
where r= (x, y, 2) =(x2—%X), Ve—V1, 22-21). (16.6) 





2¢ r dr 2¢ c 


Using Q of Eq. (13.3) 


| fwonars] =f wowars 


° 
¥ 4 ¢ 


v Vv Vx dJ 
+f(-| (a7 +as7) J+ J-- ((@1 + @2)T) 
2c 2c? 2¢ rdr 


vx? dJ 
Pts, 
2c? rdr 


On the left side the ’ means that a// quantities are 
evaluated in K’. Combining this with Eq. (16.6) 
and using the transformation formula for ¢*¢ it 
is found that 


| fw*ts—ewwyars] 
-[(v 


¢ 
Vv 


re) | 
x 


2c a 





= f wL7-owars. (16.7) 
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Since dV,dt is Lorentz invariant, so is Eq. (15). This equation does not take into account Q and 


The use of products ¢x for ¥ corresponds to small 
accelerations. Reducing Eqs. (16.1), (16.3) it is 


represents only the result of using Eq. (17). The 
expression for Q given by Eq. (16.3) is a general- 
found that if ization of a similar term? in the energy for electro- 
(E+J-+coyp:tcasps+(8:+8:)Mct}y=0, (17) magnetic interactions between point charges. In 

‘ _ : ‘ the latter case J = +e?/r and the term in Q repre- 
sents the combined effect of the magnetic energy 
and the retardation in the electrostatic potential. 
The term in Q has to be handled with circum- 


then the so-called ‘‘large’’ components of satisfy 
an equation of the Pauli type. This equation 


contains a non-Hermitean operator and is of the . 
spection because it represents the last term of an 


form 
expansion. It is known” in the electromagnetic 
Ey+ I ro J 2 case that that this term is only good for first- 
ied pe Pi~Pr- order calculations of the energy and that such 
: calculations lead to sensible results using directly 
+ p2Jp2—p2*J J¥+---=0, Eq. (16.1) and regarding it as a perturbed Eq. 


hoe ; (17). It may, therefore, be expected that the 
all other terms containing only Hermitean oper- , . . 
oe : : 3 ; same procedure applies to the Q of Eq. (16.3). 
ators. The occurrence of a non-Hermitean part is -. : . 
rhe result of the calculation may be expressed in 


analogous to that already encountered in Eq. (8). ‘ , ' 
_ : * aay “, terms of the Pauli function ¥™. If 
Analogously to Eqs. (8’), (9) it is convenient to 
transform the wave function by — 
' Vw } uy; Wut T 
yo = W+(p,°+ po”) Vv 8 .M?c? (17.1) 
-fv FY as: a VM agdr (17.2 
soas to have the density represented by (W()* Ww" 
The two-component wave function ¥“ satisfies (u, v, uw’, v’) =1, 2, 3, 4; (a, B, a’, B’) =3, 4 
os as h?AJ at p;* pi’ then the operator Y” may be called the equivalent 
- ~2 oh oh. pa ake . - aa . a . 
ie 2Mc?+ J 4M%e? IM 8 Mie Pauli operator. Using it with the Pauli function 
a @4 U- c* . - . . . 
in first-order energy calculations is equivalent to 
h : _ using Y with the Dirac functions satisfying Eq. 
4. = 7’ on a . : : 
Mie 2 LV iJ X pi joi) ¥ 0. (17) (17). Since the two parts of Q entering Eq. (16.3 
' occur also in other possible equations, their 
(AJ =A, J=A2J). equivalent Pauli operators are given separately : 
i. 1 | 2h ; dj 
4(a;a2)J |? = 4] p,po+—fr(po— p,) +h?(3 —2e,02)f+h?{ r? —r°e,02+(re,) (ree) | 
8.M?c? 1 rdr 
+ 2hf((r X pe joi — [Lr X ps jos) | , (17.3) 
dj" 1 8h 2h rdf 
— 3(a,r) (aor) = ——_—_) 4fx°xp,"p2° +—fr(P2— p1) +— r(p2— pi) 
rdr 8 M2? | 1 t dr 


— 2hf(L 1 X po joi — [tr X pi joo) +h2f(12+ 20,02) 


dsdf df 
+h -) +0 [9r? + re.02— (t01) (rez) J). (17.4) 


dr\rdr rdr 


where f=dJ/rdr, t=%—fe=(x!, x?,.x°). (17.5) 


A short calculation shows that the operators of Eqs. (17.3) and (17.4) are Hermitean. The above 
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expressions can be expressed more elegantly by means of symmetrical formulas involving J, p:, Po, 
pir), etc. It was thought, however, to be more useful to give the expressions in the above form 
because the occurrence of pi, pz only on the right, uniquely defines the expressions. Combining Eqs. 
(17’), (17.3), (17.4) one has as an equivalent Pauli equation in ¥™, 


‘ - pi? +p2? pii+ po hf 
—E+2Mc?—J+———--— 





+——([rX (2p2—p:) Jo.— [4 X (2pi— pz) Jor) er 


J 








2M 8M%c? = 4M?*c? {*c* 

f , F Shf h_ rdf h*f 

Ke PE oe Pd Se 
2M? «AM *C? " 4iMYe? dr 8 Mt? 
he df hr?’ d/df 
K (15+40,0,) -——— —[5r?+r°e,02— (re:)(re2) | —-———- — (—) w%=0. (17.6) 
4M°*c? rdr 8M*c? dr\ rdr 

For J=-—e*/r this equation simplifies into Eq. form (16.3) and not directly through the require- 


48) of reference 13. In that case terms in oj, 


(o:r)(@er) combine to give a contribution to the 


energy 
(he, 2Mc)*[ oor ~3 3(re;)(re2)r-* ], (1 7.7) 


which represents the magnetic interaction energy 
between spins having each a magnetic moment 
he/2Mc. The interaction energy of Eq. (17.7) 
vanishes to the first order for the deuteron. On 
the other hand" Eq. (17.6) contains terms which 
do not vanish for the S states of the deuteron and 
give a net effect in e,02 of amount 


— (h?/6M°*c*)(e,02)AJ (H?). (17.8) 


This interaction gives different energies in the 
singlet and triplet states, which corresponds to 
experience. The sign of Eq. (17.8) is such that for 
simple attractive potentials varying smoothly 
and monotonically with r, the energy is higher for 
triplets than for singlets. Thus ordinary, non- 
exchange, purely attractive interaction energies 
of the type of Eqs. (16.1), (16.3) are not likely 
to correspond to experience even though they 
contain a spin-spin dependence of the interaction 
energy. A repulsion in the region of smaller r 
produces a reversal of the sign of the effect of the 
expression (17.8). It is conceivable that such an 
interaction corresponds to reality, particularly 
since a combination of repulsions at small dis- 
tances with attractions at larger distances is 
helpful for securing stability of heavy nuclei 
without the aid of exchange forces. 

The presence of terms in oie, and (e;r) (esr) is 
secured above through the introduction of the 


ments of invariances. A form such as that of Eq. 
(17.6) remains invariant if arbitrary terms in 
spin-spin interactions are added to the left side, 
provided these terms are of order v*/c?. There is 
thus nothing binding about (17.8). The only 
argument advanced here is that of the attractive- 
ness of using a simple Diracian form such as that 
of Eq. (16.3). The verification of invariance given 
by Eqs. (16.4) to (16.7) is now seen to contain 
little more than a direct inspection of Eq. (17.6) 
and its comparison with possible spin-orbit inter- 
actions of Eq. (15.7) as well as possible classical 
forms given by Eq. (13.2). It will thus not be 
necessary to go through a formal proof in 
discussing Diracian forms 

The use of Eq. (16.2) is somewhat analogous to 
that of Eq. (10). It should be noted, however, 
that the term in 8:82/J without Q does not give an 
invariant contribution to the Lagrangian.” This 
can be verified either in the Dirac or the Pauli 
form. 

Reducing Eq. (16.2) with Q=0 to the Pauli 
form one finds 


| Sn cae 22 
E-2Me?+J+ |] -—-+—— 
| iL 2M 8M*%? 


p?J+2piJpit+ Jp? 


h 
—- Ce.xw|| w=0. (18) 
4M%c2 


Comparison with Eq. (15.4) shows that the spin- 
orbit interaction is given correctly by Ea. (18). 
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Comparison with Eq. (13.2) shows that the sum 
of the coefficients of p:°J, piJpi, Jp:* is twice 
what it should be. The equation can be corrected 
for invariance either by arbitrarily modifying 
these terms or by looking for Diracian forms 
which will automatically give a correct behavior. 
By means of Eqs. (17.3), (17.4) such forms can be 
reduced to equations of the Pauli type. 

It was seen that the classical Eq. (14.2) has a 
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simpler physical interpretation than Eq. (14.1 
A Diracian form agreeing with Eq. (14.2) in the 
sense of the correspondence principle is 
| E+c(eipi)+c(a2p2) + (61+82) Mc? +8182/ 

T $(a:a2)J+ 3 (air) (aor) (dJ rdr)\y =(0. (18.1) 


Its equivalent Pauli form is found to be, by 


means of Eqs. (18), (17.3), (17.4), (17.5), 


, 2 p:* D3 p?J+2piJpit+Jp?? h - J 
—~E+2Me!— J+ Fj) ——-——__+ ——— +——T[v.J x p; ]o; ]————(p:p>) 
| im1L2M 8%? 8 M22 4M22 2M%c? 


41M?c? dr 


£ivi “C 


Here the spin-orbit interactions are entirely of 
the Thomas type and there are no spin-spin 
interactions. Thus Eq. (18.1) may be looked at as 
a generalization of Eq. (10). It represents an 
interaction in which particle 2 is moving ap- 
proximately in the field of a scalar due to particle 
1. Eqs. (17.3), (17.4), (17.6), (18.2) apply only 
to interactions of a nonexchange type. 


4. APPLICATIONS 


According to Rumbaugh and Hafstad'* there 
is a level of Li’ at about 400 kv above the 
normal level. This level is probably the ?P, 
part of the *P normal level of Li’. According 
to Hafstad and Tuve" the resonance of the reac- 
tion C!®?+H!—N'"-+-hy shows indications of fine 
structure of the order of 80 kv. According to 
Herb, Parkinson and Kerst'*® there are indica- 
tions of fine structure in the resonances for y-ray 
production of Li, B, F, Al, Na, Be when these 
elements are bombarded by protons. The order of 
magnitude of the difference between adjacent 
resonance levels is about 200 kv. The spin-orbit 
interactions must be arranged so as to account 
for this order of magnitude. Using ordinary 
interaction potentials the order of magnitude of 
the expected splitting of a single particle moving 
in a p orbit is 


(19) 


- rc 
i ra [Jett Jav, 
M*c?r,? 


AE=— 


3h? =") h? 
4M*c* av 


rdr 


f h ; rdf 
~——— ae (37+ )x(p:—p))- 


h? rdf rd/df 
—| 157+10 +— (—) |} v=o. 
8M*c? dr dr \rdr/ I! 


(18.2) 





where the subscript av indicates the average 
value. Here J.¢; is the effective value of the inter- 
action energy due to all the other particles con- 
sidered as producing an effective central field 
The nuclear radius is ro. For [Jers ay =10 Mev 
and ro=e*/mc* the expected AE is then of the 
order of 50 kv. This agrees with the observations 
on N?!* but is too small for the 400 kv splitting in 
Li’. These are fitted better with [ Jes: Jay = 20 Mev, 
ro=e?/2mc* which is not unreasonable. Eq. (19 
corresponds to the use of Eq. (15.4) with 6=1, 
because in such a case the spin-orbit interaction 


consists entirely of terms of type 


[vidJiXei:)p; and > Ji; 


i+ i+i 


gives directly the potential in which particle 7 is 
moving. Thus Wigner forces are readily recon- 
cilable with present experimental indications as 
to the magnitude of fine structure using essen- 
tially interactions of the type of Eq. (18.1). The 
inversion of the fine structure in Li’ is also in 
agreement with this type of interaction energy. 

For Majorana forces Eq. (15.7) with a= 1 gives 
an interaction which is somewhat analogous to 
the Wigner case, just discussed, as is seen from 
Eq. (15.8’). With this form the following con- 
tributions to the fine structure are found for a 
proton interacting with a complete shell of 
2(2L’+1) neutrons 








he 
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(2ZL+1)(2L’+1)h? 
AE=-— — J (2@+1) 008 6 
327? M?c?L(L+1) 


xX P1[3Q:102+1101'Q2+1r20102" ] 

+sin? 6P,'[(4+L(L+1))Q:02+27r101'Q2 

+ 272010.’ +rir2Q1'Q2' ]} 
XRiv(A)Ri(2)PiuJidridte. (19.1) 


Here L =azimuthal quantum number of proton; 
L’'=azimuthal quantum number of complete 
shell; P; = Legendre function of argument cos 8, 
where @ is the angle between the directions of 
points 1 and 2; Rz-, rQ are the radial functions of 
proton and neutron respectively. Derivatives 
with respect to r and cos @ are indicated by ’. The 
normalization is such that 


fRerar=1 f erar- 1. (19.2) 


Eq. (19.1) applies only to cases in which a proton 
is outside a complete shell. It does not take care 
of perturbations due to other states of the shell, 
which, according to Feenberg and Wigner,!® are 
important in nuclear structure. The effect of such 
perturbations is not taken into account in this 
paper. 

For Li’ estimates of expected fine structure can 
be made by supposing that there are two protons 
in s states and one proton in a p state. Two 
neutrons are supposed to be in s states and two 
in p states. The pairs of like particles in s states 
form complete shells. Their effect can be com- 
puted by means of Eq. (19.1). The two neutrons 
in p states are considered here to be coupled into 
a'S term. They effectively form a complete shell 
of two particles. Eq. (19.1) applies to their effect 
on AE, provided the result of the calculation is 
divided by 6/2=3. This can be verified by ex- 
plicit calculation. In addition there is an inter- 
change effect due to the fact that the proton in 
the p state is a particle identical with the two 
protons in the s shell. In the estimates, the radial 
wave functions of all particles with the same L 
were taken to be the same, independently of 
whether they are protons or neutrons. The radial 


functions used were such that 


ite (19.3) 
Q=N ge -" 2 NQ?=8r5 2/3 ae 


and the interaction potential was taken to be 
J=Ae~™™. (19.4) 


The contributions to AE were found to be: 


3h? A (uv)*/*(u+4a) 
(AE),,= —— (19.5) 
4M*c? [a(ut+v)+(ut+v)?/4 >? 
due to the s shell of neutrons, 
3h? A(uv)*/"a@ 
(AE),.= (AE), -—— - (19.6) 
2M?c? [a(ut+yr) + py }>/? 
due to the s shell of protons, and 
h? Aav'(4av—3v?) 
(AE) »=-— —— (19.7) 


4M?*c? 2av+v? }/? 


due to two p neutrons in a 'S configuration. 

Using A =38 mc?, a=16 Mmc*/h®, and trying 
to fit the apparent experimental value AE = —0.8 
mc*, one obtains approximate agreement using 
p=v=40 Mmc*/h*. No precise fit was attempted. 
The above values of u and » correspond to values 
of (r?),y' of 1.7 10-8 cm and 2.2 X%10- cm for s 
and p particles, respectively. The maximum in 
the number of p particles per unit thickness of a 
thin spherical shell (per Ar=1) would be at 
2.0107" cm. Since the calculation is not pre- 
cise, these numbers may be considered satis- 
factory. It may be premature to expect an exact 
agreement inasmuch as even in atoms calcula- 
tions of fine structure using Hartree and Fock 
fields do not always give correct results. 

The values of » and » are seen to agree quite 
well with the calculations of Feenberg and 
Wigner,’® who obtain from mass defect values 
of u~=v=2ac of 26 and 32 in units Mmc*/h’. 
Increasing the value of A from 38mc* to 63mc* 
the value of u= v, as estimated from fine structure 
of Li’, drops to 27 Mmc*/h?. The agreement thus 
obtained may be partly accidental. 
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The first of a series of mass-spectrograph studies of the 
products of ionization of polyatomic molecules by low 
speed electrons is described. The following are the positive 
ions detected in methane together with their appearance 
potentials in volts: CH,* (13.1+0.4); CH;* (14.4+0.4); 
CH,* (15.7+0.5, 22.940.8); CH* (23.340.6); C* (26.7 
+0.7); H* (22.740.5, 29.4+0.6); H.* (27.940.5); and 
H;* (25.341.0). The following negative ions were also 
observed: CH.~; CH~; C~; and H~. From the higher 
appearance potential of H* and the appearance potentials 
of C* and C~ three independent upper limits are calculated 
for the heat of dissociation of CH, to normal atoms which 
lead, respectively, to the following three upper limits for 
the heat of sublimation of carbon: L(C)=6.1+0.6 volts 
or 141+14 kcal.; L(C)=5.7+0.7 volts or 131416 kcal.; 
L(C)=6.0+0.9 volts or 139+21 kcal. Thus of the three 
values of L(C) recently proposed from data on predissocia- 


INTRODUCTION 


NE of the important physical tools for 

studying molecular structure is the mass 
spectrograph wherein ions are formed by impacts 
of electrons of known energy with molecules. 
With this instrument one can determine in most 
cases what ions are produced from molecules of a 
particular variety through approximate measure- 
ments of the mass-to-charge ratios of the ions. 
Also, one can determine for ions of each kind: 
the relative probabilities for formation under 
different conditions in the ion source; the 
appearance potential or the minimum energy 
required to produce them; and, possibly, the 
kinetic energies with which they are formed. The 
first two types of information are of interest 
chiefly in themselves and can seldom be com- 
pared with results obtained by other means either 
experimental or theoretical. The appearance 
potential of undissociated ions or the ionization 
potential of the molecule under investigation is 
also important per se and often checks or supple- 
ments spectroscopic or theoretical results. In- 


1 The results on positive ions formed in methane were 
discussed briefly in a paper by W. Bleakney, E. U. Condon 
and L. G. Smith which was presented at the Symposium on 
Molecular Structure held at Princeton, N. J. under the 
auspices of the American Chemical Society (Dec. 31, 1936- 
Jan. 2, 1937) and which has been published in the Journal 
of Physical Chemistry. 
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tion in CO, these results are inconsistent with the value 
L(C) = 168.8 kcal. but are quite consistent with the value 
L(C)=123.6 kcal. Correlations of known ionization 
energies with the best obtainable data on the relative 
energies of the possible stable states of a carbon and four 
hydrogen atoms give reasonable interpretations of all 
observed appearance potentials. Upper limits for the 
ionization potentials of CH;, CH», and CH are deduced. 
These are: J(CH;)=9.9 volts; J(CH.)=12.0 volts; 
J(CH)=11.7 or 16.2 volts. Two general conclusions are 
drawn. The first is that most ions are only formed when 
the parent molecules are broken into as many parts as 
possible. The second is that interpretations of all but four 
appearance potentials are possible on the assumption 
that the total excitational and kinetic energy of products 
formed in an electron impact from a CH, molecule be 


less than two volts. 


terest in the appearance potentials and kinetic 
energies of other ions lies not so much in the 
quantities measured as in the dissociation and 
ionization energies which can in some cases be 
calculated therefrom by use of data from other 
scources. Since these ionization and dissociation 
energies comprise the most important contri- 
butions of a mass spectrograph to the field of 
molecular structure, a brief review of how they 
are obtained and of the chief sources of uncer- 
tainty involved in their calculations will probably 
not be out of place here. 

The appearance potential of an ion which lacks 
one or more of the atoms of the parent molecule 
is the sum of three quantities: a dissociation 
energy; an ionization energy; and the total 
minimum excitational and mutual kinetic energy 
with which the products can dissociate. (The 
ionization energy in the case of a negative ion is 
an electron affinity—with negative sign—and if 
both a positive and a negative ion are produced, 
it is the algebraic sum of an ionization energy and 
an electron affinity.) Provided one knew what all 
the products of dissociation are, how much 
kinetic energy each has, and in what excited state 
each is formed, he could calculate the dissociation 
energy from the appearance potential and a 
knowledge of the ionization energy or vice versa. 
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Unfortunately, if more than one atom of the 
parent molecule is missing in the ion, there are at 
least two possibilities for the states of aggregation 
of the missing atoms, and the reaction occurring 
can only be determined through considerations of 
energy relations. Furthermore, the kinetic ener- 
gies of unionized products of dissociation cannot 
be calculated from that of the ion if more than 
one of the former is produced from a molecule. 
Finally, a priori knowledge of the states of 
excitation of the products is never at hand. 
Conclusions drawn from combinations with other 
data of the appearance potential and kinetic 
energy of an atomic ion formed from a diatomic 
molecule are subject only to uncertainties in the 
electronic states of the products of dissociation. 
In the case of the appearance potential of an ion 
lacking but one of the atoms of a parent poly- 
atomic molecule, conclusions are also subject to 
uncertainties with regard to vibrational energy 
of the ion. In general, where there are more than 
two products of dissociation of a single molecule, 
not only electronic and vibrational energy but 
kinetic energies of the products as well as their 
states of aggregation are uncertain. Here the 
best one can hope to do is to determine the 
reaction occurring and from this to determine a 
value of either the dissociation or the ionization 
energy involved plus the total excitational and 
kinetic energy of the products. It is evident that 
such a value is an upper limit for the true 
dissociation or ionization energy (or a lower limit 
for the true electron affinity). Additional factors 
complicating the interpretation of results ob- 
tained by the electron impact method are: 
experimental errors in these results (which are 
seldom less than 1 percent and often as high as 4 
or 5 percent); and uncertainties in quantities 
obtained from other sources. The latter uncer- 
tainties, however, add importance to any con- 
clusions which may be drawn even if it is 
possible to determine only upper or lower limits 
for the true values. 

The results for diatomic molecules obtained by 
the electron impact method have in many cases 
corroborated and supplemented spectroscopic 
and thermochemical data. Also they have given 
experimental evidence in support of the Franck- 
Condon principle and the repulsive electronic 
states of such molecules predicted by quantum 
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mechanics. Results for the triatomic molecules 
that have been studied (H2O, HS, N2O, NOs, 
SO2, COz and CS.) have also in some cases 
checked and added to thermochemical and 
spectroscopic data. However, little attempt has 
been made to interpret in the light of other data 
results for the more complex molecules which 
have been studied (NH3, CeNe, CH, CH;Cl, 
Cs3Hs, CsHi9, CsHis, CoHe, CsHs). Good reviews 
of all this work have been prepared by Smyth,’ 
and more recently by de Groot and Penning,’ and 
by Sponer.* 

The mass spectrograph analysis of methane, 
results of which are reported in this paper, is the 
first of a series of such analyses to be undertaken 
in this laboratory. The results for methane yield 
some information regarding unknown and un- 
certain quantities. It is to be hoped that with the 
accumulation of further data for other molecules 
more such information will be obtained. 


APPARATUS AND PROCEDURE 


The mass spectrograph used to study ionization 
processes in methane is of the type described by 
Bleakney® and by Tate, Smith, and Vaughan.* 
A detailed description of the instrument, the gas- 
handling system, the electrometer-tube amplifier 
used to measure the ion current, and the device 
for recording automatically this current as a 
function of the potential applied to accelerate 
either the ions or the electrons will be found in a 
recent paper.’ 

Fig. 1 shows a typical mass spectrum obtained 
with the automatic recorder of the positive ions 
with m/e values from 12 to 16 formed in methane 
by electrons with an energy of 50 electron volts. 
When this record was obtained the magnetic field 
was kept constant at about 550 gauss (as in all 
work on methane), and the potential applied to 
the electron gun was fixed at 50 volts, while a 
variable potential accelerating the ions was 


applied by the drum potentiometer which is a 

2H. D. Smyth, Rev. Mod. Phys. 3, 347 (1931). 

3 de Groot and Penning, Handbuch der Physik, Vol. 23/1 
(J. Springer, Berlin, 1933). 

4H. Sponer, Molekulspektren (J. Springer, Berlin, 1935, 
1936). 

5 W. Bleakney, Phys. Rev. 40, 496 (1932). 

J. T. Tate, P. T. Smith, and A. L. Vaughan, Phys. 
Rev. 48, 575 (1935). 

7™W. W. Lozier, P. T. Smith, L. G. 
Bleakney, Rev. Sci. Inst, 8, 51 (1937). 
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Fic. 1. Mass spectrum of the heavier positive ions formed 
in methane. Ordinates represent ion currents and abscissae 
analyzer voltage or e/m values. I~ is the electron current 
recorded simultaneously with the ion current. The small 
peak to the left of the CH,* peak corresponds to ions with 
m/2=17 which can only be CH,* ions. The ratio of the 
heights of these two peaks is ~1/100 which is approxi- 
mately the abundance ratio C4/C®, 


part of the recorder. The height of a peak is 
assumed to be proportional to the intensity of the 
corresponding ionization or number of ions of the 
given kind formed per second per unit volume.® 
The pressure of methane in the mass spectrograph 
when Fig. 1 was obtained was ~4X10-> mm 


of Hg.® 


’The deflection of the galvanometer is strictly pro- 
portional to the ion current arriving at the collector. This 
fact coupled with the fact that peaks for ions with m/e as 
high as 6 are flat topped indicates that even for ions with 
m/e=16 this assumption should not be far from the truth. 

* Pressures were estimated by calculation from the 
absolute values of electron and total positive ion currents, 
the ey of the ion source, and the value given by 
J. T. Tate and P. T. Smith (Phys. Rev. 39, 270 (1932)) for 
the cross section of common gas molecules for ionization 
which was assumed to be correct for methane. Such 
estimates of pressures should not be in error by more than 
factors of two. 
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Fig. 2 illustrates the method of determining 
the appearance potential of an ion (in this case 
H+ formed from CH,). Curve 1 was taken with 
the ion-accelerating potential fixed so that H* 
ions passed through the exit slit of the analyzer 
and with the variable potential from the drum 
applied in series with a fixed potential to ac- 
celerate the electrons. Curve 2 was taken in 
similar fashion with the analyzer set to record 
Het formed in helium gas which was 
admitted into the apparatus simultaneously with 
the methane. The curve for He* serves to correct 
the recorded electron energies for contact po- 
It will be observed 


ions 


tentials in the instrument. 
that the upward break of curve 2 from 
horizontal portion is more abrupt than that of 
and that to the right of the break, curve 2 
while curve 1 has 
10 Hence in de- 


its 


curve 1 
nearly a 
considerable upward curvature. 


straight line 


is 


} t +— 
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Fic. 2. Appearance potential curves for H* ions formed 


in methane and Het ions formed in helium. Ordinates 
represent ion currents. I~ is the electron current recorded in 
each case simultaneously with the ion current. 


10 This effect is characteristic of the appearance potential 
curves for all ions from methane and has been noted by 
others who have measured the appearance potentials of 
molecular ions. The effect is presumably due first to the 
fact that at room temperature a gaseous molecule may 
exist in one of a larger number of energy states than a 
gaseous atom, Secondly, it is due to the fact that a molecu- 
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termining the lowest appearance potential of H* 
ions it was not feasible to attempt to fit curves 1 
and 2 together, nor could one extrapolate to zero 
current any linear portion of curve 1 to the right 
of the upward break. It seemed best in this as in 
all other cases to use the method of selecting the 
appearance potential for each ion as the lowest 
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amplifier and the relative pressures of methane 


and the calibrating g 


curve for ions from methane had roughly 


yas so that in each case the 


the 


same slope to the right of its upward break as the 


curve for ions formed 
With the total 
5x<10-§ mm of Hg, 


pressure 


in the 
of 


the pressures 


calibrating gas. 
impurities about 


of methane 


electron energy at which the corresponding ion used in determining the appearance potentials of 
current was detectable. This method is open to positive ions ranged from about 5 X 10~* to about 
several objections, the chief of which is that the 510-5 mm of Hg, while the pressure was about 
value obtained depends somewhat on the sensi- 1X10-* mm of Hg when we determinations . 
tivity used in obtaining the corresponding curve. were made for negative ions.’ An equipotential, 
This source of error was reduced as far as_ oxide-coated cathode was used in obtaining all 
possible by adjusting the sensitivity of the appearance potentials. Because of the smaller 
eT range of the speeds with which electrons leave 
lar ion may be formed in one of a much larger number of = cuch a source. the upward breaks in curves like 
closely spaced energy states within a range of a few volts ens * 
above the lowest state. those of Fig. 2 were more sharply defined when 
TABLE I. Jons formed in methane. 
. 
© or TOTAL 
IONIZATION OBSERVED ADDITIONAL ENERGY 
(ELECTRON | APPEARANCE CALCULATED MINIMUM REQUIRED OR IONIZATION 
OBSERVED ENERGY POTENTIALS PROBABLE PROCESS ENERGY ENERGY DEDUCED 
Ion 50 VOLTs) (VOLTS) CH, VOLTS (VOLTS) ’ 
Positive ions 
H* 3.0 22.7+0.5 H*++CH.+H (1) 21.6 1.1 ‘ 
or Ht++CH+H, (2 20.6 2.1 
29.4+0.6 H*++C+3H (3 28.6 0.8 
H,** 0.3 27.9+0.5 H.+*+C+2H (4 26.0 1.9 
H,* * 0.005 25.3+1.0 H;*+C+H (5) 24.1— D(H3*) 1.2+D(H;*) 
or H;*++CH- (6) 20.6 — D(H;*) —E(CH 4.7+D(H;*)+£(CH 
0.6 | 26.7+0.7 C*++4H (7) 26.3 0.4 
CH* 1,7 | 23.3+0.6 CHt+3H (8) 11.6+J7(CH) I(CH)=11.7 
or CH*++H.+H (9) | 7.1+J(CH I(CH)=16.2 ‘ 
CH, 4.2 15.7+0.5 CH.*+H: (10) | 3.7+I](CH2) I(CH2)=12.0 
22.9+0.8 CH.*+2H (11) | 8.1+J7(CH: =2.8 
CH;* 39.5 14.4+0.4 CH;*++H (12) | 4.5+](CHs) I(CH;)= 9.9 
CH,* 50.7 13.1+0.4 CH,* (13 I(CH, I(CH,) = 13.1 
. 
Negative ions 
H | 0.04 6.1+0.3 H-+CH; (14 3.8 2.3 
| 6.9+0.8 H-+CH+H, (15 6.4 0.5 
7.3+1.0 H-+CH2+H (16) 7.4 ~0 
18.0+0.6 H-~+CH*t+H, (17) =18.1 ~0 
or H~+CH;* (18) =13.7 4.3 
0.003 27.4+0.6 C-+Ht+3H (19 26.7 0.7 
Ci 0.005 10.2 +0.3 CH-+3H (20 11.6—E(CH) E(CH)=1.4 
24.5+0.5 CH-+H*+2H (21 25.1—E(CH) “0.8 
or CH~+H,;* (22) 20.6—D(H;*) —E(CH 3.9+D(H;*)+E(CH 
CH, 0.002 8.9+0.3 CH,.~+2H (23 8.1—E(CH2) 0.8+E(CH2) e 
23.4+0.6 CH.-+H+t+H (24 21.6—E(CHz2) 1.8+E(CH2) 
or CH,~+H,* (25 19.0—E(CH,) 4.4+ E(CHz2) 
* The results for H2* and Hs;* ionization in methane are of interest in connection with attempts to produce helium of mass three by bombard- ’ 
ment of deuterium with canal rays from an intense, high voltage discharge in deuterium. (Phys. Rev. 46, 81 (1934); Phys. Rev. 47, 800(A) (1935 


Samples of the treated gas were allowed to 


stand over hot copper oxide 


and a liquid-air trap to remove the hydrogen in the form of water. In 


the residual gases admitted to the mass spectrograph later used for obtaining the present results on methane, 


ions were found with 


m/e=3 whose 


appearance potential was approximately equal to the ionization potential of He (24.5 volts). These were at first believed to be He** ions, but it 
was later found that the rate of diffus sion through 
It was then suggested that the ions might be formed from some hydrocarbon gas like methane. 
probably not HD* ions formed from deuteromethane 
have been H;* ions from ordinary methane. 


since their 


appearance potential was too low; 


a capillary of the gas from which they were formed was smaller than that of ordinary helium 
The present results indicate: 
and (2) that only very few of them could 


(1) that they were 
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such a source was used than when a tungsten 
filament supplied the electrons. 


EXPERIMENTAL RESULTS 


A list of all ions, both positive and negative, 
which were observed, their approximate relative 
intensities for an electron energy of 50 volts, and 
their appearance potentials will be found in the 
first three columns of Table I. The errors quoted 
in column three are estimated probable errors. 
Helium served as the calibrating gas in deter- 
mining the appearance potentials of the negative 
ions and H* ions, while argon and a slight 
impurity of hydrogen were used in the cases of 
the other ions. In all cases the calibrating gas 
was in the instrument along with methane. 
No multiply charged ions could be detected at 
electron energies as high as 80 volts. 

The positive ions were all observed when either 
a tungsten filament or an oxide-coated cathode, 
operating at considerably lower temperature, was 
used as an electron source. Their relative in- 
tensities were the same within experimental error 
in the two cases indicating that few, if any, ions 
were formed by ionization of products of thermal 
dissociation at either cathode. In Table II are 
given the results of three independent measure- 
ments of the relative intensities of the heavier 
positive ions showing the constancy of these 
quantities at widely different pressures and 
electron currents. Constancy of the relative 
intensities of H;3+, H2+, and H?* ionizations at 
different pressures is shown by Fig. 3 in which are 
plotted the intensities of H;+ and H.* as func- 
tions of the intensity of CH,* and the intensity 
of H+ as a function of that of H2+. Measurements 
of the intensity ratios H;+/CH,*, He+/CH,*, and 
H*+/CH,* at electron currents ranging from 3.1 
to 60, from 0.52 to 32, and from 0.52 to 60 
microamperes, respectively, gave results which in 


TABLE II. Percents of total ionization at various pressures and 
yj 
electron currents. 





APPROXIMATE ELECTRON 
PRESSURE CURRENT | 

(MM OF Hg) (MICROAMP.) | ¢ H.* | CH:* CH:2*|CH*t; Ct 
1x 10-6 0.53* | 52.3) 39.2) 3.4 | 1.3 | 0.6 
4x 10-5 4.3* |52.0| 38.8] 3.8 | 1.6 | 0.6 
2x 10-6 61.07 49.0/ 41.0} 4.3 | 1.9 | 0.7 


* Electrons supplied by oxide-coated cathode. 
t Electrons supplied by tungsten filament. 
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Fic. 3. Intensities of H;+ and H,* as functions of CH,* 
and the intensity of H* as a function of that of H.*. The 
linearity of these curves shows that H*, H,*, and H;* ions 
are formed from CH, molecules as primary products. 


no case disagreed with the ratios obtained from 
the curves of Fig. 3 by more than a factor of two. 
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From all these facts we may be sure that each 

~ positive ion is formed as the result of a single 

impact between an electron and a CH, molecule. 

The heights of the negative ion peaks were not 

measured as functions of pressure or electron 

current, nor were these peaks looked for when a 

tungsten filament was being used as a source of 

electrons. However, the absence of secondary 

processes in the formations of any of the positive 

ions together with the fact that the ion current 

vs. electron energy curves for all the negative ions 

observed in methane (Fig. 4) are of the same 

character as those obtained by various investi- 

gators for other negative ions which are formed 

in single electron impacts leads one to assume 

tT i that all negative ions from CH, molecules are 
formed in the same manner 


[os ee a . ~ 
COMPARISONS WITH OTHER ELECTRON IMPACT 
RESULTS FOR METHANE 








Th | | The ionization potential of CH, was measured 
a | prior to 1929 by several workers using instru- 
| MmLe | ments which could not separate ions according to 

-4 = e their mass-to-charge ratios. These results are 

~~ ° . om re 

ell aes OE summarized in Table III. Using a mass spectro- 

graph, Hogness and Kvalnes" found that when a 
tungsten filament supplied electrons CH,*, CH;*, 
l | | CH,*, CH*, C+, H.*, and H* ions were produced 
=>, ao in methane, but that with an oxide-coated 
‘wae cathode as source of electrons only CH,* and 
sae a0 Sie ae : d : 
= oe oon CH;* appeared in about equal numbers at 
: +» ' a ° . - - 2 
| = minimum electron energies of 14.5 and 15.5 volts, 
r | \ respectively. They therefore ascribed the forma- 
wa 
| | _ — , a a 
| raB_eE III. Early results on the ionization potential of CH. 
) } 
| | | A 
| IONIZATION 
=| | 1‘ -OTENTIAL 
| —- + OBSERVER REFERENC! VoLts 
| | } | | —- _ - — — - — - 
-— | | Mayer Ann. d. Physik 45, 1 (1914 13.5 
! Hughes and Dixon Phys. Rev. 10, 495 (1917 9.5 
| | Hughes and Klein Phys. Rev. 23, 450 (1924 13.9 
| —~h1_ ietsch Dissert. Berlin p. 534 (1926 12.5+0.5 
| Glockler J. Am. Chem. Soc. 48, 2021] 14.4 to 15.2 
1926 
Pietsch and Wilke* Zeits. f. Physik 43, 342 (1927 14.58 +0.05 
Morris Phys. Rev. 32, 942 (1928 14.4 





* Pietsch and Wilke also found a second critical potential for ior 
tion at 15.40+0.05 volts which was later interpreted by Hogness ar 
Kvalnes, reference 11, as the appearance potential of CHs?*. 
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et | ee) eS Fic. 4. Ionization efficiency curves for negative ions 
| formed in methane. Ordinates represent ion current. I~ is 
the electron current recorded in each case simultaneously 
a ee a ee 1 io . - : > 
with the ion current. The electron energy increases from 
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tions of the other ions to ionization of the 
products of thermal dissociation at the tungsten 
filament. Hipple and Bleakney,” analyzing 
methane with a mass spectrograph of the same 
type as, but less sensitive than, that used in the 
present work, observed the ions CH,*, CHs*, 
CH,*, CH*, C+, and H* at the minimum electron 
energies 13.7, 14.7, 15.7, 23, 27 and 31 volts, 
respectively. Pressure, electron current, and 
pumping speed variations showed these all to be 
primary products. In addition, an H2* peak was 
observed about 10~ as high as the CH,* peak. 

Because of the use of lower pressures made 
possible by improved vacuum technique and 
greater sensitivity of the device for measuring ion 
currents, it is believed that the author’s values 
for the appearance potentials of CH,* and CH;* 
are more reliable than those of Hogness and 
Kvalnes or those in Table III. All the results of 
Hipple and Bleakney agree within experimental 
error with those of the author excepting their 
values of the relative intensity of H,* ionization 
and the appearance potential of H* ions. Evi- 
dently H* ionization in their apparatus was too 
weak to enable them to detect its appearance at 
22.7 volts, and their 31—volt value probably 
corresponds to the author’s value of 29.4 volts 
for the second appearance potential of H* ions. 

The value 12.72+0.05 volts for the ioniza- 
tion potential of CH, reported by Smith and 
Bleakney™ was determined by a method to which 
several objections have been raised. The value 
given in Table I was obtained by the method 
described above. 


INTERPRETATION AND DISCUSSION 


In the discussions below the following abbrevi- 
ations will be used: 


A(X*) =Appearance potential of X* ions formed 
by electron impact in methane; 
I(X) =Ionization potential of X; 
E(X) =Electron affinity of X; 
D(X)=Heat of dissociation of X to normal 
atoms in the gaseous state; 
Q(X) =Energy evolved in the formation of a 
gram mole of X from its elements in the 
8]. A. Hipple, Jr. and W. Bleakney, Phys. Rev. 47, 
802(A) (1935). 


3L.G. Smith and W. Bleakney, Phys. Rev. 49, 883(A) 
(1936). 


states in which they exist at room tem- 
perature and atmospheric pressure ; 
L(C) =Heat of sublimation of carbon (energy 
absorbed by a gram mole of carbon in 
changing from diamond to a gas con- 
sisting of atoms in the *P normal state). 


Here X stands for the chemical symbol of 
a molecule, a radical, or an atom. All these 
quantities will correspond to reactions taking 
place at the absolute zero of temperature.“ 
Whenever energy values have a direct bearing on 
thermochemical data they will be given both in 
volts (meaning electron volts per molecule) and 
in kcal. (meaning kcal. per mole). The conversion 


4 


factor is: 1 volt =23.054 kcal. 
STATE ENERGY 
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IS.1 VOLTS ‘348 KCAL. 
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CH,+2H 8.1 « 187 ° 
CH+H+H, - ff... TI 164 ° 
C+2H, £ 62 * :142 « 
CH,+H * 45 = :104 

CHg+Hg ——— 3.7 * 64 C* 
CH, s__ 0 - fe) e 


Fic. 5. Energies of the possible stable states of combi- 
nation in which a carboy atom and four hydrogen atoms 
may exist referred to the state CH, as zero. 


Fig. 5 shows the energies of the possible stable 
states of combination in which a carbon atom 
and four hydrogen atoms may exist, referred to 
the state CH, as zero. These energies are in 
accord with what seem to be the most reliable 
data available, though considerable uncertainty 


4 Dissociation and ionization energies obtained spectro- 
scopically, since they involve only considerations about two 
energy states of a single molecule, correspond to reactions, 
from the standpoint of thermochemistry, at absolute zero. 
According to the Franck-Condon principle an electron 
striking a molecule cannot change directly either the 
translational or rotational kinetic energies of the nuclei. 
Hence an appearance potential determined by the electron 
impact method is simply the energy difference between a 
vibrational level of an upper electronic state and one of a 
lower electronic state and is thus the minimum energy 
required for a reaction at absolute zero. 
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exists as to their validity. Levels c’, d’, e’, and e” 

must lie below levels c, d, e, and e’, respectively, 

by D(H2) =4.46+0.04 volts or 102.9+0.9 kcal.'® 

The position of the level e is determined by 

means of the following cycle: 

Caiamonat+2H2-CHy, +Q(CHs,) 

4H-—2H, +2D(H2) 
Cont diamond + L(C) 


4D C H,) 


Coast 4H —C H, 


This shows that D(CH,)=Q(CH,4)+2D(H2) 
+L(C), which yields 
294.5 9 keal. 
D(CH,)= 224.5 +1.9 kcal lac) (1) 


9.84+0.08 volts | 


on substitution of the above value of D(H2) and 
the value Q(CH,) = 18.66 kcal. or 0.81 volts.'® 
There is considerable uncertainty concerning 
the true value of L(C). Bichowsky and Rossini!’ 
after reviewing all data published prior to Jan. 1, 
1934 set the value at 170—a kcal., ‘“‘where a is an 
unknown whose value at the present writing may 
be as much as +60 kcal.’”’ Recently data on 
arog 9 in CO have set three possible 
values : L(C) = 168.8+0.5 kceal., 123.6+0.5 kcal., 
or sot ne kcal. These data cannot, however, 
decide which is the correct value since they do 
not reveal in what electronic states the C and O 
atoms are formed in a predissociation process. 
Goldfinger and Jeunehomme," using these values 
tocalculate vapor pressure vs. temperature curves 
for carbon, conclude that experimental curves, 
inaccurate as they may be, quite definitely rule 
value and seem to favor the 
their final conclusion is that 
value 


out the lowest 
highest. However, 
the balance of evidence is in favor of the 


L(C) =123.6+0.5 kcal. or 5.36+0.02 volts. (2) 
Adopting this value one obtains from Eq. (1) 


D(CHs) = 348.1+2.0 kcal. or 15.10+0.10 volts. 
(3) 


A125, 2 


18 Richardson and Davidson, Proc. Roy. 
(1929). 

16 This quantity was obtained from the value, reference 
18, for the corresponding reaction under constant pressure 


Soc. 


at 18°C by correcting for the thermal energies of the 
reactants. 

17 Bichowsky and Rossini, Thermochemistry of Chemical 
Substances (Reinhold Publishing Corporation, N. Y., 
1936). 

18 P. Goldfinger and W. Jeunehomme, Proc. Faraday 
Soc. 32, 1591 (1936). In this paper will be found an 


excellent bibliography of recent papers on this subject. 
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The location of level } in Fig. 5 at 4.5 volts 
above a is in agreement with the estimate of 
Norrish'® from thermal and spectroscopic data. 
It also agrees with the results of Trenner, 
Morikawa, and Taylor”? which show that the 
energy required for the reaction CH,—~CH;+H 
at constant pressure and room temperature is 
almost certainly greater by a small amount than 
the heat of dissociation of hydrogen measured 
under the same conditions. Correcting the latter 
quantity’ for the thermal energies of the re- 
actants CH, and CH;+H, one concludes that 
the level 6 must be slightly more than 102.4 kcal. 
or 4.44 volts above zero. Further evidence on 
this point is afforded by s*! theoretical 
calculations of the energies of formation, meas- 
ured from the minima of the potential surfaces, of 
CH,, CH3;, CHe, and CH from normal carbon and 
hydrogen atoms. From the normal frequencies” 
of these molecules the author has calculated their 
zero-point energies to be: 1.17 volts, 0.82 volts, 
0.46 volts, and 0.18 volts, respectively. If now 
for the energy of formation of CH, the value 
16.27 volts (quantity (3)+1.17 volts) is used in 
Voge’s equations, a new value for the energy of 
formation of CH; may be obtained which, when 
corrected by the zero-point energy of CH; and 
subtracted from (3), gives 4.40 volts or 101.4 
kcal. as the height at which level 6 should be set 


Voge’s 


above zero. 

Level d must be below level e by D(CH) for 
which the value 3.5 volts‘ has been adopted. 
Level c was located by adding the zero-point 
energy of CH: to Voge’s result for the energy of 
formation of CHz based on the above value of 
D(CH) and subtracting the result from D(CH,) 
(3). 


) Positive ions 

H*.—The dissociation energy plus the 
kinetic and excitational energy involved in the 
higher of the two appearance potentials of H* 
ions is 29.4—13.5=15.9+0.6 volts, where 13.5 
=I1(H). The dissociation process therefore 


19R. G. W. Norrish, Trans. Faraday Soc. 30, 106 (1934). 

20N. R. Trenner, K. Morikawa and H. S. Taylor (paper 
soon to appear). 

21H. H. Voge, J. Chem. Phys. 4, 581 (1936). 

* Frequencies of CH, and CH as well as those of the 
methyl halides from which those of CH; were estimated are 
given by Sponer, reference 4. The frequencies of CH: are 
given by Sutherland and Dennison (Proc. Roy. Soc. 
A148, 250 (1935)). 


total 


volts 





IONIZATION AND 
could be that represented by a transition from 
level a (Fig. 5) to any one of the higher levels 
except c’ and e’’. If this higher level is e the total 
excitational and kinetic energy which an H?* ion, 
three H atoms, and a C atom must have when 
formed is 15.9—15.1=0.8+0.6 volts; if it is any 
level lower than e at least 4.3 volts must be 
assigned to kinetic and excitational energy of the 
products of dissociation. It seems reasonable, 
therefore, to assume that the reaction wherein 
H* ions are formed at a minimum energy of 29.4 
volts is: CH,-H*++C+3H. Subject to the 
validity of this assumption, the following upper 
limits may be set: 


| D(CH,) =15.9+0.6 volts or 366+14 kcal. | 


| L(C) = 6.1+0.6 volts or 141+14 kcal. | 4) 


This conclusion, if correct, cannot be reconciled 
with the value 168.8+0.5 kcal.'* for L(C). 
Subtracting 7(H)=13.5 volts from 22.7 volts, 
one sees from Fig. 5 that the dissociation process 
involved in the formation of H* ions at 22.7 volts 
minimum energy must be the three 
represented by transitions from a to c, d’, or b. 
The total kinetic and excitational energy of the 
products in these cases would be: 9.2—8.1=1.1 
volts ; 9.2 —7.1=2.1 volts; or 9.2—4.5 =4.7 volts. 
Since the latter value seems high it is probable 


one of 


that these ions are produced either in the reaction 
CH,—-H++CH.2+H or the reaction CH,—-H* 
+CH+Ha2. 

H,*.—It from Fig. 5 that the 
dissociation energy involved in A(H¢*) is that 


is evident 


required for one of the reactions whereby a CH, 
molecule is transformed to one of the four states 
e’, d’, e’’, or c’. Since the dissociation plus kinetic 
and excitational energy is 27.9—15.4=12.5+0.5 
volts (where 15.4 volts =7(He)) it seems likely 
that the dissociation is to level e’ rather than toa 
lower level, and hence that H,* ions are formed 
in accordance with the CH,—H;+ 
+C+2H. The products must possess at least 
12.5—10.6=1.9+0.5 volts of kinetic energy or 
energy of excitation. 

H;*+.—If we let D(H;*) stand for the unknown 
energy absorbed in the reaction H3;+-+H,.+H* 
we have as the energy of the dissociation of CH, 
to a state involving three H atoms plus the total 
kinetic and excitational energy involved in 


A(H;+): A(H;+) + D(H;+) + D(H) — J(H) 


reaction 
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= 16.3+1.0 volts +D(H;*). In accordance with 
previous reasoning this state must be either d or e 
of Fig. 5. This would indicate that H;* ions are 


probably formed in the reaction CH,—-H;* 
+C+H, the energy required to excite or 
give kinetic energy to the products being 


{16.3+D(H;*)} —15.1=1.2+1.0 volts+D(H;°*). 
However, it appears that in the case of these ions 
one should consider the possibility that products 
of dissociation formed along with a positive ion 
may be negatively charged as well as neutral. A 
consideration of the negative ions (see below) 
shows that possibly an H;* ion is formed together 
with a CH 


be considered in the interpretations of the ob- 


ion (but that negative ions need not 


served appearance potentials of any other type 
of positive ion). 

C+.—There are three transitions a—e”’, a—¢e’, 
and a—6< (Fig. 5), the energy required for which 
involved in 
from 


could be the dissociation energy 
A(C*). Subtraction of /(C)=11.2 
26.7+0.7 volts yields 15.5+0.7 volts as the sum 
of this dissociation energy and any kinetic or 


volts 


excitational energy with which the products may 
be formed. It therefore seems most reasonable to 
assume that C+ ions are formed according to the 
reaction CH,—-C++4H, the total kinetic and 
excitational energy of C++4H being 15.5—15.1 
=0.4+0.7 volts. If this assumption is correct the 
observation A(C+) =26.7+0.7 volts leads to the 
conclusion : 


| D(CH,) =15.5+0.7 volts or 357 +16 kcal., | 
| L(C) = 5.7+0.7 volts or 131416 kcal. | 


Here again we have fairly strong evidence against 
the value L(C) =168.8+0.5 kcal." 
CH*.—There are two possible reactions which 
could account for the production of CH* ions 
from CH, molecules: (a) CH,—-CH*++3H; and 
(b) CH,—-CH*+H-+Hsg. If (a) is correct J(CH) 
= A(CH?*) — {D(CH,) — D(CH)} = 11.7 + 0.6 
volts (according as the products of dissociation 
are formed with excitational and 
kinetic energy), while if (b) is correct 7(CH)=11.7 
+D(H:)=16.2+0.6 volts. No data on J(CH) 
which might enable one to choose between (a) 
and (b) exist.” 
3 Ina previous paper, reference 1, the value 150 kcal. 
was calculated for L(C) from the difference between 


A(CH*) and the appearance potential of CH* ions 
formed in acetylene. The reasonableness of this value was 


or without 
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CH,*+.—Since there are two appearance po- 
tentials for CH2* ions and only two reactions 
whereby they may be formed it seems quite 
probable that 15.7+0.5 volts is the minimum 
energy required for the process CH,—>CH2*+H,: 
while 22.9+0.8 volts is that required for the 
process CH,—CH,.*++2H. Thus, from the ener- 
gies required for the corresponding dissociations 
without ionization (Fig. 5), we may conclude 
that J(CH2)=15.7—3.7=12.0+0.5 volts, and 
hence that the products of the reaction requiring 
22.9 volts must dissociate with at least 22.9 
— (8.1+12.0)=2.8 volts of kinetic and exci- 
tational energy. 

CH;*+.—Simple considerations show that 14.4 
+0.4 volts is the minimum energy required to 
produce CH;* ions according to the reaction 
CH,—-CH;++H. This means that 7(CH;)=14.4 
—4.5=9.9 volts according as the CH;* ions and 
H atoms are formed with or without energy. 
This is not in very good agreement with either 
Fraser and Jewett’s value of 11.1+0.5 volts*® 
obtained by direct measurement on free CH; 
radicals or with Mulliken’s theoretical estimate 
of 8.5 volts. 


(b) Negative ions 


The curves shown in Fig. 4 are characteristic of 
ion current vs. electron energy curves for negative 
ions formed in various substances. The interpre- 
tation which is at present accepted of such 
curves”’ is that a peak at low electron energies 
indicates the formation of a negative ion through 
dissociation of a molecule and capture of the 
impacting electron, while the rising of a curve to 
a broad maximum at higher energies corresponds 
to the dissociation of molecules into pairs of 
positive and negative ions. No explanation of the 
apparent upward breaks at ~47 volts, 50 volts, 
and 45 volts in the curves of Fig. 4 for H~-, C-, 
and CH-, respectively, has been attempted. It is 
taken to be evidence in favor of reaction (a). However, 
on the assumption that reaction () is correct and that 
CH*t ions are produced from C:H2 molecules in the 
process C;H;—-CH*+CH one obtains by a similar calcu- 
lation L(C) = 128 kcal. 

* The conclusions with regard to CH,* ions in a previous 
paper, reference 1, were less definite because of doubt as to 
the reality of the appearance potential at 22.9 volts. This 
doubt has since been largely dispelled. 

*°R. G. J. Fraser and T. N. Jewett, Phys. Rev. 50, 1091 
(1936). 

26 R.S. Mulliken, J. Chem. Phys. 1, 492 (1933). 

27 W. W. Lozier, Phys. Rev. 46, 268 (1934). 
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quite possible that these are due to some peculi- 
arity of the apparatus and not to the formations 
of negative ions according to reactions different 
from those requiring lower energies. 

H-.—Taking E(H)=0.7 volt," we see from 
Fig. 5 that the following are the four possible 
processes requiring least energy according to 
which H~ might be formed by dissociation and 
capture : 


— 3.8 volts; 


(a) CH,—-CH;+H- 

(b) CH,-CH+H:2+H- —6.4 volts; 
(c) CH,—-CH:+H+H- —7.4 volts; 
(d) CH,-C+H+H:+H~-—9.9 volts. 


The observation that H~- ions are formed in 
three different processes all requiring less energy 
than (d) suggests rather strongly that (a), (b) 
and (c), respectively are the reactions occurring 
at the observed minimum energies 6.10.3 volts, 
6.9+0.8 volts, and 7.3+1.0 volts. This interpre- 
tation requires that the total kinetic and excita- 
tional energies of the products formed be at 
least 6.1—3.8=2.3 volts in reaction (a), 6.9—6.4 
=0.5 volt in reaction (b), and 7.3—7.40 volts 
in reaction (c). 

If we take 7(CH) = 11.7 voltsand J(CH2) = 12.0 
volts Fig. 5 indicates that the 18.0-volt appear- 
ance potential of H~ can only be the minimum 
energy for one of the reactions CH,j-H-~+CH;* 
or CH,;—-H-+CH*+Hs2:. The energy required 
for the first of these would be 14.4 volts (that 
needed for the reaction CH;—H+CH;*) — E(H) 
or 13.7 volts plus the kinetic and excitational 
energy of H~+CH;*. Hence the latter energy 
would be 4.3 volts. The energy required for the 
second reaction, aside from the kinetic and 
excitational energy of H~+CH*+He, would be 
7.1 volts+J(CH)—#(H) or 18.1 volts. Thus if 
H~- ions are produced according to this process 
the products can be formed in their normal 
states with practically no mutual kinetic energy. 
Because of the relatively strong CH;* ionization 
and unresolved ion current at m/e=13 produced 
at minimum energies than 18.0 volts, 
neither the CH;* ions, if formed in the first 
reaction, nor the CH? ions, if formed in the 
second reaction, at 18.0 volts minimum energy 
could have been detected. 


= Hylleraas, Zeits. f. Physik 65, 208 (1930). 


less 
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C~.—The striking feature of the curveshownin 
Fig. 4 for C~ ions is the absence of a detectable 
peak at low electron energies. From this it 
appears that a C~ ion is only formed along with 
a positive ion. This ion may be H*, H2*, or Hs*. 
The assumption that the reaction CH,—-C- 
+H+H;* occurs at a minimum energy of 
27.4+0.6 volts seems distinctly unreasonable on 
the grounds that the additional H;* ionization 
would have been detected in the appearance 
potential curve for that ion. It is quite possible 
that 27.4+0.6 volts may be identified with the 
27.9+0.5 volt appearance potential of H,* ions. 
The fact that at an electron energy of 50 volts 
H;* ionization is about 80 times as intense as C~ 
ionization shows that most of the H»* ions are 
formed according to some reaction (a) other 
than (b) CH,—-C-+2H+H,;*. From the dis- 
cussion of A(H,*) it is evident that probably (a) 
is the reaction CH,-C+2H+H,*. The mini- 
mum energy required for (b) should, except for 
differences in the kinetic and _ excitational 
energies of the products, be less than that 
required for (a) by E(C) =1.9+0.7 volts.2® How- 
ever, in the appearance potential curve for H2* 
there is no abrupt change of curvature at 
electron energies above 27.9 volts showing that 
if both (a) and (b) take place the minimum 
energies required are about the same and hence 
that C-+2H+H,* must possess about 1.9 volts 
more energy when formed than C+2H+H¢?. 

A simpler and more likely explanation of the 
formation of C~ ions is provided by the assump- 
tion that they are accompanied by Ht? ions. 
The intensity of H+ ionization is ~10* times 
that of C~ ionization at 50 volts electron energy. 
Hence we should not expect to detect the appear- 
ance of the relatively small H* ionization pro- 
duced in the process CH,—-C-+2H+H? at 27.4 
volts minimum energy. Without allowance for 
excitational and mutual kinetic energies of the 
products the minimum energy for this process 
should be D(CH,)+J(H)—£(C) =26.7+0.7 


*? The author estimated the value E(C) =1.9+0.7 volts 
by plotting the square roots of the ionization potentials of 
N, O*, and F** vs, their atomic numbers and extrapolating 
to atomic number six. Glockler (Phys. Rev. 45, 111 
(1934)), by a similar extrapolation, has obtained the value 
1.37 volts for E(C). Use of this value instead of 1.9 volts 
lowers the total excitational and kinetic energy of products 
as well as the upper limits (6) for D(CH,) and L(C) by 
0.5 volts or 12 kcal. 


No 
~ 
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volts. Hence on this assumption we may conclude 
that the products possess at least 27.4—26.7 
=0.7+0.9 volts of energy, or that D(CH,) 
=A(C-)+E(C)—J(H) and hence that” 


| D(CH,) =15.8+0.9 volts or 364+21 kcal., | (6) 
IL(C) = 6040.9 or 139421 “ .{\? 


The latter conclusions like (4) and (5) provides 
fairly strong evidence against the value 168.8 
kcal.'* for L(C). 

CH-.—The first appearance potential of CH- 
ions may readily be explained as the energy 
required for the process CH,—>CH~+3H. This 
yields the result that E(CH)=11.6—10.2=1.4 
volts according as the products are formed with 
or without energy. The other reaction (CH, 
—CH-+H-+H.,) which could account for these 
ions demands that the products be formed with 
an energy of 2.1 volts+H(CH). The second 
appearance potential of CH™ ions could be the 
minimum energy for the reaction CH,—-CH- 
+H;*. This interpretation is indicated by the 
agreement of A(CH~)=24.5+0.5 volts with 
A(H;*)=25.341.0 and by the fact that the 
intensities of H;+ and CH™ ionizations are 
approximately equal at all electron energies at 
least as high as 60 volts. Again using D(H;*) to 
represent the energy of the reaction H;+—H,z 
+H*, we see from Fig. 5 that only 20.6 volts 
— E(CH)—D(H;*) should be required to break 
a CH, molecule apart into a CH™~ ion and an 
H;* ion in their normal states without mutual 
kinetic energy. Hence if this be the correct 
interpretation the amount of the impacting 
electron’s energy which must go into kinetic and 
excitational energy of the two ions is 3.9 volts 
+E(CH)+D(H;*). 

The process whereby a CH™- ion and other 
products of a CH, molecule could be formed with 
the least total kinetic and excitational energy is 
CH,—CH-+2H+H?. Here we see from Fig. 5 
that this energy would be (24.5—25.1) volts 
+E(CH) and is =0.8+0.6 volts if we take 
E(CH)=1.4 volts as indicated in the previous 
paragraph. 

CH,;-.—When a CH:~ ion is produced by 
dissociation of a CH, molecule and capture of the 
impacting electron, the other atoms of the 
original molecule may dissociate as 2H or as Hs. 
The energy required to excite or give kinetic 
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PABLE IV. Values calcuiated for D(CH,) and L(C). 
BASIS OF CALCULATION D(CHs)*< L(C)-= 
A(H*) =29.4+0.6 volts 15.9+0.6 volts or 366+14 kcal. 6.1+0.6 volts or 141+14 kcal. 
A(C*) =26.7+0.7 volts 15.5+0.7 volts or 357+16 kcal. 5.7+0.7 volts or 131+16 kcal. 
A(C~) =27.4+0.6 volts 15.8+0.9 volts or 364+21 kcal. 6.0+0.9 volts or 139+21 kcal. 


energy to the products formed would be 0.8 
volts+E(CHez) in the first case or 5.2 volts 
+HE(CHg2) in the second case. Since the latter 
value seems high it is probable that the reaction 
occurring is CH,—CH,~+2H. The formation of 
a CH>- ion at energies greater than 23.4+0.6 
volts must proceed according to one of the two 
reactions CH,->CH.-+H+H?* or CH,y—-CH;- 
+H,+t. Aside from the excitational and mutual 
kinetic energies of the products, the minimum 
energies for these two reactions are respectively : 
21.6 volts—E(CH2); and 19.0 volts—E(CH2). 
Thus the total kinetic and excitational energies 
of the products must be: 1.8 volts+E(CHg); or 
4.4 volts+£(CH,). It is unlikely that the positive 
ions formed in one or the other of these processes 
would be sufficiently numerous to be distinguish- 
able from those formed in other ways. 


SUMMARY AND CONCLUSION 


A summary of the conclusions which were 
drawn in the above discussions from the values of 
dissociation energies given in Fig. 5 will be found 
in the last three columns of Table I. The values 
in column five are calculated by adding ionization 
energies to these dissociation energies. They 
represent minimum energies for the correspond- 
ing reactions exclusive of the total kinetic or 
excitational energy products must have when 
formed by electron impact. Values of the latter 
quantities are given in column six where possible 
and are obtained by subtraction of the values of 
column five from the corresponding appearance 
potentials in column three. 

It should be emphasized that the quantities of 
Fig. 5 are not certain and that hence some of the 
interpretations may be incorrect. One cannot be 
completely sure at present of any reaction except 
(12) and (13). However, it is a fact of considerable 
significance that a reasonable explanation for 
every observed appearance potential exists. At- 
tention should be drawn specifically to the fact 
that in all cases except A(CH,*)=22.9 volts, 


A(H~-) =6.1 volts, and probably the appearance 
potentials of H;* and CHe2°-, interpretations are 
possible which require that the total excitational 
and kinetic energy of products formed by electron 
impact from CH, molecules be less than two 
volts. 

It will be noted that the ions H.*, H;*, C+, C-, 
CH-, and CH;- are not produced unless in each 
case the energy of the electrons is at least slightly 
greater than the largest of the calculated mini- 
mum energies for possible processes. This may be 
true of CH* ions as well. The higher of the two 
appearance potentials of H* and the higher ap- 
pearance potential of CH,* also satisfy this con- 
dition. The only type of ion for which an appear- 
ance potential satisfying this condition was not 
detected is H~. These facts seem to indicate that 
at least in most of the cases mentioned, the 
possible processes which require the largest calcu- 
lated minimum energies are the ones which 
actually take place: i.e., that most ions are only 
formed when the parent molecules are broken 
into as many parts as possible. This conclusion 
gives one some reason to prefer the reactions (5), 
(8), (21), and (24) to the reactions (6), (9), (22), 
and (25) (column four of Table I). 

This conclusion also adds weight to the rela- 
tions (4), (5), and (6) which are gathered together 
in Table IV. These relations are based on the 
assumptions that the appearance potentials in 
column one of Table IV are minimum energies for 
the reactions (3), (7), and (19) in column four of 
Table I. The values quoted are in each case less 
than or equal to the true values of D(CH,) and 
L(C) according as the total excitational and 
kinetic energies of the products is greater than or 
equal to zero. Thus as a further conclusion, we 
may say that of the three values of L(C) de- 
termined from a predissociation limit in CO", the 
results for methane provide good evidence against 
the value 168.8 kcal. and are completely con- 
sistent with the value 123.6 kcal. They are also 
consistent with the lower of the two possible 
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values (136 kcal.) obtained by Lozier”’ from data 
on the products formed by electron impact in 
CO. 
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Thermionic Emission into Dielectric Liquids 


Epwarp B. BAKER AND Howarp A. BoLtz 
The Detroit Edison Company, Detroit, Michigan 


(Received September 28, 1936) 


The presence of thermionic emission at ordinary temperatures for cathodes in contact with 
dielectric liquids is demonstrated. The current voltage relation is shown to obey the law 
T=Ie(*¢?D) kT), where J is the current, F the field, and D the dielectric constant of the liquid. 
This is a modification of the Schottky law in vacuum. The effects of adsorbed hydrogen and 
oxygen, of surface irregularities, and of work function are discussed. The formation of a wax by 
electron bombardment of the liquid is also reported. 


A. INTRODUCTION 


ANY theories have been advanced to ex- 

plain the mechanism of current conduction 
in pure dielectric liquids in the region of high field 
strengths. A good resume of these may be found 
in Nikuradse’s book.' The most generally ac- 
cepted of these is the process of ionization by 
collision, in analogy to gaseous conduction. It is 
assumed that ions are initially present in the 
liquid which are the current carriers at low field 
strengths. These are supposed to be formed by 
some ionizing agent such as x-rays or cosmic 
rays, or from ordinary electrolytic dissociation. 
At intermediate field strengths this ion current 
saturates and becomes independent of field until 
such a value is reached that the ions (or electrons) 
attain sufficient velocity to form new ions by 
collision. New ions are formed at such a rate that 
the current becomes an exponential function of 
the voltage J=J e*” which is found to agree 
fairly well with experiment, considering J» and C 
as adjustable constants. Most published data 
available to the authors are not sufficiently 
accurate, however, to distinguish this relation 
from J =J,e*’", for example 


In contrast to the above view, the results 





' A. Nikuradse, Das fliissige Dielectrikum, pp. 135-136. 


obtained in our investigation prove that, in the 
case of the liquid investigated (pure toluene), at 
least, the only important of current 
conduction is by electrons emitted from the 


means 


cathode thermionically. 

This mechanism was first proposed in tangible 
form by H. Edler and O. Zeier? who found the 
temperature variation of the current to corre- 
spond to the Richardson equation for thermionic 
emission. Although they did not derive the 
current voltage relation given in this paper they 
found the current to vary less rapidly with 
voltage than J = J,e*’. 

L. Inge and A. Walther’ also disagree with the 
mechanism of ionization by collision and invoke 
instead the cold cathode effect, which would 
require a current voltage relation of the form 
I=I,V*e-*'". They plot log I/ V? 
against 1/V and obtain not a straight line, but a 
curve which they believe to be approaching a 
straight line in its high voltage portion. However, 
if their data are replotted as log J against »/ V, a 
straight line is obtained in verification of the 
formula J = J,e*’" of the present paper. 


essentially 


2H. Edler and O. Zeier, Zeits. f. Physik 84, 356 (1933). 
31. Inge and A, Walther, Tech. Phys. U.S. S. R. 1, 539 
(1935). 
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B. THEORY 

I. In vacuum 

The emission of electrons from a hot cathode in 
a vacuum is governed by the well-known 
Richardson law, J[=AT“%e~**/*? where J is the 
current density, ¢ is the work function in electron 
volts and eg =x is the amount of energy required 
to remove an electron from the metal to a 
vacuum. If a sufficient extracting field is applied 
to the cathode the current becomes saturated 
except for the Schottky correction, which intro- 
duces a dependence upon the field F of the form 


[=A T*e-* kT ge(eF)! kT — Toe’? } kT 


This Schottky relation, suitably altered to in- 
clude the effect of a medium possessing a 
dielectric constant greater than unity, provides 
the limiting form at high fields of the current- 
voltage relation for electronic conduction in a 
dielectric liquid. Before introducing this modifi- 
cation, it will be convenient to make a few 
remarks about the Schottky law in vacuum which 
will also be pertinent when the dielectric liquid is 
present. 
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Becker and Mueller‘ have shown that the 
electrons emitted from the cathode surface must 
have a normal component of velocity sufficient to 
project them beyond the attracting influence of 
the surface fields, into the region where these 
fields are smaller than the external extracting 
field. They have shown that these fields may be 
divided into three classes according to their 
range of action. First, are the powerful but 
extremely short range fields of the metal atoms 
that are responsible for the major portion of the 
work function. Second, are the image forces 
which predominate out to about 10-° cm. At 
still greater distances, forces due to nonuni- 
formities of the work function over the surface 
are most important. 

At the actual surface of the metal the atomic 
fields are of the order of 108 volts/cm which is the 
magnitude required for the cold cathode effect. 
Moving out from the cathode, the intrinsic field 
falls off rapidly and it is possible to neutralize it, 
or to produce a surface of zero field strength, by 
applying an external field. The distance of this 
surface of zero field from the cathode is dependent 
upon the strength of the applied field. On the 
cathode side of this surface, the net field is in the 
sense that tends to return electrons to the 
cathode, and on the other to send them to the 
anode. An electron consequently must have 
sufficient initial velocity to reach this surface 
against the adverse field in order to be emitted. 

In Fig. 1, the surface fields are plotted against 
distance on a logarithmic scale. In the range in 
which the image field predominates, the curve 
becomes a straight line corresponding to the 
function F;=«/42?, where z is the distance of the 
electron from the cathode, and e« is its charge. 
Above several hundred angstroms the patch 
fields are most important and below about 10A 
the atomic fields of the metal atoms assume 
control. Each end of the curve is shown to 
broaden out, representing the possible variations 
in value of the field at various points over a plane 
surface parallel to the plane cathode, due to 
patchy work function, and in the case of the 
short range fields, to surface irregularities as well. 

A cross section showing the zero field surfaces 





4J. A. Becker and D. W. Mueller, Phys. Rev. 31, 431 
(1928). 








for various applied fields is pictured above the 
graph. If the metal surface, A, is not perfectly 
plane and of uniform work function, only those 
zero field surfaces, such as D lying within the 
effective range of the image force will be plane. 

The geometrical surface irregularities that 
exist even on a highly polished electrode cause 
deviations from the Schottky law at high applied 
fields. Points cause most of the emission from a 
surface of uniform work function, if the applied 
field is high enough to cause the surface of zero 
field (C) to be appreciably distorted by their 
presence. Thus, the more nearly flat the cathode 
is, the higher must be the applied field, or the 
closer to it must be the surface of zero field, in 
order to have the effect of points of any conse- 
quence. In investigating the cold cathode effect, 
it must be flat even on an atomic scale in order to 
avoid the influence of points. 

If the applied field is F, and the surface is 
perfectly plane, we have the ordinary Schottky 
relation. If, however, a point exists on the 
surface, the applied field at the tip of the point is 
not F but poF, where fp is a field factor (>1) 
depending on the sharpness of the point. At the 
surface of zero net field, the applied field is pF, 
where p is a function of F, equal to pp at the 
point and becoming unity at a sufficient distance 
from the cathode. The Schottky relation for the 
current from the point is 


T= [geen P'ik7, 
If the point is sufficiently sharp, i.e., if p is large, 
the contribution of current from the point may 
overshadow that for all the remainder of the 
surface. This modified relation may be approxi- 
mately obeyed if a number of points of about the 
same size exist simultaneously. 

Since p is a function of F, a straight line will 
not be obtained if log J is plotted against 4/F as 
is commonly done to illustrate the Schottky 
effect. Instead, a curved line having increasingly 
greater slope with increasing field will result. 
However, if p is a slowly varying function of F, 
the curvature is slight and the curve may appear 
straight over the portion investigated, but will be 
too steep. Since the knee of the curve log I vs. 4/ F 
obtained experimentally (caused by the approach 
to the limiting Schottky line) has the opposite 
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curvature, a point of inflection is produced and 
the appearance of straightness enhanced. 

The Schottky law is derived on the assumption 
that the image law force F;=«/4s* is the only 
long range force of importance. This is ac- 


complished® by showing that an external field, F, 
produces a change Ag=(eF)! in the effective 
work function by neutralization of the image 
force at the critical distance z.= (€/4F)! from the 
cathode surface. At lower applied fields the effect 
of the patch fields overshadows the image law 
and the Schottky law does not apply. 


II. Modification introduced by the presence of a 
dielectric liquid 
If we replace the vacuum by a dielectric liquid 
we may expect the phenomenon to be altered by 
the following considerations: 


1. The work function will be lowered because the 
amount of energy necessary to transfer an elec- 
tron from the metal to a vacuum is greater than 
that required to transfer it to a dielectric liquid. 
This difference in energy corresponds to the 
energy necessary to remove a “‘free’’ electron 
from the liquid to a vacuum. By a free electron in 
a dielectric liquid we mean an electron not in 
“chemical’’ combination with molecules of the 
liquid but simply immersed in it, the polarizing 
forces its field exerts on the surrounding mole- 
cules of the liquid being responsible for the 
difference of energy in question, which we may 
term the “heat of immersion of an electron.’’ We 
may expect that this amount of energy will be 
independent of the nature of the metal, and as a 
result that the relative order of the metals in the 
photoelectric series will be the same as in vacuum. 
We may also predict that the amount of the 
lowering of the work function will be at least 
roughly dependent upon the dielectric constant 
of the liquid, being greatest for liquids of highest 
dielectric constant. 

2. The temperature used must be lower since the 
dielectric liquid would boil or decompose at the 
temperatures used in studies of thermionic emis- 
sion in vacuum. Since 7 and ¢ occur in reciprocal 
relation in the exponent of Richardson’s formula, 
we may expect that a lowering of both T and ¢ 
simultaneously will permit appreciable emission 


5]. Becker, Rev. Mod. Phys. 7, 108 (1935). 
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at room temperature, if ¢ is lowered to such an 
extent that the ratio is not greatly altered. 

3. The saturation of the current will be much less 
complete at high fields partly because the temper- 
ature is lower; the exponent in the Schottky 
correction is consequently larger, and the de- 
pendence on the field of more importance. In 
addition, the Schottky correction must be modi- 
fied by the fact that the image force is weakened 
by the presence of a liquid of dielectric constant 
ereater than one. The image force in a dielectric 
of dielectric constant D is F =e 42°D, the dis- 
tance from the cathode of the surface of zero field 
is t.-=(e/4FD)', and the revised Schottky cor- 
rection is 

I 1.91(FD)! 


log io = 


T=IpetsF2)/'T or 
I T 


with Fin volts per cm.* This is the form that the 
current voltage characteristic takes at high fields. 
The slope of the line logig J=c,/F+log Jo is 
much greater than in a vacuum. 

As a rule, somewhat higher fields are required 
to obtain the straight line portion of the curve 
than is the case in vacuum. For a given field 
strength, the surface of zero field as determined 
by Z.=(e/4FD)! lies closer to the cathode, and 
the effect of points on the surface is much more 
pronounced. Stated in other words, in order to 
avoid the effect of points, a much smoother 
surface must be employed as cathode in a 
dielectric liquid, than is necessary in vacuum. 

4. Collisions of electrons with, and attachment to 
molecules of the liquid may interfere with the 
direct application of the theory of emission in 
vacuum to that in a liquid. Evidently, however, 
these are of minor importance for toluene since 
the data agree with the modified Schottky 
equation. At low fields collisions and attachment 
to molecules may be of importance but it would 
be difficult to distinguish their effects from effects 
caused by patchy fields. 

The authors have observed evidences of chem- 
ical action caused by electron bombardment. If 
high field strengths are used and electrode con- 
ditions are such that high currents are flowing, 
enough of a wax-like substance is deposited on 

® One may obtain this expression by introducing the new 


value of the image force in Becker’s derivation of the 
Schottky correction (reference 5). 
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Fic. 2. The electrodes and quartz support. 


the anode to be readily discernible in about a 
week’s time. This substance is similar in ap- 
pearance, and in its complete insolubility in 
organic liquids, to the substance called X which 
is often observed in underground power cables 
having oil and paper as a dielectric. Its presence 
there has always been ascribed to electron 
bombardment, resulting from discharges in gas 
pockets in the insulation. This work, however, 
shows that it may also be produced in the liquid 
phase alone by electrons produced thermionically, 
and therefore that its presence does not presup- 
pose the existence of gas pockets, or voids. 

It is conceivable that the molecules of some 
liquids might have sufficient electron affinity 
seriously to impair or completely to suppress the 
electron emission altogether by attaching the 
electrons before they reach the surface of zero 
field. In such an event the ions thus formed 
would be discharged again at the cathode. 


C. APPARATUS 


The electrodes used were platinum-iridium and 
nickel, turned to the form shown in Fig. 2. The 
faces were ground flat and polished on a metal- 
lographic wheel. The electrodes were then tied 
with wire of the same metal to the small quartz 
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frame shown, being aligned by means of grooves 
ground in the quartz. The bottom electrode 
shank was drilled so that the assembly could rest 
on a tungsten point serving as a high voltage lead 
in. The upper electrode was drilled perpendicular 
to its axis near the end to receive a stiff platinum 
wire serving as a measuring lead and as a means 
of carrying the assembly. 

The quartz ring shown was platinized. It made 
contact with a platinized ring in the apparatus 
proper, and served as a leakage guard. A guard to 
preserve a uniform field is neither necessary nor 
desirable. 

The electrode distance was determined by 
means of a traveling microscope after the 
assembly was in place. 

The apparatus proper is shown in Fig. 3. It 
consists of a boiler containing toluene and 
metallic sodium, a central mixing chamber and 
the electrode chamber. It was designed to permit 
ready accessibility to the electrodes and the 
addition of soluble impurities to the toluene 
without allowing the entrance of air or moisture. 
By continued distillation from sodium the 
toluene is constantly repurified and desiccated, 
and the electrode chamber is continually washed 
out with fresh toluene. 

Since the present work is concerned only with 
the demonstration of the presence of thermionic 
emission and the agreement with the modified 
Schottky law, the absolute value of the work 
function itself is of importance only insofar as it 
effects the general order of magnitude of the 
currents obtained. Consequently, it is not neces- 
sary to work in high vacuum, which would be 
extremely difficult in the presence of an organic 
liquid. A nitrogen atmosphere was maintained, 
however, partly to avoid oxidation of the toluene 
during distillation, and partly because adsorbed 
nitrogen is known’ to have a comparatively 
small effect on the work function of a cathode in 
vacuum. The effect of other adsorbed gases will 
be treated in connection with the experimental 
results. 

The toluene vapor leaving the boiler passed 
over metallic sodium held in a bend in the tube 
leading from the boiler to further aid in removing 
moisture, part of which may distill before it can 


7J. H. de Boer, Electron Emission and Adsorption 
Phenomena (Cambridge Press, 1935), p. 147. 





Fic. 3. The electrode chamber and still. 


react with the sodium immersed in the toluene. 
The vapor condensed in the two water con- 
densers, which served to protect from the vapor 
the ground glass joints, by means of which 
the electrodes and foreign substances were 
introduced. 

A stopcock attached to a cross arm connecting 
the two entrance tubes led toa P.O dry ing tube, 
a pyrogallol bulb and a rubber bag reservoir. A 
stream of nitrogen flowed through this train 
whenever a joint was opened, preventing the 
influx of air and moisture. 

A siphon connected the mixing reservoir with 
the electrode chamber permitting either a con- 
tinuous flow of toluene from one to the other or 
an accumulation of toluene in the mixing cham- 
ber, which may later be run through the electrode 
chamber. This is accomplished by tilting the 
apparatus, which introduces a few gas bubbles 
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Fic. 4. Current-voltage characteristics for A oxygen on 
platinum-iridium, and B, C, D hydrogen on platinum- 
iridium. Temperature was 45°C. 


in the siphon, stopping its action until the level 
in the mixing chamber is high enough to start the 
siphon again. Foreign materials may be added to 
the accumulated toluene, which is then siphoned 
to the electrode chamber. 

In some cases the electrodes were heated in 
vacuum by means of radiofrequency induction, 
cooled and then quickly transferred to the 
apparatus, exposing them to the air only for a 
short time. In this way the nature and amount of 
adsorbed gases could be controlled sufficiently 
well for the present purpose. 

The electrical equipment has been described 
previously. 

The toluene used was C.P. grade further 
treated with mercury, concentrated sulphuric 
acid, sodium hydroxide, washed with water, and 
distilled from metallic sodium before introduction 
into the apparatus. 


D. EXPERIMENTAL RESULTS 
The electrodes first used were made of plati- 
num-iridium, because they were easy to clean by 
washing with chromic acid, then with distilled 


Sci. Inst. 7, 50 


SE. B. Baker and H. A. Boltz, Rev. 
(1936). 
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water and finally heating in a gas flame. Ex- 
tremely low conductivity was observed even at 
high field strengths the first few times they were 
treated in this manner. Thereafter the conduc- 
tivity became progressively higher with each 
treatment although there was no reason to 
suspect any change in the composition of the 
toluene. It was suspected that the electrodes had 
absorbed hydrogen from the flame, and heating 
them in a vacuum confirmed this. 

Evidently an adsorbed layer of hydrogen 
lowers the work function in the same manner as 
in vacuum.’ Now it is well known that it is 
extremely difficult completely to outgas platinum 
by heating in vacuum, but that it is easily 
possible to replace an adsorbed layer of hydrogen 
by one of oxygen, by heating platinum in oxygen, 
and conversely oxygen may be displaced by 
heating in hydrogen. Nitrogen, however, has 
little tendency to be adsorbed, or if adsorbed, 
does not appreciably affect the work function of 
platinum. 

These facts offer a method of control of the 
condition of the surface, even if the electrodes be 
transported through air after the desired con- 
dition is attained. It was found possible to 
produce low or high conductivity at will by 
having a layer of oxygen or of hydrogen on the 
surface. Fig. 4 shows current-voltage charac- 
teristics for platinum-iridium electrodes plotted 
as log J vs. 4/ F. Curve A is for adsorbed oxygen 
and curves B, C and D for adsorbed hydrogen. 
Comparing curves A and D, we observe a ratio 
of the currents of the order of 10°, an enormous 
difference, which incidentally proves that all 
but an entirely negligible portion of the current 
represented by curve D is due to electrode 
condition, since the purity of the toluene was the 
same for A and D. This eliminates the possibility 
of the current being due to ions generated in the 
liquid by ordinary electrolytic dissociation, even 
at fields as low as 10 kv/cm. 

Adsorbed layers of hydrogen and oxygen on 
platinum are present,'® at least for the most part, 
in ionic form. In the case of hydrogen the ions are 
positive. This positively charged layer lowers the 
effective work function by tending to attract 
electrons out of the metal surface. Together with 


9 Cf. de Boer, reference 7, pp. 149-151. 
10 de Boer, reference 7, p. 150. 
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their image charges in the metal these ions form 
dipoles which on the cathode are oriented in the 
same direction as the applied field. In addition 
they are to a certain extent mobile on the 
surface," even at room temperature. On a per- 
fectly flat cathode, i.e., where the applied field is 
strictly constant, there will be no tendency for 
these dipoles to be displaced by the field. If, 
however, hills and valleys exist, they will be 
driven to the valleys, toward the regions of 
lowest field. The oxygen dipoles, since their 
orientation cannot be greatly influenced by the 
applied field, will be driven to the points. When 
the applied field is removed, they will again 
distribute themselves more uniformly over the 
surface by virtue of their mutal repulsion and by 
thermal diffusion. 

Since an adsorbed layer of hydrogen lowers the 
work function, and since the hydrogen ions are 
driven to the valleys on the cathode, where the 
applied field is Jess than that calculated from the 
applied voltage and the distance with the as- 
sumption that the cathode is plane, we may 
expect that the slope of the Schottky line will be 
less than for a plane surface, i.e., »<1. That this 
is indeed the case may be observed in curves 
B, Cand D of Fig. 4. 

In the case of adsorbed oxygen the ions are 
driven to the points, but since oxygen increases 
the work function, the points are thereby 
rendered relatively inactive and the remainder 
of the surface being partially cleared of oxygen 
by the field, contributes most of the current, by 
virtue of its lower work function. This fact 
enables us to obtain a closer agreement with the 
theoretical slope, since the remainder of the 
surface, devoid of the points, is practically flat. 
This could not be done in the case of platinum 
because the approach to the Schottky line lies at 
such high field strengths (curve A) that break- 
down is too imminent. In fact, while the current 
from oxygen on platinum is very low, it is 
erratic due to the sporadic removal of oxygen by 
the applied field, which suddenly exposes points 
of relatively low work function. Thus while an 
adsorbed layer of oxygen serves as a protection 
from large electron currents, its removal by the 
field may perhaps be a cause of breakdown for 


1 de Boer, reference 7; also E. K. Rideal, An Introduction 
to Surface Chemistry (Cambridge Press, 1930). 
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Fic. 5. Current-voltage characteristics for A and C hydro- 
gen on nickel and B nickel. Temperature was 45°C. 


electrodes in this condition, and electrodes 
usually are in this condition unless specially 
treated. 

It may be emphasized at this juncture that the 
effective emitting area of the cathode was not 
entirely under control in these experiments. A 
displacement of a Schottky line parallel to itself 
can occur either as a result of a change in 
emitting area or a change of work function. The 
slope of the line, however, is independent of 
these factors, and depends only on the dielectric 
constant, the temperature and absolute con- 
stants (except if #1). Thus, while the fact that 
the curves exhibit a straight line portion plotted 
in this manner is in itself good proof of the theory, 
the additional fact that the slope of this line has 
exactly the correct value, independent of the 
distance d between the electrodes, cannot pos- 
sibly be attributed to chance agreement. Other 
attributes of the current, such as its instan- 
taneous response to voltage changes (except with 
large amounts of adsorbed gases, which diffuse 
somewhat slowly at room temperature) are also 
in accord with the theory. 

The work function of platinum is about 6 
volts, that of nickel about 5 volts. We may 
therefore expect greater emission from nickel 
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nickel. Temperature was 45°C. 


than from platinum. In addition, we may expect 
the work function of nickel to be less influenced 
by adsorbed hydrogen than that of platinum, as 
is the case in vacuum, since the adsorbed ions are 
less firmly bound on nickel than on platinum 
because of its lower work function.” 

In Fig. 5, curve A was obtained with a nickel 
electrode that had been heated in hydrogen. It 
shows the smaller slope characteristic of emission 
from valleys. Curve C was obtained by partially 
outgassing the electrodes and still shows the 
effect of hydrogen. Curve B was obtained by 
extreme outgassing over a period of 20 hours. 
Even though transferred to the apparatus 
through the air the nickel was partially bare and 
the curve exhibits the steeper slope due to the 
presence of points not covered by oxygen. 

2 de Boer, reference 7, p. 151; R. W. Gurney, Phys. Rev. 
47, 479 (1935). 
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In Fig. 6 curves A, B and C are examples of 
nickel outgassed to remove hydrogen, and then 
soaked in oxygen at room temperature. They all 
show the correct Schottky slope. 

Curve D in Fig. 6 was obtained by adding 
iodine to the toluene. In this case the effect of 
adding iodine may be seen by direct comparison 
with curve B, which represents the condition 
just previous to the addition. Since iodine is an 
electronegative element, it should behave in a 
manner analogous to oxygen and therefore raise 
the work function, and decrease the current. 
Instead, the current was increased considerably. 
This increase may possibly be attributed to the 
electrolytic deposition of hydrogen on the cath- 
ode supplied by the hydroiodic acid liberated by 
the action of iodine on toluene, which takes place 
very slowly at ordinary temperatures. The fact 
that the slope of the high field portion of the 
curve is less than the theoretical supports this 
view. The current represented by curve D, once 
established, was independent of the iodine con- 
centration, being unchanged by repeated washing 
out of the electrode chamber with fresh toluene. 
It had all of the properties of an electron current. 
This curve was included particularly to illustrate 
the large electron currents that are possible even 
at low fields. 

Since we have shown that these electron 
currents not only impair the insulating qualities 
of dielectric liquids but can cause chemical 
deterioration of them as well, the possibility of 
increasing the work function of electrodes such 
as the conductor in oil insulated cable, becomes 
evident. It is curious, in this connection, that 
the presence of oxygen should appear to be 
beneficial, since it has long been considered to be 
harmful because of the danger of oxidation of the 
organic insulating liquid. 

In conclusion the authors wish to express their 
appreciation of the efforts of the personnel of the 
research laboratory whose cooperation made this 


work possible. 
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In the present note the question of the separation of the 
kinetic energy of the center of gravity from the total kinetic 
energy of the nucleus is considered. This separation is 
necessitated in the calculation of nuclear energy levels by 


the use of the Hartree method. It is found that under 


certain conditions the wave function for the whole nucleus 
in absolute coordinates can be replaced by a product of the 
same function of the coordinates of the particles relative 


N the calculation of the energy levels of light 

nuclei from the individual particle model 
(Hartree approximation), it is the usual pro- 
cedure to use wave functions in which the nuclear 
particles are coupled to a fixed point (absolute 
coordinates) rather than to each other. The 
kinetic energy T of the nucleus calculated with 
such wave functions will contain the spurious 
kinetic energy of the center of gravity which 
must be deducted. For this purpose several 
authors have replaced the mass M of a single 
particle by the reduced mass ./(1—1/N) where 
N is the total number of particles in the nucleus.! 
However, such a procedure is correct only for a 
two particle problem or for the case that all the 
particles are in the 1s shell. It is our purpose in 
this note to investigate this point in a somewhat 
more general case. This investigation, we shall 
find, is greatly facilitated by the fact that the 
wave functions ordinarily used, i.e., those in an 
oscillator potential, actually contain the center 
of gravity only in a trivial factor. By separating 
off this factor one obtains directly wave functions 
containing only relative coordinates. 

For the oscillator potential 


U=—Ur+} Mer’, (1) 
w being the frequency of a classical oscillator of 
mass M (proton or neutron mass) in the field U, 


the individual particle wave functions in abso- 
lute coordinates are: 


y=ce-™, (2) 


(2’) 


1s level 


2p level y=c'(x, y or z)e7*”"’, 


1Cf., e.g., S. Fliigge, Zeits. f. Physik 96, 459 (1936). 


to the center of gravity and a 1s wave function for the 
center of gravity. These conditions are that oscillator wave 
functions be used for the individual particles and that the 
nucleus contain only one partly filled shell in both protons 
and neutrons with all shells of lower energy 


number filled. It then follows that the kinetic energy of the 


quantum 


nucleus is equal to the value calculated in absolute coordi- 
nates minus the energy of a single 1s particle 


and in general for ‘‘energy quantum number’” &, 


Here ¢ and c’ are normalization constants, P;_, is 
a polynomial of degree k—1 and 


a= Mw 2h. (3) 


For the wave function WV of the whole nucleus 
we may take a product of two antisymmetric 
(determinantal) wave functions, the 
and one for the 
terminant being constructed from the spin and 
space wave functions (2) of the individual parti- 
cles. In all cases the wave function VY will contain 


one for 


neutrons protons, each de- 


oN. 
the factor e~**1" 
transformed into center of gravity coordinates R 


This factor may easily be 


and relative coordinates 9; with respect to that 
point: 
- 


R=Sr,/N (4) 


1 
and oi(é:, ni, £3) =1i(xs, vi, 23) —R(X, Y,Z), (5) 


writing the Cartesian components of the vectors 
in the brackets. The relative coordinates 9; are 
of course not independent but fulfill the relation 


N 
Loi=0 (6) 
1 


2 By the energy quantum number & we mean the quan- 
tum number which defines the energy of a single particle 
in the oscillator field. For the field (1) this energy is 
E=— Up+hw(k+ 3). The relation between & and the 
principal and azimuthal quantum numbers, » and /, is 
k=2n—l—1. See H. A. Bethe and R. F. Bacher, Rev. 
Mod. Phys. 8, 82 (1936), §32. 
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From (5) and (6) it follows that 
rV,.2 —aN R? —ar,.2 ‘es 
ee) rs -_ aNR e ad) pi : (7) 


We have now to consider the transformation of 
the determinants in the wave function which 
remain after the factorization of the gaussian 
term (7). For the sake of simplicity we shall first 
consider the case in which, for each kind of 
particle, the 1s shell is complete and the 2) 
shell is partly filled. In this case each determinant 


will have, in general, the form 


a By ooo Betis * V181 lle 
Qe Bors> VoBo eee 


Dir=|- - : (8) 


ae0ie *°° 


in which @ and 8 are the spin wave functions 
specifying the spin component in a given direc- 
tion sz. Now in D we may add to any column a 
constant multiple of any other column without 
changing the value of the determinant. Thus we 
may add to the x,a; column —X times the a; 
column and the former column becomes £;a; by 
(5). Similarly the y,8; column can be replaced by 
ni8; and so on without changing the value of the 
determinant. It follows therefore that 


D(r;) =D(9i) (9) 
and from (7) W(r;) =e~***(9,). (9’) 
It should be noted that in the derivation of (9) 
it is essential that the 1s shell be complete.* It is 
also true that the separation effected in (9’) 
depends on the properties of the oscillator wave 
functions in an essential way and this result will 
be valid only for such wave functions. 

The extension of this result to the case of 
higher shells is obvious. The relations (9) and 
(9’) will be valid in this case if only, for each 
kind of particle, all the shells of lower energy 
quantum number? be filled.‘ It is sufficient in 
each case to consider only a single monomial in 
the polynomial P; since the determinant with 
P, can be written as a sum of determinants in 


3 Of course the 1s shell for one kind of particle need not 
be complete if there are no 2 particles of this kind. The 
relation (9) is trivial in the case of only 1s particles. 

‘The number of filled neutron shells need not be the 
same as the number of filled proton shells. 


AND M. E. ROSE 


which P, is replaced by the various monomials. 
In the transformation of a column containing 
x*y*z? one gets £n*¢* plus additional monomials 
of lower degree. These additional monomials will 
coincide within a constant factor (powers of 
X, Y and Z) with elements in one of the columns 
belonging to a shell of smaller k and therefore 
give a vanishing contribution to the determinant. 

We may now proceed to the calculation of the 
nuclear kinetic energy. The transformation of the 
kinetic energy operator in absolute coordinates, 
N 
DA,,, gives, as is well known, the kinetic energy 
operator for the center of gravity, Ar/N, plus an 
operator Q» depending only on the relative 
coordinates. The operator Q%» is not the sum of 
the Laplace operators A», but will contain in 
addition cross derivatives of the form divg,- grado, 
because the o; are not independent (Eq. (6)). 
Therefore it is not simple to calculate the kinetic 
energy 7”, exclusive of the energy of the center of 
gravity, with this operator and the wave function 
V(o;). However, T’ can be easily calculated as 
the difference between the total kinetic energy in 
absolute coordinates, 7, and the kinetic energy 
of the center of gravity, Tr. T is the sum of the 
kinetic energies of the individual particles in the 
oscillator potential. For the kinetic energy of the 
center of gravity we have, using (3), 

nh? f'dRe-*X ®* Ape—@N F* 


T r= -——_ — 
2MN 


= 3hw; (10) 


JS dRe 2aN R? 


that is, just the kinetic energy of a single 1s 
particle. 

For the kinetic energy 7” of a nucleus with all 
shells (for both particles) of energy quantum 
number K or less filled and with g particles in the 
shell of quantum number K+1 we find 

K 


thw {g(2K+3)—3+ ¥ 2k(k+1)(2k+1)} 
k=1 


Tx’ 


=lhw{q(2K+3)-34+K(K+1)(K+2)} 


=1hw{(2K+3)N—3 
—}3K(K+1)(K+2)(K+3)}, (11) 


in which use has been made of the fact that the 
total number of particles of both kinds in the 
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shells with quantum number & is 2k(k+1). The 
case of greatest interest is that in which the 2 
shell is partly filled and the 1s shell complete. In 
this case K = 1 and the kinetic energy is 


hw hw 
T,’ =—(5q+9) =—(5N-—11). (11°) 
4 4 
The difference between the result for the 


kinetic energy obtained here and that resulting 
from the use of the reduced mass is 


ATx’=Tx N— hw, (12) 
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which for K =1 is 


AT’ =hw(N—4)/2N (12’) 
This correction varies between 25 and 50 percent 
of the energy ?/w for the nuclei of mass number 
between 6 and 16, and will have the effect of 
raising the energy levels as formerly computed 
by about 2 to 4 MV for these nuclei. 

It is with pleasure that one of us (M. E. Rose 
acknowledges his indebtedness to The American 
for a grant from the 


Philosophi« al Sc ciety 


Penrose Fund. 
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Spherical Symmetry of Self-Consistent Atomic Fields 
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HE method of the self-consistent field (s. c. f. 
for an N electron atom is equivalent to the 
assumption 


¥=Yill)---Yw(V) (1) 


in the variational principle. This leads (neglecting 
spin terms) to the equations 


Vivit(e,—2V)¥;=0 (j=1,---,N), (2a) 


" 
> f Wi(Px) | ?/re-dt,, (2b) 
k=1(k+)) 


where 
V (r;) = —Z/rj;+ 


Z being the nuclear charge, and atomic units 
being used. On the other hand, all calculations 
which have actually been made so far with the 
s. c. f. have proceeded on the assumption that the 
field for each electron is spherically symmetrical. 
This is effected by making the additional as- 
sumption 


v;(r;) = (24) ei Pm 51,(0)Rlr;) (3) 


in the variational principle, where we denote by 
Pni(@) the normalized tesseral harmonic P;'™ 
(cos 6), m; and 1; being the magnetic and azi- 
muthal quantum numbers of the jth electron. 


The assumption (3) is, however, incompatible 
(2b 


s electrons. Excluding this special case, therefore, 


with Eqs. (2a), unless all the electrons are 
(3) imposes a restriction on the form of the wave 
function which is additional to that implied by 
(1). Consequently, although the assumption of 
spherical symmetry is entirely reasonable from 
the physical point of view and simplifies con- 
siderably the numerical work, it would appear 
that the generalized Eqs. (2) are more accurate 
than those ordinarily used, and this seems to have 
been the generally held opinion. If this were really 
the case, one might solve the more general equa- 
tions by approximate methods and estimate the 
magnitude of the corrections that would thereby 
result. A closer examinaticn, however, leads to 
the conclusion that if the assumption (1) is made 
(or’ the generalizations considered below), the 
central field assumption is, in effect, a necessary 
one if the solution obtained is to fulfil certain 
general requirements. 

Let us first of all examine the nature of the 
solution of Eqs. (2). There is a degree of arbi- 
trariness about the solution, but we may safely 
assume that the central field approximation is a 
good one. Let us, then, see what type of solution 
results when we solve (2) by successive approxi- 
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mations, starting with a definite central field 
solution (3) (in which the functions R; are those 
obtained from the usual Hartree equations) as 
initial approximation. This may be accomplished 
by continually substituting the lower order ap- 
proximation in the expression (2b) for V,(r). 
Thus the second approximation is given by sub- 
stituting (3) (where R; is regarded as known) in 
(2b), and then solving Eqs. (2a). Excluding the 
case of all /’s being zero, this makes V; a function 
of r;, 0;, and the /’s and m’s of the different elec- 
trons, but independent of ¢;. Eq. (2a) then has 
solutions of the form e"*’f,,(7;, 0;), where m is an 
arbitrary integer. The eigenvalues ¢; are non- 
degenerate, and it is clear that we must take one 
of the solutions with m=m;.' Thus the second 


approximation is 


¥;(rj) =e'™s?i fm, (r;, 03) (G=1,---,N) (4) 
and proceeding successively in this way, we see 
that the final solution is also of the form (4). We 
indicate by the superscript m that f,,, depends on 
the m’s of all electrons, and not merely on m; 
itself. We may expect that fnm,”(r;, 0;) will not 
differ very much, however, from the initial 
central field approximation P,,,1,(8);R;j(7;). 

But the solution (4) makes not only the eigen- 
functions, but also the eigenvalues «; and the 
energy of the whole atom a function of the mag- 
netic quantum numbers m (as well as of the n’s 
and /’s). This result is clearly incorrect, and, in 
fact, renders the solution indeterminate, since m. 
may have any integral value between +/;. To 
obtain a result for the energy which is independ- 
ent of the m’s, it is therefore necessary to average 
over the m’s. But the assumption (3) (with the 
R’s assumed independent of the m’s) is com- 
patible with the equations (2) if the potentials 
V, are averaged in this way; for each term in the 
summation in (2b) is then of the form 


lk 
1/2 f (1/ru)(2h +1) 1 7 [Pmitk(Ox) ]? 


ms tk 


. [Ri(rz) }*dr 


1 Which eigenvalue is to be taken could be fixed by 
consideration of the initial approximation (thus intro- 
ducing the principal quantum numbers of the different 
electrons). We might, e.g., employ the method given by 
the author (Trans. Roy. Soc. Canada 30, 75 (1936)) for 
solution of (2a) (with V; known) by successive approxi- 
mations, starting with (3) as initial approximation. 


and the summation with respect to m, has just the 
value (2/,+1)/2,? so that we have 


LV,(t;) Jav= — 2/1; 


+1/4er > [Ri(re) ]?/rje-drx, 


k=1(k2+)9) 


which is a function of r; only, so that the assump- 
tion (3) is permissible, and we obtain, in fact, just 
the Hartree equations for the radial functions. It 
may be noted that such an averaging over mag- 
netic quantum numbers is necessary even when 
the restriction (3) is made before applying the 
variational principle, in order to obtain the 
Hartree equations.* It is true that we might 
perform the averaging after solving Eqs. (2) (in 
the form (4)) rather than before; but it is very 
doubtful if this would improve the result, and the 
process of solution would, of course, be much 
more laborious. 

The dependence of the solution on the mag- 
netic quantum numbers is, of course, due to the 
spatial degeneracy of the central field solution 
(3), and we might try to correct for this, while 
still retaining essentially the s. c. f. method, by 
assuming for the wave function a linear com- 
bination of products of the form (cf. (4)) 

” 
IT e'i*igm;(r;, 0;) (5) 


7=1 


for the different possible values of the m’s, i.e., 
for m;= —l;, ---,1;. The functions g,,, (not neces- 
sarily identical with the f,, considered above) 
would then be determined by the variational 
principle. In this case, we must certainly impose 
the condition that the wave functions belonging 
to a given state of the whole atom must have the 
correct transformation properties with respect to 
change of axes.‘ We may first suppose the gm, in 
(5) expanded in the form 


gm;(7;, 0;) = > A mp;(1j) Pmjp;(9;). 


pj=|mj] 


Now it is impossible to find linear combinations 


2See, e.g., Sommerfeld, Wellenmechanischer Ergén- 


zungsband (1929), 103. 

3 Condon and Shortley, Theory of Atomic Spectra (1935), 
356. 

4 See, e.g., Wigner, Gruppentheorie, etc. (1931). 





of products of the form 


N 
II e*™ ii Pm sp; 6;) 


7=1 


with the same coefficients (in the linear combina- 
tion) for different values of the p; which transform 
in the correct manner under a rotation of axes; so 
that the only possibility of finding linear com- 
binations of functions of the form (5) which have 
the correct transformation properties and which 
tend continuously to the correct wave functions 
as the electronic interaction be imagined to tend 
to zero, is to take 


gm;(r;, 8;) =const. X R;(r;) Pm; 1;(8;), (6) 


where R, is independent of the m’s. We then have 
wave functions of the form 


. 
vVu"=IT Rr): DL CL, M; my, ---, my) 
j7=1 Myo ++, My 
. 
xI e'i¢i Pm; 1;(0;) (7) 
gi 


belonging to a given L state of the whole atom, 
where M goes from —ZL to L, and in the summa- 
tion m, goes from —/, tol,, etc. The constants C 
are then independent of the electronic inter- 
action, and are, in principle, determinable from 
group-theoretical considerations alone.‘ 

It thus appears that any wave functions which 
are built up as linear combinations of functions 
such as (5) without the ‘“‘central field’’ assumption 
(6) do not have the correct transformation 
properties, and are therefore to be rejected. Ap- 
plication of the variational principle with the 
wave function (7) yields equations for the radial 
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functions R which differ slightly from those of 
Hartree in that the field for each electron is 
different for the different states arising from the 
same electronic configuration. The corrections 
obtained in this way, however, are of about the 
same order of magnitude as the exchange cor- 
rections, so that there is no point in taking them 
into account without at the same time including 
the exchange corrections.® 

In all the above we have neglected exchange 
degeneracy and the Pauli principle; but even if 
we take these into account, the conclusion re- 
garding the necessity for central fields remains 
valid. For the correct wave function may be 
obtained by first solving the wave equation when 
exchange degeneracy and spin are neglected, and 
then forming linear combinations of this solution 
and its various permutations, the coefficients 
being appropriate spin functions. But these spin 
functions are again independent of the electronic 
interaction and are determinable from group- 
theoretical considerations alone. It is therefore 
sufficient to study, in the first place, those ap- 
proximate wave functions which neglect spin and 
exchange degeneracy, and our previous conclu- 
sion remains true. The generalization of (7) to 
include exchange and the Pauli principle then 
yields, when the variational principle is applied, 
precisely the radial equations of the s. c. f. with 
exchange.® 

5 Cf. Slater, Phys. Rev. 35, 210 (1930). The difference 
between the Hartree equations and those considered here 
is due to the fact that in the former the wave function 
assumed in the variational principle does not have the 
proper transformation properties. This is not so, however, 
in the equations of the s. c. f. with exchange (see, e.g., 
Hartree and Hartree, Proc. Roy. Soc. A154, 588 (1936 
Tnis point is discussed more fully for two-electron con 


figurations in a paper (not yet published) of the writer, 
where a generalization of the s. c. f. equations is also given. 
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THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may be secured by 


addressing them to this department. Closing date 


s for this department are, for the first issue of the 


month, the eighteenth of the preceding month; for the second issue, the third of the month. The Board 


of Editors does not hold itself responsible for the 


opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


On the Anomalous Diffraction Gratings. II 


In my former letter on this question! I have pointed out 
that, if one considers the presence of Sommerfeld’s surface 
waves on a grating, one may explain the part (b) of the 
“typical phenomenon” of Wood and Strong. The introduc- 
tion of these waves was physically justified but quite 
arbitrary from the point of view of the mathematical 
formalism. Owing in particular to critical discussion with 
Professors Laue and Sommerfeld, to whom I am very 
indebted, I have now attempted to overcome the mathe- 
matical difficulties. 

Sommerfeld’s method, to find the field of an oscillator 
on a metallic surface, consists in developing this field in a 
Fourier integral of cylindrical waves. A particular trans- 
formation of this integral shows that the superficial part 
of the field is represented by a superficial damped wave. 
Following this method I have represented the complex of 
reflected, refracted and diffracted waves in the theory of 
the metallic gratings as a Fourier integral of plane waves. 
Rayleigh’s theory follows then automatically without a 
naive application of the principle of Huygens. 

It was only on this point that I noted that the bright 
bands, part (5) of the “typical phenomenon,” are already 
given by developing the method of Rayleigh and Voigt. 
These bands are due to the same particularity in the mathe- 
matical formalism as Sommerfeld’s superficial waves; they 
have their center at the same wave-length as that which 
corresponds to the resonance of Sommerfeld’s waves when 
we represent the grooves of the grating as a series of har- 
monic resonators. I may also say that, although I had not 
seen the simple mathematical aspect of the experimental 
phenomenon, the explanation which is given in my former 
work was physically correct. 

In Rayleigh’s work the metal is thought to be perfectly 
conducting; then the intensity of a continuous diffraction 
spectrum, as a function of the wave-length, becomes 
infinite at Rayleigh’s wave-lengths, and so does its first 
derivative. When the conductivity is not infinite this 
singularity breaks down into two parts: (a) the deriva- 
tive remains infinite at Rayleigh’s wave-lengths, (0d) 
the intensity becomes finite and has a wide and displaced 
maximum. 

In conclusion I wish to summarize the exact meaning of 
the results obtained, which seem now to be satisfactory. 
We may calculate the complex amplitudes (functions of 
the wave-length) of the diffracted waves in a continuous 
spectrum as a series of ascending powers of the depth of the 
grooves. One of the terms of second order: (a) is dis- 
continuous at Rayleigh’s wave-lengths, (6) has a wide 
for slightly larger wave- 

? 


maximum—in its absolute value 


lengths. The qualitative and quantitative influence of thes« 
irregularities on the resulting intensity of the spectrum 
depends very largely upon the profile of the grating as 
pointed out in the former letter. 

A complete account will be published in the Annalen der 
Phystk. 

Uco Fano 
Torino, Italy, 


January 1, 1937. 


1 Fano, Phys. Rev. 50, 573 (1936). 


The Infrared Absorption of the Hydrogen Halides in 
Organic Solvents 


The appearance of the interesting papers by West and 
his colleagues on the infrared and Raman spectra of simple 
gases in nonionizing solvents! ? has suggested to the writer 
the propriety of recounting his observations on the infrared 
absorption of hydrogen bromide and hydrogen iodide in 
organic solvents, made in continuance of work on the 
spectra of hydrogen chloride solutions.*: ¢ In the case of the 
aromatic solvents, chlorobenzene, nitrobenzene, and 
m-nitrotoluene, the solubility of HBr and HI is so slight 
that no definite conclusions could be drawn concerning 
the absorption of the solute, although it appears that the 
fundamental frequencies are lower than in the gas. In 
benzene solution the fundamental bands of HBr and HI 
appear at 3.94 and 4.54, respectively, the absorption 
maximum not differing appreciably from that of the free 
gases. The observations on the iodide solutions had to be 
made immediately after the solutions were prepared, since 
the solute decomposes very rapidly. The HBr solutions 
offered similar difficulties, but the decomposition took 
place more slowly. In the case of every solution studied, 
the absorption coefficient of the dissolved gas was greater 
than that of the free gas. The absorption of HCl was found 
to be more intense than that of HBr and HI, a result to be 
expected from the studies made on these compounds in 
the gaseous state. The studies have been extended to 
include solutions of HCl, HBr and HI in ether, in which 
solvent a high degree of ionization is known to occur. In 
each ether solution the bands characteristic of the solute 
were entirely absent, a result in accord with West's 
observations. As yet no work has been done on solutions of 
hydrogen fluoride. 

DupLEY WILLIAMS 

University of Florida, 

Gainesville, Florida, 
January 28, 1937. 
1W. West and P. Arthur, J. Chem. Phys. 5, 10 (1937). 
2W. West and R. T. Edwards, J. Chem. Phys. 5, 14 (1937). 


3 E. K. Plyler and D. Williams, Phys. Rev. 49, 215 (1936) 
*D. Williams, Phys. Rev. 50, 719 (1936) 
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Some Generalizations of the § Transformation Theory 


According to a hypothesis of Pauli, worked out in detail 
by Fermi and others,' the emission of an electron in the 
process of 8 transformation must be accompanied by 
simultaneous emission of a neutrino in order to satisfy 
the conservation of spin and energy (m—p+e*t+v or 
p—n+e-+v). We should like to discuss here the possibility 
of two other similar processes: The emission of a pair of 
electrons (n—-n+e++e~ or p>p+e*+e-) and the emission 
of a patr of neutrinos (n—n+v+v or p->p+r+yr). Such 
processes evidently do not correspond to nuclear trans- 
formations though they could occur together with y-radia- 
tion if the nucleus is excited 

These processes may be of importance for the explanation 
of forces between heavy particles and also their magnetic 
moments. In principle, such an explanation was given? on 
the basis of (e, v) pair emission (ordinary 8 transformation) ; 
however, for both effects the calculated results (accepting 
U-K interaction) were too small by a factor of about 10". 
Considering the interaction forces and magnetic moments 
as due to (e*, e~) or (v, v) pair emission, however, we can 
obtain the correct values, since we have considerable 
freedom in choosing the probabilities for these new proc- 
esses. We shall ascribe the above two effects to (e*, e 
emission because, in order to give a reasonable explanation 
of magnetic moments of heavy particles, the momentary 
emission of pairs possessing sufficiently large magnetic 
momenta is needed. 

In the case of (et, e~) emission (and also (», ») emission 
the charge of the heavy particle remains unchanged, 
although the spin may change. This seems to correspond 
to the fact that forces between heavy particles depend on their 
spins but not on their charges. Indeed, experimental evi- 
dence shows that the (m, p) and (p, p) forces are equal if 
the spins of two particles are antiparallel. The (m, p) force 
for parallel spins is an attraction and we conclude that the 
same is true for the (p, p) force. Thus, apart from Coulomb 
forces, attraction between any pair of heavy particles results. 
If we furthermore assume that the total potential energy 
of a nucleus is a sum of interaction energies of separate 
pairs, then, owing to the absence of valency saturation for 

our type of exchange-phenomena, heavy nuclei would 
collapse to a radius comparable with the range of forces 
between constitutional particles.‘ Therefore, we must 
assume that the interaction between two particles must be 
affected (reduced) by the presence of other particles. lf several 
pairs may be exchanged simultaneously between two or several 
heavy particles,® such an effect can be easily understood. 
It has been shown by perturbation calculations that a 
perturbing potential which leads to a correct magnitude 
of interaction gives rise to high order perturbations which 
converge rather slowly;* from this it follows that the above- 
mentioned simultaneous exchange of several pairs may 
possess a considerable probability. Assuming that the 
probability of (e~, e+) transformation is 10" times larger 
than the probability of the ordinary (e, v) transformation, 
attractive forces and magnetic moments can be described 
quantitatively. This assumption would also lead to (e~, e*) 
pair emission from excited nuclei, which process would 
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compete with the y-ray emission whenever the excitation 


energy exceeds 2 mc*. If the (e~, e*) pair is emitted with 


one mv total kinetic energy, the probability would be 10° 
sec.-! (10% 8 decay constant for this energy) which, for 
heavy nuclei, is comparable with the probability of the 
Dirac pair formation by y-rays in the vicinity of th 
emitting nucleus.’ 

It seems reasonable to expect that the ratio of the proba 
bilities of (v, ») and (», e) emissions is of the same order of 


magnitude as for (», e) and (e*, e~) emission. This would 


lead to comparable probabilities for (vy, v) pair emission 
and for the emission of quanta of gravitational radiation.* 
G. GamMow 
E. TELLER 
George Washington University (GG and ET), 
Carnegie Institution of Washington (GG), 


Washington, D ° 
January 19, 1937. 


1E. Fermi, Zeits. f. Physik 88, 161 (1934); G. Uhlenbeck and K 
Konopinski, Phys. Rev. 48, 7 (1935 

2? D. Iwanenko and I. Tamm, Nature 133, 981 (1934); C. Wick, Rend 
R. Nat. Acad. 21, 170 (1935); C. Weisziicker, Zeits. f. Physik 102, 572 
1936). 

M. A. Tuve, L. R. Hafstad, and N. P. Heydenburg, Phys. Rev. 50, 

806 (1936); G. Breit, E. Condon, and R. Present, Phys. Rev. 50, 825 
1936). 

*W. Heisenberg, Rapports de Congres Solvay, 1933 

§D. Iwanenko and A. Sokolow, Nature 138, 684 (1936); L. and G 
Nordheim, Phys. Rev. (in press 

*H. Bethe and R. Bacher, Rev. Mod. Phys. 8, 82 (1936 

7 R. Oppenheimer and M. Plessett, Phys. Rev. 44, 53 (1933) 

8 A. Einstein, Berlin, Sitzber. Akad. Wiss. 688 (1918). 





Isotopic Constitution of Neodymium 


In a previous paper an example of a mass spectrum ob- 
tained from a mixture of rare earth elements was given.! 
This included masses at 148 and 150 which could not be 
identified with any known isotopes. Dr. Aston has sug 
gested? that they are new isotopes of neodymium. I have 
recently analyzed the ions from a spark between fairly pure 
neodymium electrodes, and find that the masses at 148 and 
150 belong to this element (Fig. 1). The faint mass at 141 


Nd 142 144 146 148 150 
| | | | | 
; 


Fic. 1. Mass spectrum of neodymium 


— 


is probably due to a trace of praseodymium. No other 
impurities were found. 

Samarium ions were obtained from sparks to an elec- 
trode made by filling a nickel tube with a mixture of 
samarium oxide and aluminum. The isotopic constitution 
was found to be the same as that observed by Dr. Aston 
with the abnormal intensity relationships among the 
isotopes reported by him. 


A. J. DEMPSTER 
Ryerson Physical Laboratory, 
University of Chicago, 
Chicago, Illinois, 
February 1, 1937. 


1 Proc. Am. Phil. Soc. 75, 735 (1935) 
2 F. W. Aston, Nature 137, 613 (1936). 
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Nebulae as Gravitational Lenses 


Einstein recently published! some calculations concern- 
ing a suggestion made by R. W. Mandl, namely, that a star 
B may act as a “gravitational lens” for light coming from 
another star A which lies closely enough on the line of sight 
behind B. As Einstein remarks the chance to observe this 
effect for stars is extremely small. 

Last summer Dr. V. K. Zworykin (to whom the same 
idea had been suggested by Mr. Mandl) mentioned to me 
the possibility of an image formation through the action of 
gravitational fields. As a consequence I made some calcula- 
tions which show that extragalactic nebulae offer a much 
better chance than stars for the observation of gravitational 
lens effects. 

In the first place some of the massive and more concen 
trated nebulae may be expected to deflect light by as much 
as half a minute of arc. In the second place nebulae, in 
contradistinction to stars, possess apparent dimensions 
which are resolvable to very great distances. 

Suppose that a distant globular nebula A whose diam- 
eter is 2£ lies at a distance, a, which is great compared with 
the distance D of a nearby nebula B which lies exactly 
in front of A. The image of A under these circumstances 
is a luminous ring whose average apparent radius is 
B=(yoro/D)*, where yo is the angle of deflection for light 
passing at a distance ro from B. The apparent width of the 
ring is A8=t/a. The apparent total brightness of this 
luminous ring is g times greater than the brightness of the 
direct image of A. In our special case g=2/a/ED, with 
1=(yoroD)}. In actual cases the factor g may be as high as 
q= 100, corresponding to an increase in brightness of five 
magnitudes. The surface brightness remains, of course, 
unchanged. 

The discovery of images of nebulae which are formed 
through the gravitational fields of nearby nebulae would 
be of considerable interest for a number of reasons. 

(1) It would furnish an additional test for the general 
theory of relativity. 

(2) It would enable us to see nebulae at distances greater 
than those ordinarily reached by even the greatest tele- 
scopes. Any such extension of the known parts of the uni- 
verse promises to throw very welcome new light on a 
number of cosmological problems. 

3) The problem of determining nebular masses at 
present has arrived at a stalemate. The mass of an average 
nebula until recently was thought to be of the order of 
My =10° Mo, where Mo is the mass of the sun. This esti- 
mate is based on certain deductions drawn from data on 
the intrinsic brightness of nebulae as well as their spectro- 
graphic rotations. Some time ago, however, I showed? 
that a straightforward application of the virial theorem to 
the great cluster of nebulae in Coma leads to an average 
nebular mass four hundred times greater than the one 
mentioned, that is, My’=4xX10"Mo. This result has 
recently been verified by an investigation of the Virgo 
cluster.2 Observations on the deflection of light around 
nebulae may provide the most direct determination of 
nebular masses and clear up the above-mentioned 


discrepancy. 
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A detailed account of the problems sketched here will 
appear in Helvetica Physica Acta. 
F. ZwIckKy 
Norman Bridge Laboratory, 
California Instit of Technology, 


Pasadena, California, 
January 14, 1937. 





\. Einstein, Science 84, 506 (1936 
2 F. Zwicky, Helv. Phys. Acta 6, 124 (1933 
Sinclair Smith, Astrophys. J. 83, 23 (1936 


Emergence of Low Energy Protons from Nuclei 


In some experiments recently described! the emission of 
protons in alpha-particle induced transmutations has been 
studied. In several cases the interesting fact was noticed 
that protons of relatively low energy were emitted in 
considerable numbers. Thus for each of the reactions 


Al? +Het—Si* + H!, 
Ps! +Het—+S* +H}, 
C5 +Het+A*® +H}, 
Ca*®?+ He*—Sc*+ H}, 


a group of protons of maximum range 20 cm or less is found 
and the yield is in general large (more than one-third of the 
total number of protons emitted). In each case protons of 
range 10 cm are observed with no apparent diminution of 
the probability of emission. The question arises as to how 
these low energy protons get out of the composite nucleus. 

In recent experiments in this laboratory the excitation 
curve for the emission of neutrons from argon under alpha- 
particle bombardment has been plotted and the nuclear 
radius found to be 7.3 10~" cm which is in accord with 
Bethe’s revised radii for the radioactive elements? and may 
be taken as a basis for calculation of the nuclear radii of 
Si**, S*, A*8, Ca*® and Sc*. Other evidence (e.g., scattering 
experiments) indicates, if anything, smaller radii than those 
found in this way. In Table I are given the radii so calcu- 
lated, together with the heights of the corresponding proton 
barriers and the range of a proton just able to surmount 
them. It will be seen that in every case the experimentally 
observed ranges are smaller than necessary to scale the 
barrier. It therefore appears that we can draw one of two 
significant conclusions from the experimental data. Either 
barriers to emerging protons are abnormally low or the com- 
posite nucleus containing the final product element and the 
proton has a finite lifetime sufficiently long to enable the 
proton to leak through the barrier. The latter view, which is 
in accordance with Bohr’s conception of transmutation,* 


TABLE I. 

PROTON RANGI EXPERI 

NUCLEAR BARRIER TO SCALE MENTALLY 

PRODUCT RADIUS HEIGHT BARRIER FOUND 
NUCLEUS <x 10cm) Mev) cm) RANGE 
Si” 6.7 3.0 14.0 <10 
S* 6.9 3.3 16.5 <10 
A® 7.2 3.6 19.0 <10 
Ca® 7.4 3.9 22.0 14 
Sc* 7.5 4.0 23.0 <10 
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seems to give the simplest explanation consistent with the 


remainder of our knowledge about nuclei. It is of some 
interest to estimate the lifetime of the composite Ti* 
nucleus formed in the last reaction. This is found to have a 
minimum value of approximately 10~?° second, but it is 
probably much longer as it may be that protons of much 
smaller range than have been detected are actually emitted. 


ERNEST POLLARD 
Sloane Physics Laboratory, 
Yale University, 
New Haven, Connecticut, 
February 2, 1937. 

1E. Pollard and C. J. Brasefield, Phys. Rev. 50, 890 (1936); 51, 8 
(1937). 

2H. Bethe, Phys. Rev. 50, 977 (1936 

3 The idea is developed in some detail by W. D. Harkins, Phys. Rev 
51, 52 (1937). 


Cyclotron Operation without Filaments 


With the use of a transmission line! to supply high fre- 
quency power to the dees of a cyclotron, it is possible to 
operate with a continuous glow discharge in the vacuum 
chamber. Under these conditions it is unnecessary to use 
filaments to create ions for acceleration since the high 
frequency discharge supplies them. In this way we have 
obtained a steady ion beam current of 5 microamperes. 
The operating conditions are very stable and the maximum 
beam current may be drawn for any desired length of time. 

Fig. 1 shows a typical curve relating pressure in the 
vacuum chamber and ion current drawn from the cyclotron 
under the above conditions. The end of the curve shows the 
pressure at which the discharge ceased for that particular 
run. The optimum pressure is of the order of 10° mm Hg. 
At lower pressures it is necessary to use filaments to obtain 
a beam from the cyclotron. An increase in pressure from 
about 10-7 mm Hg up to the above point increases the 
beam current when using filaments. With filaments the 
maximum steady beam we have been able to get is 8 
microamperes. 

When operating without filaments, the pressure is in- 
creased to about 10° mm Hg where the h.f. potential 
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Fic. 1. Relation between gas pressure and ion beam current. 
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strikes the discharge. The excess gas is then pumped out 
until the beam is a maximum and steady operating condi- 
tions are-obtained. The beam current under these condi 
tions, as well as when using filaments, increases with an 
increase in power transmitted to the dees. 

The operating conditions with a continuous glow dis 
charge in the vacuum chamber are much more stable than 
with no glow discharge, since a small change in gas pressure 
does not affect the tuning of the dee circuit appreciably. 
When bursts of gas occur the line impedance match is 
destroyed and the power delivered to the dees decreases 
for a short time, but the pumps soon reestablish the proper 
pressure conditions whereupon the power delivered to the 
dees and the beam current go up to normal again. 

P. GERALD KRUGER 
G. K. GREEN 
F. W. STALLMANN 
Department of Physics, 
University of Illinois 


Urbana, Illinois, 
January 27, 1937 


1Green and Kruger, Phys. Rev. 51, 57 (1937 


Asymmetric Zeeman Effect Patterns in a 
Complex Spectrum 


Zeeman effect analysis of an atomic spectrum usually 
shows each energy level to be split symmetrically about its 
field-free value. Asymmetry has hitherto been recognized 
neglecting hyperfine structure) only in LS coupling spectra 
and certain especially simple two-electron configurations." 
In these special cases the wave number contribution 04 Mg 
to a sublevel due to the external field can be made to attain 
the same order of magnitude as that due to spin-orbit 
interaction if the latter is sufficiently small (Paschen-Back 
effect). Here g is the weak field Zeeman effect splitting 
factor, M is the magnetic quantum number, and oyc is the 
Larmor precession frequency. 

This limitation of the known asymmetrically split levels 
to simple spectra, is due to the accidental lack of known 
level groups sufficiently close to show the interaction, in 
the more complex spectra that have been magnetically 
analyzed. In addition to the well-known symmetric 
splitting term oy Mg, each of the sublevels with the same 
M value, of any two neighboring levels of the same parity 
whose J values differ by not more than one, may be ex- 
pected to have in the expression for its wave number, a 
repulsion term with the resonance denominator D=vg—v» 
+0nM(ga—gv). Here ga, gp are the weak field g values and 
hcva, hev, the field-free energies of the levels. This term 


equals 
Kox?M[D if Ja-Js =0; 
Koy?(22—M?)/D if |Ja—Jo| =1; 
0 if |Ja—J»| >1. 


Here K is a constant characteristic of the pair of levels, 
and in the case |J,—J,| =1, J is the larger of the two 
values J, Jy. The expressions are approximate. Levels with 
|Ja—Js| >1 may be expected to show some repulsion, 
associated with the fact that J is not a good quantum 
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number in a magnetic field; but if the repulsion is suffi- 
ciently small it may be hard to detect, and experiment 
may show sublevels with the same M crossing. « 

The spectrum of neutral tungsten? is rich in levels about 
which nothing is known except parity, v, J, g, and transition 
intensities. Lacking a priori knowledge of the value of K, 
we have used the magnitude of oy?/D in all cases where 

Ja—J»| =0 or 1 as a measure of the expected nonlinear 
contribution to the energy of a sublevel.* Although there 
are several pairs a good deal closer, the pair 21,448.662(1.3), 
21,453.76,(2.56)—the data are given in the form v,(g) for 
each level—shows by far the largest expected nonlinear 
contribution within the range of field strengths ordinarily 
attainable for Zeeman effect studies. For the M=—1 sub- 
levels of this pair, o4*/D equals 0.26 cm™ when oy=1 
cm (i.e., in an external field of 21.4-10* gauss) and 1.6 
cm in twice as great a field. For og =4 cm™, D would be 
zero, i.e., the sublevels would cross if the repulsion term 
were absent. The level 21,453 happens to be one of the very 
few in the whole spectrum to which one may at present 


assign some partly meaningful quantum numbers other 


than J. It acts in some respects as though it deserved the 
name 5d‘6s6p *P;. “‘Multiplet separations” in tungsten are 
of the order of thousands of units. It is not easy to see how 
21,448 and 21,453 could be very closely related dynamically. 

Of the 657 tungsten lines whose Zeeman effect patterns 
are listed in the highly consistent work of Jack,‘ many 
lines are listed with asymmetries of intensity, but only 
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seven with asymmetries of position. Three of the asym- 
metric pattern measurements may be explained easily as 
due to faintness or overlapping of complicated patterns; 
one belongs to a doubly classified line; one is still unclassi- 
fied. The remaining two (v22,003 and v18,128) belong, 
respectively, to the levels 21,448 and 21,453. Both patterns 
are incomplete. They appear to afford definite evidence of 
interaction, with asymmetries on the order of 0.1 cm™ in a 
field for which 04 =1.072 cm™. Of course this has spoiled 
the g value determinations, so the values 1.3, 2.56 listed 
above are not the true weak field g values. Quantitative 
analysis has been postponed pending the availability of 
better patterns. An attempt is being made to obtain such 
patterns. 

We are indebted to Professor Breit for a discussion of the 


problem. 
J. E. Mack 
University of Wisconsin, 
Madison, Wisconsin. 
O. LApPortT! 


University of Michigan, 
Ann Arbor, Mich 
January 22, 1937. 
1 For a summary see Condon and Shortley, Theory of Atomic Spectra, 
§5, reference 16. 

20. Laporte and J. E. Mack, forthcoming publication. 

3 Equally justifiable would have been the procedure of using og2M/D 
where Jg—J»p=0 and oy*(J?—M?2)/D where Jeg—Jp=1. This would 
have favored the pairs 37,297.56, 37,309.13s and 40,463.16, 40,476.37 
Actually the (unresolved) patterns of the lines from these levels are 
listed as being symmetric. In these cases at least, K must be of the order 
of 1/J? or smaller. 
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R. Jack, diss. Géttingen 1908; Ann 1909). 
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